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(57) ABSTRACT 

In a coded orthogonal frequency division multiplex 
(COFDM) system n-bit data words are encoded as 
2 m -symbol code words (binary, quaternary, octary, etc.). The 
code words are selected for desired low peak-to -me an 
envelope power ratio (PMEPR) characteristics of transmis- 
sions over a COFDM channel, from a set of cose ts of a linear 
sub-code of a code having a specified generator matrix. The 
code words thus identified by the procedure described can, 
even for values of m in excess of 3, simultaneously limit the 
PMEPR to 3 dB, provide specified error control 
characteristics, be implemented in a feasible manner using 
analytical circuit techniques (e.g. with combinatorial logic), 
and include sufficiently many different code words to enable 
data to be transferred at useful rates. Other selections of code 
words can be made, enabling a higher maximum PMEPR or 
a reduced error detection capability to be accepted in order 
to obtain a higher code rate. 

34 Claims, 10 Drawing Sheets 
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METHODS AND APPARATUS FOR 
ENCODING AND DECODING DATA 

TECHNICAL FIELD 

This invention relates to methods and apparatus for 
encoding and decoding data, and to communications sys- 
tems incorporating such methods and apparatus; in 
particular, though not exclusively, the invention relates to 
such methods, apparatus and systems for transmitting data 
using multi -carrier frequency division multiplexing, such as 
coded orthogonal frequency division multiplexing 
(COFDM). 

BACKGROUND ART 

The increasingly widespread use of powerful computer- 
based equipment and appliances is resulting in rapidly 
increasing requirements for high-speed, high-capacity 
communications, based on both wireline and wireless tech- 
niques. One way of accommodating this increase is to design 
communications systems for operation at higher and higher 
symbol rates. However, an increase in symbol rate typically 
involves a corresponding decrease in internsyrabol spacing, 
and this in turn can lead to greater inter-symbol interference 
(ISI) and consequential errors in symbols as received. In the 
case of wireline links ISI can arise, for example, from signal 
reflections at imperfectly matched interfaces between dif- 
ferent portions of the communications link. For wireless 
links a significant source of ISI is multi-path propagation 
involving signal reflections from objects (such as buildings) 
in the space traversed by the signals; mobile wireless 
communications can be especially susceptible to ISI from 
this source, because movement of a mobile unit using the 
link can cause continuous and complex variation in the 
multiple propagation paths extending to the mobile unit. 

A technique which is being actively studied and devel- 
oped to increase communications system capacity and coun- 
teract ISI is multi-carrier (or multitone) frequency division 
multiplex operation, such as COFDM operation. In a 
COFDM system multiple symbols are transmitted 
simultaneously, each symbol controlling the modulation 
(e.g. by phase-shift keying) of a respective one of multiple 
carrier signals having closely-spaced frequencies. By choos- 
ing the carrier-frequency spacing to have an appropriate 
relationship to the rate at which symbols are transmitted on 
each carrier, it is possible to ensure that each individual 
carrier can be demodulated without interference from the 
carriers adjacent to it (in which case the carriers are said to 
be orthogonal to one another). The set of symbols transmit- 
ted simultaneously on respective carriers during a symbol 
period (a code word) is chosen to encode a group of data 
symbols (a data word) and includes redundancy (i.e. the 
code word has more symbols than the data word); this 
redundancy is used upon reception for detection of errors in 
the symbols as received, and, if desired, for correction of 
those errors. 

COFDM systems counteract ISI in part because the 
transmission of multiple symbols in parallel enables a 
desired data rate to be attained with a lower symbol rate than 
would be necessary if a serial system were used to transmit 
individual symbols in succession. A lower symbol rate 
implies a longer duration for each symbol and thus a longer 
period between symbol transitions, which reduces the 
impact of ISI. The error detection and correction capabilities 
provided by the inclusion of redundancy enhance the resis- 
tance to data corruption, by enabling some errors which 
occur in received symbols owing to signal fading and other 
causes to be detected and corrected. 



'3,859 Bl 

2 

However, present COFDM systems encounter difficulties 
in obtaining the maximum potential benefit for a particular 
transmitted signal power, owing to limitations arising from 
the typical peak-to-mean envelope power ratio (PMEPR) of 

5 the actual transmitted signal. Because this signal is effec- 
tively the sum of several signals at closely spaced frequen- 
cies and with respective phases which change at each 
symbol transition, the overall signal tends to exhibit tran- 
sient pronounced peaks separated by substantially longer 

10 intervals of much lower though still varying amplitude. The 
transmitting equipment must be adjusted to reproduce the 
peaks without clipping or other distortion, so for much of the 
time, in the intervals between peaks, that equipment is 
operating at a signal level well below its maximum capa- 

15 bility. As a result the geographical range of the equipment is 
significantly lower than would be expected for its nominal 
power rating, or alternatively a much more powerful equip- 
ment must be provided to cover the desired range and then 
operated inefficiently. 

20 It is known that this problem can in principle be alleviated 
by careful selection of the code words which encode the 
possible data words (for example in 'Block coding scheme 
for reduction of peak to mean envelope power ratio of 
multicarrier transmission schemes* by A. E. Jones, T. A 

25 Wilkinson and S. K. Barton, Electronics Letters, Dec. 8, 
1994, vol.30, no.25, pp. 2098-2099). However, in selecting 
code words to reduce the PMEPR it is important also to 
ensure the code has good error correction properties. 
Furthermore, implementation of a high-throughput commu- 

30 nications system desirably involves the use of long code 
words and of rapid encoding and decoding procedures; the 
use of simple look-up tables, though readily implemented, 
does not lend itself to rapid operation using compact (and 
low-power, inexpensive) circuitry, especially during 

35 decoding, if the number of valid code words to be checked 
is large, as is typically the case with long code words (e.g. 
more than about 16 symbols). It is therefore preferable for 
encoding and decoding operations to be definable in terms of 
some analytical procedure which for example performs 

40 combinatorial logic on the data word or code word symbols 
to convert one to the other, or at least restricts the size of any 
look-up table which is used. 

An additional difficulty is presented by the sheer quantity 
of possible code words from which a selection may be made 

45 to define a code for use in a typical system. As the number 
of symbols in a code word increases the number of possible 
code words increases much faster, even for binary symbols; 
for symbols having just three (ternary) or four (quaternary) 
possible values, the number of possible codes available 

50 becomes so huge even for quite short code words that a 
simple exhaustive search through all possible code words is 
not feasible in a realistic time even just to minimize the 
PMEPR. Attempting to satisfy error correction and imple- 
mentational requirements as well merely adds to the diffi- 

55 culty of such a search, especially as it will not even be 
known whether there actually exists a code word set which 
is of the desired size and which simultaneously satisfies 
preferred criteria relating to PMEPR, error detection or 
correction and ease of encoding and decoding. This difli- 

60 culty is further exacerbated by the desirability of having a 
large number of code words available for use in the chosen 
code, in order to be able to encode a large number of 
different data words (i.e. data words containing relatively 
many symbols) and thus attain an acceptably high data 

65 transfer rate. 

It is known that certain kinds of codes may have some 
desirable properties. Thus the paper * Synthesis of Power 
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Efficient Multitone Signals with Flat Amplitude Spectrum 1 
by B. M. Popovic, IEEE Transactions on Communications, 
vol.39, no.7, July 1991, pp. 1031-1033 indicates that any 
binary or polyphase complementary sequence can be used to 
construct a multitone signal with a crest factor (analogous to 5 
PMEPR) less than or equal to 6 dBV. Examples of these 
sequences have been described, as in the case of binary 
Golay complementary sequence pairs ('Complementary 
series' by M. J. E. Golay, IRE Transactions on Information io 
Theory, vol.IT-7, April 1961, pp. 82-87). However, Popovi 
c's paper is directed to situations such as multipath fading 
measurements and multitone jamming, in which a single 
code word can be used continuously in generating the 
required signal. Accordingly Popovic gives no guidance on 
how multiple such sequences could be identified just to 
attain the 6 dBV limit on PMEPR in a practical code for 
conveying information; and no mention whatsoever is made 
of the other requirements mentioned above for a code 20 
suitable for practical implementation. Indeed Popovic states 
that "the construction of multitone signals having minimum 
peak-to-peak amplitude is an old problem that still has no 
analytical solution". ^ 

It should be noted that Golay complementary sequence 
pairs and codes defined using them are different from and 
should not be confused with so-called Golay codes; Golay 
complementary sequence pairs and Golay codes were first 
defined, independently of one another, by the same 30 
researcher and for that reason alone both are referred to by 
his name. 

It is an object of this invention to provide methods and 
apparatus for encoding and decoding in which a substantial 35 
number of different code words are identifiable in practice, 
despite the potentially huge search space of possible code 
words, and which provide a predictable low limit on 
PMEPR, have specified error detection and correction prop- 4Q 
erties and are feasible for implementation in high-speed 
circuitry. 
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where z has integer values from 0 to 2 m -l, y has integer 
!5 values from m to 0, x indicates bitwise multiplication of 
factors expressed as (m+l)-bit binary numbers, and the 
result of division is expressed as a single digit. 

According to another aspect of this invention there is 
provided a method for transmitting data words as 2 m -symbol 
code words using multi-carrier frequency division 
multiplexing, such as coded orthogonal frequency division 
multiplexing (COFDM), comprising the steps of: 

encoding each data word in accordance with its value and 
code words selected in complete cosets, for desired 
peak-to-mean envelope power ratio characteristics of 
transmission, from a set of one or more cosets of a 
linear sub-code of a code having a generator matrix as 
follows: 



((2.0+ l)x2 m J/2 m 
((2.0+ l)x2 m_t )/2 m " 1 

{(2.0+ l)x20/2* 

((ZO+lJx^/a 1 
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((2.1 + l)x2' n )/2 ff, 
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((2.1 + l)x2°)/2° 
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According to one aspect of this invention there is pro- 
vided a method of encoding data words as 2 m -symbol code 
words for transmission using multi-carrier frequency divi- 
sion multiplexing, such as coded orthogonal frequency 
division multiplexing (COFDM), wherein data words are 50 
encoded in accordance with their value and code words 
selected in complete cosets, for desired peak-to-mean enve- 
lope power ratio characteristics of said transmission, from a 
set of one or more cosets of a linear sub-code of a code 
having a generator matrix as follows: 



((2^ + l)x2 m )/2 m 
((2^ + l)x2 M - 1 )/2 m - 1 

((l« + l)x2')/2' 
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((2.(2 M -l)+l)x2 m )/2 m 
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where z has integer values from 0 to 2 m -l, y has integer 
values from m to 0, x indicates bitwise multiplication of 
55 factors expressed as (m+l)-bit binary numbers, and the 
result of division is expressed as a single digit; and 

decoding each code word using an analytical decoder. 
According to a further aspect of this invention there is 
provided a method of encoding data words as 2 m -symbol 
code words for transmission using multi-carrier frequency 
division multiplexing, such as coded orthogonal frequency 
division multiplexing (COFDM), wherein m is at least 4 and 
data words are encoded in accordance with their value and 
code words comprising a set of one or more complete cosets 
of a linear sub-code of a code having a generator matrix as 
follows: 
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(do+Dxrvr* 
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((Z2 + l)x2")/2" ■•■ ((2.(2"-l)+l)x2 ,,, )/2 m 
((Z 2 + l)x2"- l )/2 w - 1 ■•■ ((Z(2 m -L) + l)x2 OT - , )/2 m - 1 

((2.t + l)x2>)/2> •■• ((2.(2"- 1)+ l)x2')/2> 

((Zz + 1) x 2 1 ) / 2 1 • ■ • ((2.(2" - 1)+ l)x2 l )/2 1 
((Zz + 1) x 2°) /2° • ■ • ((2.(2" - 1) + I) x 2°J/2° 

where z has integer values from 0 to 2 m -l, y has integer 
values from m to 0, x indicates bitwise multiplication of 
factors expressed as (m+l)-bit binary numbers, and the 
result of division is expressed as a single digit. 

According to a yet further aspect of this invention there is 
provided a method of encoding data words as 2 m -symbol 
code words, each code word symbol having 2 J different 
possible values, wherein m is at least 4 and j is at least 3, and 
data words are encoded in accordance with their value and 
code words comprising a set of one or more complete cosets 
of a code having a generator matrix as follows: 

((Z0 + l)xr)/2" ((2.1 + l)x2 m )/2" 

((Z0 + l)x2 m_i )/2 m " 1 ((Zl + l)x2 m_1 )/2 m ~ 1 ••• 

((Z0 + l)x2 J ')/2 jr <(Zl + l)x2*)/2' 

((Z0 + l)x2 1 )/2 1 ((Zl + l)x2 1 )/2 l 
((Z0 + l)x2°)/2° ((Zl + l)x2°J/2 <l 



((Zz + l)x2 M )/2'" ■•■ ((Z(2 w, -l) + l)x2 ,n )/2 m 

((2^ + l)x2 M - 1 )/2 ffl - 1 ■•■ ((2.(2 m -l)+l)x2 m - 1 )/2 m - 1 

15 

((Z* + l)x2>)/2> ••• ((2.(2" - 1) + l)x2- y )/2 J ' 

((Zz+l)x2')/2 l ••■ ((Z(2 m -l) + l)x2 1 )/2 1 
20 ((Zz + l)x2°)/2° - ((Z(2"-l) + l)x2°)/2 0 

where z has integer values from 0 to 2 m -l, y has integer 
values from m to 0, x indicates bitwise multiplication of 
factors expressed as (m+l)-bit binary numbers, and the 
25 result of division is expressed as a single digit. 

According to a farther aspect of this invention there is 
provided a method of identifying a 2 m -symbol code word 
within a set of code words having a generator matrix 

30 ((Z0+ l)x2 m )/2" ((Zl + l)x2")/2" 

((Z0+l)x2 JW - 1 )/2' n - 1 ((Zl + l}x2 M - 1 )/2'"- 1 ■•■ 

((2.Q+l)x2 y )/2 y ((Z1 + 1JX20/2' 

35 

((ZO+Dx^)/^ ((Zl + l)x2 1 )/2 1 
((Z0+l)x2°)/2 t) ((Zl + l)x2°)/2° 

40 

((Zz+ l)x2 m )/2" ((Z(2**-l)+l)x2")/2" 
((Z 2 + l)x2" , - l )/2 m - 1 ((2.(2 m -l) + l)x2' n - 1 )/2' n - 1 



((Zi + l)x2")/2" ■•■ ((Z(2 ffl -l)+l)x2 m )/2 m 45 
((Zz + l)x2 m - l )/2 m_l ■•■ ((Z(2 w -l) + l)x2 B - l )/2 w " 1 

((Zz + 1) x 2') / V ■ ■ ■ ((2.(2" - 1) + L ) x V) IV 

50 

((Zz + 1) x 2 1 ) / 2 1 • ■ ■ ((2.(2" - 1 ) + 1 ) x 2 l ) /2 1 
((Zz + 1) x 2°) /2° • ■ • ((2.(2" - 1 ) + 1) x 2°)/2° 

where z has integer values from 0 to 2 m -l, y has integer 55 
values from m to 0, x indicates bitwise multiplication of 
factors expressed as (m+l)-bit binary numbers, and the 
result of division is expressed as a single digit. 

According to another aspect of this invention there is 
provided a method of encoding data words as 2 m -symbol 60 
code words for transfer over a communication path, each 
code word symbol having 2/ different possible values, 
wherein j is a variable positive integer, and data words are 
encoded in accordance with their value and a code having, 
for any value of j, code words contained within a set of one 65 
or more cosets of a linear sub-code of a code having a 
generator matrix as follows: 



((2.z + l)x2')/2 y ((2.(2" - 1) + 1) x2>)/2' 

((l£ + l)x2')/2 l ••■ ((2.(2" - 1) + DX2M/2 1 
((lz + l)x2°)/2° •-■ ((Z(2"-l) + l)x2°)/2° 

where z has integer values from 0 to 2 m -l, y has integer 
values from m to 0, x indicates bitwise multiplication of 
factors expressed as (m+l)-bit binary numbers, and the 
result of division is expressed as a single digit, and each code 
word symbol has 2/ different possible values where j>l, 
wherein coefficients for rows of said matrix are derived by: 

applying a plurality of fast Hadamard transforms or 
equivalents thereof to a 2 m -symbol input word; and 

combining results of said transforms or their equivalents 
to identify said code word. 

According to another aspect of this invention there is 
provided a coder for encoding data words as 2 m -symbol 
code words for transmission using multi-carrier frequency 
division multiplexing, such as coded orthogonal frequency 
division multiplexing (COFDM), comprising means for 
receiving data words and for encoding them in accordance 
with their value and code words selected in complete cosets, 
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for desired peak-to -mean envelope power ratio characteris- 
tics of said transmission, from a set of one or more cosets of 
a linear sub -code of a code having a generator matrix as 
follows: 

s 

((Z0 + l)x2 m )/2 ,n ((ll + Uxr')/?" 
((10 4-l)x2 m - l )/2 m - 1 ((ll + Dxr -1 )/^ -1 

((2.0 + UX2O/2' ((Zl4-l)x20/2 J ' ■» 10 

(O-O + Dx^)/^ ((Zl^ljx^)/^ 
((10 4-l)x2°)/2° ((2,1 + l)x2°)/2° 

15 

((Zz + l)x2 m )/2" •-• ((2.(2™ -1)+ I)x2" , )/2 B * 
((2.zfl)x2 w -']/2^ 1 ((2.(2 w, -l) + l)x2 m - l )/2^ i 

((2.i4-l)x2J')/2' ■•■ ((2.(2" , -i)+l)x2^/2J' 20 

{(2.1 + 1) x 2 ! ) / 2 1 * ■ ■ ((2.(2 m - 1 ) + 1 ) x 2 l ) /2 1 
((2.z + 1) x2°)/2° • ((2.(2*' - 1) + 1) x 2°)/2° 

25 

where z has integer values from 0 to 2 m -l, y has integer 
values from m to 0, x indicates bitwise multiplication of 
factors expressed as (m+l)-bit binary numbers, and the 
result of division is expressed as a single digit. 3Q 

Circuitry implementing this invention is relatively simple 
and compact. Thus, for example, identifying code words for 
use in encoding does not require exhaustive searches (of 
possibly impractical duration) through enormous numbers 
of code words, even for relatively long code words. 35 
Likewise, decoding does not need enormous look-up tables 
containing all valid code words that may occur, which would 
again be impractical for long code words. The simplicity and 
compactness of practical implementation of the invention, as 
compared to what has previously been believed to be 4Q 
necessary (e.g. exhaustive searches and large look-up tables) 
or at the limits of feasibility, illustrate the significant 
advance which the invention provides. 

BRIEF DESCRIPTION OF DRAWINGS 

45 

Methods and apparatus for encoding and decoding data in 
accordance with this invention in communications systems 
for transmitting data using coded orthogonal frequency 
division multiplexing (COFDM) will now be described, by 
way of example, with reference to the accompanying draw- 50 
ings in which: 

FIG. 1 illustrates multipath propagation between a base 
station and a mobile communications device; 

FIG. 2 is a schematic diagram showing the effect of 
multipath propagation on propagation of a single symbol; 55 

FIG. 3 is a schematic diagram showing intersymbol 
interference arising from multipath propagation; 

FIGS. 4a to 4h are schematic diagrams showing reduction 
in intersymbol interference by the use of multi-carrier fre- 6Q 
quency division multiplexing; 

FIG. 5 shows an encoder and modulator for use in a 
COFDM system; 

FIG. 6 is a block schematic diagram of a COFDM system 
incorporating the encoder and modulator of FIG. 5; 65 

FIG. 7 illustrates the general form of a COFDM signal 
waveform without any limitation of PMEPR; 
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FIG. 8 shows an encoder and decoder incorporating the 
present invention, for encoding 5-bit data words as 8 -bit 
code words; 

FIG. 9 shows an encoder, for encoding 8 bit data words 
as 16 -bit code words; 

FIG. 10 shows a decoder incorporating the present inven- 
tion for decoding code words from the encoder of FIG. 9; 

FIG. 11 shows an alternative decoder which may be used 
for decoding code words from the encoder of FIG. 9; 

FIG. 12 shows an encoder incorporating the present 
invention, for encoding 13-bit data words as 16-symbol 
quaternary code words; 

FIG. 13 shows a decoder incorporating the present inven- 
tion for decoding quaternary code words from the encoder of 
FIG. 12; 

FIG. 14 shows an encoder incorporating the present 
invention, for encoding 18-bit data words as 16-symbol 
octary code words; and 

FIG. 15 shows a decoder incorporating the present inven- 
tion for decoding octary code words from the encoder of 
FIG. 14. 

BEST MODE FOR CARRYING OUT THE 
INVENTION & INDUSTRIAL APPLICABILITY 

FIGS. 1 to 3 illustrate in a schematic manner the problem 
of multipath propagation in the context of mobile commu- 
nications in an urban environment. A base he station 10 
communicates by wireless, e.g. at microwave wavelengths, 
with mobile communications devices such as a mobile 
telephone 12. There is a direct line-of -sight signal path 14 
between the base station 10 and the telephone 12, and, owing 
to the presence of multiple reflectors (for microwave 
signals) in the form of buildings and other structures, there 
are also multiple reflected signal paths 16. These paths are 
typically of differing lengths and thus involve differing 
propagation delays upon arrival at their destination. The 
overall signal as received at, for example, the mobile tele- 
phone 12 is the sum of the various signals arriving at the 
telephone via the different paths. 

The transmitted signal may be considered as comprising 
a succession of symbols modulating a carrier wave (e.g. by 
phase-shift keying). As shown in FIG. 2, a consequence of 
signals traversing multiple paths with different propagation 
delays is that the composite signal 18 as received for a 
symbol is not an accurate replica of the signal as transmitted 
for that symbol. In particular the received composite signal 
is typically longer in duration than the transmitted signal, 
and has a trailing portion with a significant amplitude after 
the point when an accurate replica of the transmitted signal 
would have terminated. If this trailing portion extends for a 
significant fraction of a symbol period, as shown in FIG. 3 
for the pulse nominally terminating at time T, the trailing 
portion will overlap and interfere with reception of the 
following symbol. As a result that following symbol may be 
incorrecdy received, i.e. the value detected for it in the 
receiver may not correctly represent the value as transmitted, 
even in the absence of other noise sources. Signal corruption 
of this kind is known as intersymbol interference (ISI). 

The effects of ISI can be alleviated by extending the 
duration of the transmitted symbols, as shown in FIG. 4a, so 
that each symbol as received extends beyond the trailing 
portion arising from multipath propagation delays. To coun- 
teract the reduction in symbol rate which would ensue, it has 
been proposed to transmit multiple symbols simultaneously, 
each modulating a respective one of multiple carrier waves. 
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Thus if eight carriers are used, as shown in FIGS. 4a to 4A, 
the symbol period on each carrier can be eight times that 
which would be necessary with a single carrier (as shown in 
FIG. 3) while maintaining the same overall symbol rate. If 
the frequency spacing between the multiple carriers is an 5 
integer multiple of the inverse of the symbol period, it is 
possible to recover the individual symbol stream modulating 
each of the carriers without interference from the other 
carriers. This technique is known as orthogonal frequency 
division multiplexing (OFDM). 

In addition to the use of multiple carriers, it has also been 
proposed to provide additional error detection and correction 
capabilities by the use of block encoding. Thus, as shown in 
FIG. 5, a serial stream of data symbols (binary symbols in 
this example) is notionally divided into successive groups or 
words each containing five symbols. A group of five symbols 15 
is input to a five-bit to eight-bit (5B/8B) encoder 30 which 
outputs a corresponding eight-bit code word. This code word 
is applied to a bank of eight phase shift controllers 32 which 
control a bank of eight oscillators 34 producing carrier 
signals at eight successive frequencies at intervals Af apart. 20 
Each individual bit in the code word causes a respective 
phase shift controller to impose either a positive or negative 
phase shift A<|> to modulate the carrier signal produced by the 
respective oscillator, depending on the current value of that 
bit. 25 

The modulated outputs of the oscillators 34 are combined 
in an adder 36 to produce the transmitted signal. In practice 
the functions of the oscillators 34 and the adder 36 are 
conveniently combined in a digital signal processor imple- 
menting an inverse fast Fourier transform (IFFT), as indi- 30 
cated by the dashed box 38. Thus, as shown schematically 
in FIG. 6, data to be transmitted are received in block 40 and 
converted to parallel form in a serial-to-parallel converter 
42. The parallel data are encoded in encoder 30, as described 
above, and the encoded data control generation of the 35 
multi-carrier signal in the modulator/IFFT 32/38. This multi- 
carrier signal is converted to analogue form in a D-A 
converter 44, amplified and transmitted by a transmitter 46. 

Reception is essentially a reverse of this process: the 
signal is received and amplified by a receiver 48, and then 40 
converted to digital form by an A-D converter 50. The digital 
signal is supplied to a demodulator 52 which comprises a 
digital signal processor implementing the fast Fourier trans- 
form (FFT) to separate and demodulate the component 
signals in the composite multi -carrier signal. The demodu- 45 
lated data are passed to a decoder 54 which reverses the 
encoding applied by the encoder 30, and typically also 
performs error detection and correction functions. Thereaf- 
ter a parallel-to-serial converter 56 changes the data into 
serial form to be output at block 58. 50 

The choice of code which is implemented by the encoder 
30 has significant implications for the efficient operation of 
the transmitter 46. A multi -carrier signal comprises the sum 
of several sinusoidal waves with equal amplitudes and 
equally spaced frequencies. As a result, the maximum abso- 55 
lute value of the composite signal is strongly sensitive to the 
relative phases of the component sinusoids. These relative 
phases are changed at the end of each symbol period by the 
phase shift controllers 32, according to the value of the code 
word to be transmitted during the next symbol period. 60 
Depending on the particular set of relative phases, and 
therefore on the value of the current code word, the ampli- 
tude of the transmitted signal will have occasional pro- 
nounced peaks together with intervals of much lower though 
still varying value, as illustrated in FIG. 7. Thus the peak- 65 
to-mean envelope power ratio of the transmitted signal is 
relatively high. 
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For such a signal to be transmitted without distortion, it is 
necessary for the transmitter 46 to use a linear amplifier and 
to operate it in a relatively inefficient manner, so that most 
of its dynamic range is unused except for the infrequent 
occasions when the amplitude of the transmitted signal 
reaches its peak value. 

It has been proposed to reduce this problem by appropri- 
ate choice of the set of code words which are actually used 
to encode data for transmission. The use of block coding to 
avoid the use of code words which produce extreme values 
of PMEPR has been suggested in 'Minimisation of the peak 
to mean envelope power ratio of multicarrier transmission 
schemes by block coding' by T. A Wilkinson and A. E. 
Jones, 1995 IEEE 45 th Vehicular Technology Conference , 
July 1995, pp. 825-829; although this proposal indicates that 
reductions in PMEPR are possible, it also identifies the 
difficulty of selecting a set of code words which both enables 
efficient implementation of coding and decoding and pro- 
vides desired capabilities for detecting and correcting errors 
in the demodulated signal. Possible ways of addressing his 
problem are discussed, but no actual solution is presented. A 
modification to this approach is disclosed in 'Combined 
coding for error control and increased robustness to system 
nonlinearities in OFDM' by the same authors in 1996 IEEE 
46th Vehicular Technology Conference, April 1996, pp. 
904-908. In this modification a linear block code is chosen 
to provide desired properties of error correction and detec- 
tion; it is proposed then to exploit the linearity and redun- 
dancy in the code to transform the code values systemati- 
cally to new values which reduce the PMEPR but have 
equivalent error detection and correction properties. 
Examples are given for 4B/7B, 4B/8B and 11B/15B codes, 
in which the required transformation is identified by a 
procedure involving listing all possible code word values in 
order of the PMEPR for a multi-carrier signal modulated in 
accordance with each code word. Although such an 
exhaustive-search approach is feasible for relatively short 
code words (e.g. of the order of 15 bits), the number of 
possible code words to be checked increases very rapidly 
with code word length. Thus, even for a 32-bit code word, 
in excess of 4,000,000,000 code words would potentially 
need to be listed, and this would be at best extremely 
time-consuming and in practice possibly prohibitive. The 
obstacle presented by this requirement is potentially 
increased by the absence of any guarantee that an appropri- 
ate subset of code words actually exists within the set being 
examined. 

Popovic, in the paper cited earlier, proposes the use of a 
single binary or polyphase complementary sequence. 
However, the 1996 paper mentioned above suggests that 
such sequences are not readily amenable in relation to their 
error correction/detection capability. 

Thus the problem remains of identifying a set containing 
many different code words all of which satisfy desired 
criteria regarding error correction and detection properties, 
the maximum PMEPR of multi-carrier signals modulated in 
accordance with those code words, and ease and efficiency 
of implementation of a practical encoder and decoder (for 
example including use of combinatorial logic). The non- 
trivial scale of this problem may be seen by considering the 
example of 5-bit data words to be encoded using S-bit code 
words, as required in FIG. 5. In the general case there are 
256 possible code words, of which 32 are to be selected to 
encode the 32 possible data words. In the absence of 
knowledge of how particular combinations of code words 
might affect compliance with the criteria noted above, all 
possible subsets of 32 code words out of 256 would need to 
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be considered. There are 256!/(224!.32!) possible subsets, 
i.e. of the order of 10 41 . If 8-bit data words are to be encoded 
using 16-bit code words, there are 65536!/(65280!.256!) 
possible subsets to be considered, which is of the order of 
10 726 . 

FIG. 8 shows a 5B/8B encoder 60 and a decoder 80, 
which have the following properties: 

error detection: the minimum Hamming distance (i.e. the 
minimum number of symbols whose value must be 
altered, for example by interference, to transform a 
valid code word into some other valid code word) is 2 
(more generally, for binary code words described 
below containing 2 m symbols, this distance is 2 m-2 ); 
the maximum PMEPR of multi -carrier signals modulated 
in accordance with the 8 -bit code words is 3 dB in 
terms of power level, i.e. 6 dBV; 
the block code with which data are encoded contains a 
sufficiently large number of different code words to 
enable a usefully high data transfer rate to be attained; 
the encoder 60 and more particularly the decoder 80 can 
be readily implemented using combinatorial logic. • 
Referring to FIG. 8, the encoder 60 receives five data bits 
at its inputs. Four of these inputs (for the four least signifi- 
cant bits) are connected to respective binary multipliers 62 
to 68, and the fifth (for the most significant bit) is coupled 
to a selector 70. Each multiplier 62 to 68 is arranged to 
multiply the input bit it receives by a respective fixed 8-bit 
value as follows:- multiplier 62: 0101 0101; multiplier 64: 
0011 0011; multiplier 66: 0000 1111; and multiplier 68: 1111 
1111 . In practice each multiplier may be implemented using 
eight logical-AND gates to combine the input bit separately 
with each bit of the respective fixed 8 -bit value. The outputs 
of these four multipliers are connected to a bitwise 
exclusive- OR circuit 72, which effectively performs bitwise 
binary addition (i.e. addition modulo 2 of bits of correspond- 
ing significance) of the 8-bit values it receives. 

Depending on whether the most significant input bit is 0 
or 1 in value, the selector 70 selects one of two possible 
eight-bit values: 0001 0010 or 0000 0110 respectively. The 
choice of these particular values as the ones from which 
selection may be made will be explained below. The selected 
value is supplied to a fifth input of the exclusive- OR circuit 
72, the output of which constitutes the 8-bit code word. The 
encoder 60 as shown in FIG. 8 encodes the thirty-two 
possible 5 -bit data words as the following 8-bit code words: 



TABLE 1 



5-bit data 


8-bit code 


5-bit data 


8-bit code 


00000 


0001 0010 


10000 


0000 0110 


00001 


1110 1101 


10001 


1111 1001 


00010 


0001 1101 


10010 


0000 1001 


00011 


1110 0010 


10011 


1111 0110 


00100 


0010 0001 


10100 


0011 0101 


00101 


1101 1110 


10101 


1100 1010 


00110 


0010 1110 


10110 


0011 1010 


00111 


1101 0001 


10111 


1100 0101 


01000 


0100 0111 


11000 


0101 0011 


01001 


1011 1000 


11001 


1010 1100 


01010 


0100 1000 


11010 


0101 1100 


01011 


1011 0111 


11011 


1010 0011 


01100 


0111 0100 


11100 


0110 0000 


01101 


1000 1011 


11101 


1001 1111 


OHIO 


0111 1011 


11110 


0110 1111 


01111 


1000 0100 


11111 


1001 0000 



The decoder 80 receives the 8 -bit code words at its inputs, 
which couple the bit values to an 8-bit 2nd-order Reed- 
Muller (RM) decoder 82. Such decoders are in themselves 
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known, such as the Reed decoding algorithm using majority 
logic decoding as described, for example, in 'The Theory of 
Error- Correcting Codes* by F. J. MacWilliams and N. J. A- 
Sloane, North-Holland, Amsterdam, 1986, pp.385-388. 

5 They implement simultaneous decoding, and detection and 
correction of errors affecting some number of bits of the 
input code word up to at least a limit related to the length of 
the code word. The RM decoder 82 produces two sets of 
output signals which together comprise the output value of 

10 the RM decoder. A 'second-order' set of signals 84, in this 
case three in number, identifies a 'family' or coset of output 
values, and a 'first-order* set of four signals 86 identifies a 
particular output value within that coset. These four first- 
order signals 86 directly provide the decoded values for the 

15 four least significant bits of the data word output by the 
overall decoder 80. The three second-order signals 84 are 
connected to a selector 88 which generates the value of the 
fifth, most significant, bit of the decoder's output data word. 
To this end the selector 88 compares the values of the three 

20 second-order signals 84, taken together in sequence as 
described below, with two stored 3-bit values. If the second- 
order signals 84 have the value 101 the selector provides a 
binary 0 as the value of the most significant output bit; if 
they have a value of 110, it provides a binary 1 . If the signals 

25 84 have any other value than 101 or 110, then the received 
8 -bit code word has been treated by the RM decoder 82 as 
being other than one appearing in Table 1 above, and this is 
indicated by a signal at a data error output 90 of the selector 
88. 

30 An alternative decoding technique, described below with 
reference to FIG. 11 in the context of 16-bit code words, 
only outputs valid code words (i.e. code words which can be 
output by the encoder). 

The error detection and correction capabilities of the 

35 2nd-order RM decoder 82 itself, which are well established, 
are inherited by the overall decoder 80 by virtue of the 
incorporation of the RM decoder 82 within it. Thus, as noted 
above, the code's minimum Hamming distance is 2 (=2 m ~ 2 
with m-3). The decoder 80 can detect up to 2 m ~ 2 -l errors 

40 affecting the bits constituting a code word; for m=3 the 
number of detectable errors is one. In the limiting case of an 
8-bit code word, chosen initially here for simplicity of 
explanation, the decoder 80 cannot actually be guaranteed to 
correct such an error, as the number of correctable errors in 

45 the general case is 2 m_3 -l, which for m«3 equals zero. For 
a single-bit error in the received code word, either of two 
valid 8-bit code words could in some cases be affected by a 
change in a single bit position to produce the same received 
code word. In practice the above decoding technique would 

50 arbitrarily select one of the valid possibilities, so the error 
may be (but cannot be guaranteed to be) corrected. For 
longer code words, discussed below, error correction is 
possible. 

It can be shown that a multi-carrier signal modulated in 
55 accordance with the 8 -bit code words in Table 1 has a 
PMEPR of no more than 3 dB. The decoder 80 operates 
using analytical techniques (i.e. taking advantage of system- 
atic patterns within the set of code words), and can be 
implemented, for example, using combinatorial logic at least 
60 in part; although the selector 88 as described is effectively a 
look-up table, this table does not hold the entire set of code 
words and this manner of implementation is not essential. In 
the code shown in Table 1 there are thirty-two code words 
available, which provides an acceptable data transfer rate 
65 (i.e, five data bits per code word). 

It should be noted that although the decoder 80 contains 
a 2nd-order RM decoder 82 which is in itself known, the 
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code in Table 1 does not comprise the 2nd-order Reed- 
Muller code (of length eight) RM(2,3) in itself. It is instead 
a particular, non-linear sub-code of the 2nd-order Reed- 
Muller code of length eight, and has the novel and highly 
advantageous combination of properties set out above. Such 5 
a code has not been previously recognized. Likewise the 
decoder 80 does not constitute a conventional use of the 
8-bit 2nd-order RM decoder 82, i.e. use to provide decoding 
of data encoded according to a 7B/8B 2nd-order Reed- 
Muller code. 10 

The 5B/8B code implemented by the encoder 60 and 
decoder 80 of FIG. 8 and having this desired combination of 
properties is derived in the following manner, which, as will 
be described, may readily be extended to produce codes with 
longer-length code words. The starting point for the 5B/8B 15 
code is the following generator matrix of 8-bit values which 
defines a 'base' code: 



(0000 1111) 
(0011 0011) 
(0101 0101) 

(1111 1111) 



on 



(0000 0011) 
(0000 0101) 
(0001 0001) 



(*1< 

(X 2 * 



Xi) 



20 



25 



Each row of this matrix is identified with labels X l9 X 3 , 
X' respectively for convenience of reference in the following 
description. 

The base code comprises all sixteen linear combinations 
of rows of this generator matrix, that is combinations of the 3Q 
form 



35 



where (B indicates bitwise modulo -two addition and t l9 ^ t 
a 3 and a' are coefficients which each take values of 0 or 1 . 
This base code is linear, i.e. the result of a bitwise exctusive- 
OR operation on any two of its code words is another code 
word in the code. A sub-code of this code is any subset of 
the code words in the base code (and in particular could be 
the entire code), and a linear sub-code is such a sub-code 40 
which is itself linear. In practice, a sub-code must contain a 
sufficient proportion of the code words in the base code to 
constitute a usable system. The actual proportion depends on 
a variety of parameters, in particular the length of the code 
word and the number of possible values for each symbol in 45 
a code word (two values in the case of binary). In the current 
case of binary code words of length eight, a linear sub-code 
comprising eight of the sixteen code words in the base code 
might suffice. For longer code words and/or more possible 
symbol values a smaller proportion than half (e.g. a quarter, 
eighth, sixteenth, etc.) could well suffice. 

An eight-bit generator array is created by combining the 
generator matrix with additional rows (in this case three) 
derived by bitwise multiplication of all possible pairs of the 
rows X 1( X2 and X 3 , that is X 1 *X 2 , X!*X 3 and X^*X 3 : 



50 



55 



60 



Bitwise multiplication, indicated by the symbol *, involves 
multiplying the bits at corresponding positions in the two 
values being multiplied, or equivalently performing a bit- 
wise logical AND operation on those two values. 

These additional rows are combined by a bitwise 
exclusive-OR operation (bitwise addition modulo 2) to 
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produce (in the case of code words of length eight) three 
second-order values, or second-order 'coset representa- 
tives'. A first coset representative is produced by combining 
rows (X^XJ and (X^Xg): 



(0001 0010) 



(Xi*X2) ®{X 2 *X 3 ) 



The remaining coset representatives are produced by iden- 
tifying all possible permutations of the subscripts of the row 
identifiers X 1? X^ and X 3 in the expression (X 1 *X 2 )© 
(X2*X3) which yield distinct results: 



(0001 0100) 
(0000 0110) 



{X l *X>)®(X 2 *X 3 ) 
iX l *X 2 )®(X l *X 3 ) 



Since (X 1 *X a ) has the same value as (X 2 *X 1 ), and (X 2 *X 3 ) 
has the same value as (X^XJ, it follows that (X 1 *X 2 )© 
(X2*X 3 ) and p( 3 *X^®QL l *X^ have the same value, and 
therefore are not distinct. On the same basis the permutation 
(X 1 *X 3 )©(X 3 *X 2 ) has been replaced above by the equiva- 
lent expression (X J " t X 3 )®(X 2 *X 3 ) for the sake of clarity in 
indicating which additional rows are involved in each coset 
representative. 

Finally these coset representatives are each combined by 
bitwise exclusive-OR operations with all code words in the 
base code to produce a total of 48 code words, constituting 
a set of three cosets of the base code: 



TABLE 2 


(X 1 *X 2 )©(X 2 *X3) 


(x 1 -x 3 )e(x 2 *x 3 ) 




0001 0010 


0001 0100 


0000 0110 


1110 1101 


1110 1011 


1111 1001 


0001 1101 


0001 1011 


0000 1001 


1110 0010 


1110 0100 


1111 0110 


0010 0001 


0010 0111 


0011 0101 


1101 1110 


1101 1000 


1100 1010 


0010 1110 


0010 1000 


0011 1010 


1101 0001 


1101 0111 


1100 0101 


0100 0111 


0100 0001 


0101 0011 


1011 1000 


1011 3110 


1010 1100 


0100 1000 


0100 1110 


0101 1100 


1011 0111 


1011 0001 


1010 0011 


0111 0100 


0111 0010 


0110 0000 


1000 1011 


1000 1101 


1001 1111 


0111 1011 


0111 1101 


0110 1111 


1000 0100 


1000 0010 


1001 0000 



All of these 48 code words provide the combination of 
properties described earlier and in the most general case any 
selection from among these 48 code words can be made. In 
practice it is preferable to use a number of code words which 
is an integer power of two, thereby simplifying division of 
the incoming data stream into blocks for encoding, and 
furthermore to select code words in groups of sixteen having 
the same coset representative, to facilitate implementation. 
To this end the number of coset representatives chosen is the 
highest power of two that is no greater than the number of 
coset representatives available; in this case there are three 
available, and the corresponding highest power of two is 2 1 , 
that is 2. Choosing two coset representatives provides 32 
code words; the maximum corresponding binary data word 
length for which all data words can be encoded is therefore 
five bits, as shown in Table 1 . 

In the case of the encoder 60 and decoder 80 shown in 
FIG. 8., the code words having the coset representatives 
(X 1 *X 2 )©(X 2 *X 3 ) (i.e. 0001 0010) and (X 1 *X 2 )©(X 1 *X 3 ) 
(i.e. 0000 0110) have been arbitrarily selected, as indicated 
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(0000 0000 


mi 
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(*l) 


(0000 1111 


0000 


mi) 


(X 2 ) 


(0011 0011 


0011 


0011) 




(0101 0101 


0101 


0101) 


iX 4 ) 


(1111 1111 


mi 


mi) 


(X'J 
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within the selector 70. It can also be seen that the values 
selected by the selector 70 in accordance with the value of 
the most significant bit of the input data word are themselves 
these two coset representatives. The assignment of a par- 
ticular coset representative to a particular input bit value, or 5 
indeed to a particular input bit position (e.g. the most 
significant bit), is arbitrary, but must be consistent as 
between the selector 70 in the encoder 60 and the selector 88 

in the decoder 80. The three second -order signals 84 pro- The 16-bit base code comprises all thirty-two linear com- 

vided by the RM decoder 82 correspond to each of the 10 binations of rows of this generator matrix, that is combina- 

additional rows (X/XJ, (X/X^ and (X 2 *X 3 ) in the 8-bit tions of the form 
generator array described above. The individual bits of the 

3 -bit values stored in the selector 88 correspond to each of v ~ v m v m „ _ 
these rows respectively. Each 3-bit value is associated with 

a specific value for the most significant data bit, and indi- 15 , , . . , , -> c * 

cates which pair of the three additional rows is present in the whcre a - **> **> a < ™ d a each take values of 0 or 1 ' 

coset representative assigned (in the encoder 60) to that data Th e corresponding generator array is created by combin- 

bit value. Thus the data bit value 0 has assigned to it the ing the 16-bit generator matrix with additional rows (in this 

coset representative (X 1 ,*X 2 )(B(X 2 *X 3 ); correspondingly case six) derived by bitwise multiplication of all possible 

the presence in this coset representative of the additional 20 pairs of the rows X 2 , X^, X 3 and X^: 
rows (X^XJ and (X^Xg) is indicated by the positions of 

the 1 bits in the 3-bit value 101 stored in the selector 88 for (0000 0000 0000 nil) *x 2 ) 

that data bit value. ( 0ooo 0000 oon 0011) (*i*x 3 ) 

It can also be seen in FIG. 8 that the four least significant 25 (oooo oooo oioi oioi) *x A ) 

data bits are multiplied in the encoder 60 by respective ones (oooo 001 1 oooo 0011) (X 2 *x 3 ) 

of the rows X 3 , X^ Xj and X* of the generator matrix. The (0000 0101 0000 0l01) (X 2 *x 4 ) 

first-order signals 86 provided by the RM decoder 82 (Q001 om om Qm) 
correspond to the coefficients a^ a 2 , a l7 a' for each of these 

rows. The assignment of particular input data bits to each of 30 

these rows is arbitrary, provided the same assignment is used additional rows are combined by bitwise exclusive- 

in both the encoder 60 and the decoder 80. In effect this 0R ^rations, to produce (in the case of code words of 

means that there are many possible assignments of particular length sixteen) twelve second-order coset representatives. A 

8-bit code words to particular 5-bit data words, governed by fifst cosct escntative ^ duccd by combining rows 

the assignment of individual input data bits to each of these 35 *X^) fX,*X ) and (X *X J .\ 

rows, by the choice of which data bit controls selection of a 1 ^' ^ 3 
coset representative in the selector 70, and by the choice of 

coset representatives from which that selection may be < 0001 0010 0001 1101 > (***i>eCXa.*)©<*.*4) 
made. 

As noted above, the code words in Table 2 constitute a set 40 « 4 , , . , , , ^ c . 

of cosets of a linear sub-code (in this case the whole code) ™ e othe u r , roset "P'««^ v « are produced by identifying 

having a generator matrix comprising the rows X 1( X,, X 3 * possible permutations of the subscripts of the row iden- 

and X'. Using a base code of this kind provides desired error ^ fiers X » 3" *? . m ■ ex P ression ( X i X*)© 

detection and correction properties. Selection of a linear (XSXJ&X**Xd which yield distinct results: 

sub-code ensures that encoding and decoding within each 45 

coset can be implemented compacdy, for example using (° 001 0100 0001 1011 > (X l *x 2 )®(X 2 *x A )®(Xs*x 4 ) 

combinatorial logic. A coset can in general be considered to (° 000 °U U 0011 0l01 ) (#1**3)© (*2* #3) $(*2**4) 

be a systematic shift or translation of the base code to a (0001 0100 0010 0111) (X l *x 3 )®{Xj*x 4 )e{X 2 *x A ) 

different part of the overall space of possible code words. (0000 0110 0101 0011) (^ *X 4 )®{X 2 *X A )®{X Z *x 3 ) 

Use of cosets of the linear sub-code enables avoidance of 50 (0001 0010 0100 0111) (^ *X 4 )®{X 3 *X A )®{X 2 *X 3 ) 

certain undesired code words (such as all zeroes and all (0QO1 0001 0010 1101) { x l *x 2 )(B[x l *x y )^iXy*x 4 ) 

™^l~ h .ich W0Ul ? ^ult in particularly high values of (0001 om Qm l0U) ( ^ Xl)9W , we(ft , w 



(oooo 0011 0110 oioi) (x 2 *x 3 )®{x l *x i )e(x i *x 4 ) 



PMEPR, but can also enable the error control properties of 
the base code to be retained to some extent. Choosing the 

code words from a set of such cosets allows good control of 55 < 0000 0101 0110 001,) (*•*«>•<*•*)•(*•*> 

PMEPR because, as has not previously been recognized, (° 000 0101 0011 l001 > <* 

each coset of the base code tends to contain code words (000 ° 0011 0101 1001 >. (x 2 *x 3 )B(x l *x 1 )B(x i *x<) 
having broadly similar associated PMEPRs. In the particular 

case where the coset representatives have the form identified _ . 

above, it can be shown that the PMEPR does not exceed 3 60 before ' some Potations have been replaced above by 

equivalent expressions for the sake of clarity in indicating 

tt, j- j j j- u-r** c*w - 4* which additional rows are involved in each coset represen- 

I ne encoding and decoding capabilities or this invention . r 

can be applied to longer code words as follows. To retain the talive - 

benefits of error control and ease of implementation, it is Finally these coset representatives are each combined by 

preferable for the code word length to be an integer power 65 bitwise exclusive-OR operations with all code words in the 

of two. Thus the next appropriate code word length above 8 16-bit base code to produce a total of 384 code words, 

is 16. For 16 -bit code words the generator matrix is: constituting a set of twelve cosets of the base code: 
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TABLE 3 



(x^x^efx^Xa)©^-^) (Xj-x^acx^x^epvxo (x/x^ecx^x^ecx^x^) 



0001 0010 0001 1101 
1110 1101 1110 0010 
0001 0010 1110 0010 
1110 1101 0001 1101 
0001 1101 0001 0010 
1110 0010 1110 1101 
0001 1101 1110 1101 
1110 0010 0001 0010 
0010 0001 0010 1110 
1101 1110 1101 0001 
0010 0001 1101 0001 
1101 1110 0010 1110 
0010 1110 0010 0001 
1101 0001 1101 1110 
0010 1110 1101 1110 
1101 0001 0010 0001 
0100 0111 0100 1000 
1011 1000 1011 0111 
0100 0111 1011 0111 
1011 1000 0100 1000 
0100 1000 0100 0111 
1011 0111 1011 1000 
0100 1000 1011 1000 
1011 0111 0100 0111 
0111 0100 0111 1011 
1000 1011 1000 0100 
0111 0100 1000 0100 
1000 1011 0111 1011 
0111 1011 0111 0100 
1000 0100 1000 1011 
0111 1011 1000 1011 
1000 0100 0111 0100 



0001 0100 0001 1011 
1110 1011 1110 0100 
0001 0100 1130 0100 

mo ion oooi ion 

0001 1011 0001 0100 
1110 0100 1110 1011 
0001 1011 1110 1011 
1110 0100 0001 0100 
0010 0111 0010 1000 
1101 1000 1101 0111 
0010 0111 1101 0111 
1101 1000 0010 1000 
0010 1000 0010 0111 
1101 0111 1101 1000 
0010 1000 1101 1000 
1101 0111 0010 0111 
0100 0001 0100 1110 
1011 1110 1011 0001 
0100 0001 1011 0001 
1011 1110 0100 1110 
0100 1110 0100 0001 
1011 0001 1011 1110 
0100 1110 1011 1110 
1011 0001 0100 0001 
0111 0030 0111 1101 
1000 1101 1000 0010 
0111 0010 1000 0010 
1000 1101 0111 3101 
0111 1101 0131 0030 
3000 0030 1000 3103 
0111 1101 1000 1101 
1000 0030 0131 0030 



OOOO 0110 0013 0301 
3311 1001 1100 1 010 
OOOO 0110 1100 1010 
3311 1001 0011 0101 
0000 1001 0011 1010 
1111 0310 1100 0101 
0000 1001 1300 0103 
1111 0110 0013 1030 
0011 0101 0000 0110 
1100 1010 1111 1001 
0011 0301 3311 1001 
3100 3010 0000 0 330 
0011 1010 0000 1001 
1100 0101 1111 0110 
0011 1010 1111 0110 
1100 0101 0000 1001 
0101 0011 0110 0000 
1010 1100 1001 1111 
0103 0031 1001 1111 
1010 1100 0330 0000 
0301 1100 0110 1113 
1010 0011 1001 0000 
0101 1100 1001 0000 
1010 0011 0110 1113 
0310 0000 0101 0011 
1001 1111 1010 1100 
0110 0000 1010 1100 
1001 1111 0101 0011 
0110 1111 0101 1100 
1001 0000 1010 0013 
0110 1113 3030 0031 
1001 0000 0101 1100 



(Xj-Xs^CXa-X^eCX^Xj (X/X^OCX^X^OCX^Xs) CX^X^CXa'X^apVXa) 



0001 0100 0010 0111 
1110 1011 1101 1000 
0001 0100 1101 1000 
1110 1011 0010 0111 
0001 1011 0010 1000 
1130 0100 1101 0111 
0001 1011 1101 0111 
1110 0100 0030 1000 
0030 0111 0003 0100 
1101 1000 1110 1011 
0010 0111 1110 1013 
3101 1000 0003 0100 
0010 1000 0003 3033 
3101 0111 1110 0100 
0010 1000 1110 0100 
1101 0111 0001 1011 
0100 0001 0111 0010 
1011 1110 1000 1101 
0100 0001 1000 1101 
1011 1110 0113 0010 
0100 1330 0311 1101 
1033 0003 3000 0010 
0100 3310 1000 0010 
1011 0001 0111 1101 
0113 0030 0100 0001 
3000 3101 1011 1110 
0111 0010 1011 1110 
1000 1101 0100 0001 
0111 1101 0100 1110 
1000 0010 1011 0001 
0111 1101 1011 0001 
1000 0010 0100 1110 



0000 0110 0101 0011 
1111 1001 1010 1100 
0000 0110 1010 1100 
1311 1003 0101 0011 
0000 1001 0101 1100 
1311 0130 1010 0011 
0000 1001 1010 0011 
1111 0330 0101 1100 
0011 0101 0110 0000 
1100 1010 1001 1111 
0011 0303 1001 1111 
1100 1030 0330 0000 
0033 1030 0330 1111 
1300 0303 3001 0000 
0011 1010 1001 0000 
1100 0101 0110 1111 
0101 0011 0000 0110 
1010 1100 1111 1001 
0101 0011 1111 1001 
1010 1100 0000 0110 
0101 1100 0000 1001 
1010 0011 1111 0110 
0101 1100 1111 0110 
1010 0011 0000 1001 
0310 0000 0033 0103 
3003 1133 3300 3030 
0310 0000 1100 1010 
1001 1111 0011 0101 
0110 1111 0011 1010 
1001 0000 1100 0101 
0110 1111 1100 0101 
1001 0000 0011 1010 



0001 0010 0100 0111 
1110 1101 1011 1000 
0001 0010 1 011 1000 
1110 1101 0100 0111 
0001 110 1 0100 1000 
1310 0030 3033 0333 
0003 3303 3033 0333 
3310 0010 0300 3000 
0010 0001 0331 0100 
1101 1110 1000 1011 
0010 0001 1000 1011 
1101 1110 0111 0100 
0010 1110 0111 1011 
1303 0003 1000 0100 
0010 1130 3000 0300 
3303 0003 0333 3033 
0300 0333 0001 0010 
3031 1000 1330 3303 
0300 0331 1330 3303 
3013 3000 0003 0030 
0300 3000 0003 1101 
1011 0111 1110 0010 
0100 1000 1110 0030 
3033 0311 0001 1101 
0111 0100 0010 0001 
1000 1011 1101 3310 
0113 0300 1103 1110 
1000 1011 0010 0001 
0111 1011 0010 1110 
1000 0100 1101 0003 
0111 1011 1101 0001 
1000 0100 0010 1130 



(x 4 •x^ptj-x^ep^'x*) (x^x^cx^x^ecx^x*) pc 2 *x 3 )©(x 1 *x 3 )©(x 1 *x 4 ) 



0003 0003 0030 3101 
3330 3330 3101 0010 
0003 0003 1301 0030 
3130 3330 0030 1301 
0001 1110 0010 0010 
1110 0001 1303 1101 
0001 1110 1101 1101 



0001 0001 0100 1011 
1110 1110 1011 0100 
0001 0001 1011 0100 
1110 1110 0100 1011 
0001 1110 0100 0100 
1110 0001 1011 1011 
0001 1110 1011 1011 



0000 0011 0110 0101 
1111 1100 1001 1010 
0000 0011 1001 1010 
1113 1300 0330 0303 
0000 1100 0110 1010 
1111 0011 1001 0101 
0000 1100 1001 0101 
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TABLE 3-continued 



1110 0001 0010 0010 
0010 0010 0001 1110 
1101 1101 1110 0001 
0010 0010 1110 0001 
1101 1101 0001 1110 
0010 1101 0001 0001 
1101 0010 1110 1110 
0010 1101 1110 1110 
1101 0010 0001 0001 
0100 0100 0111 1000 
1011 1011 1000 0111 
0100 0100 1000 0111 
1011 1011 0111 1000 
0100 1011 0111 0111 
1011 0100 1000 1000 
0100 1011 1000 1000 
1011 0100 0111 0111 
0111 0111 0100 1011 
1000 1000 1011 0100 
0111 0111 1011 0100 
1000 1000 0100 1011 
0111 1000 0100 0100 
1000 0111 1011 1011 
0111 1000 1011 1011 
1000 0111 0100 0100 



1110 0001 0100 0100 
0010 0010 0111 1000 
1101 1101 1000 0111 
0010 0010 1000 0111 
1101 1101 0111 1000 
0010 1101 0111 0111 
1101 0010 1000 1000 
0010 1101 1000 1000 
1101 0010 0111 0111 
0100 0100 0001 1110 
1011 1011 1110 0003 
0100 0100 1110 0001 
1011 1011 0001 1110 
0100 1011 0001 0001 
1011 0100 1110 1110 
0100 1011 1110 1110 
1011 0100 0001 0001 
0111 0111 0010 1101 
1000 1000 1101 0010 
0111 0111 1101 0010 
1000 1000 0010 1101 
0111 1000 0010 0010 
1000 0111 1101 1101 
0111 1000 1101 1101 
1000 0111 0010 0010 



1111 0011 0110 1010 
0011 0000 0101 0110 
1100 1111 1010 1001 
0011 0000 1010 1001 
1100 1111 0101 0110 
0011 1311 0101 1001 
3100 0000 1010 0110 
0011 1111 1010 0110 
1100 0000 0101 1001 
0103 0330 0013 0000 
1010 1001 1300 1113 
0103 0110 3300 1131 
1010 1001 0011 0000 

oioi iooi ooii mi 

1010 0110 1300 0000 
0303 3003 3300 0000 
3010 0130 0011 1111 
0110 0101 0000 0011 

iooi ioio mi iioo 
oiio oioi mi noo 

1003 1010 0000 0011 
0310 1010 0000 1100 
1001 0301 3331 0011 
0130 3030 3111 0011 
1001 0101 0000 1100 



(x^x^cx^ecx, *x 3 ) (x, *x 3 )© pc^x^epc^x,) pc 2 *x 3 )©(x 1 *x 2 )©pc 1 *x<) 



0000 0103 0110 0011 
1333 3010 3003 3300 
0000 0303 3003 3300 
3311 1010 0130 0033 
0000 1010 0330 3300 

mi oioi iooi ooii 

0000 1010 1001 0011 
1131 0103 0310 1100 
0011 0130 0303 0000 
3100 1003 1010 1111 
0011 0110 3030 3111 

noo iooi oioi oooo 
oon iooi oioi nn 
noo oiio ioio oooo 

0011 1001 1010 oooo 
1100 0110 0101 1113 
0303 0000 0031 0110 
1030 1113 3300 3001 
0101 0000 3300 1001 
1010 1333 0013 0110 
0103 3333 0013 3001 
1030 OOOO 3300 0310 
0303 3333 3300 0330 
1030 OOOO 0033 3003 
0330 0011 OOOO 0101 
1001 3300 3333 3010 
0130 0011 1113 3010 
1001 1100 0000 0101 
0110 1100 OOOO 1030 
3001 0031 1111 0101 
0110 1100 1113 0301 
1001 0013 OOOO 3010 



OOOO 0101 0011 1001 

nn ioio noo ono 

OOOO 0101 1100 0110 
1111 1010 0011 1003 
OOOO 3010 0033 0330 
1111 0301 3300 3003 
OOOO 3010 1100 1001 

nn oioi 0013 ono 
oon ono oooo ioio 
noo iooi mi oioi 
oon ono nn oioi 
noo iooi oooo ioio 

0013 3003 OOOO 0303 
3100 0310 1133 1010 

oon iooi mi ioio 

1100 0110 0000 0101 
0101 OOOO 0310 1100 

ioio mi iooi oon 

0101 0000 1001 0011 
1010 1131 0110 1100 
0103 3131 0130 0011 
1030 OOOO 1001 3300 
0303 3333 3003 3300 
3030 OOOO 0330 0033 
0330 0013 0303 3333 
3003 3300 3030 OOOO 
0330 0033 3030 OOOO 
3003 3300 0303 3333 
0110 1100 0101 OOOO 
1001 0031 3010 1111 
0130 3300 3030 3111 
3003 0013 0303 OOOO 



OOOO 0033 0301 1001 

nn noo ioio ono 
oooo oon ioio ono 

1111 1300 0303 3003 
OOOO 3300 0103 0330 
3333 0033 3030 1003 
OOOO 3300 3030 3003 
3333 0033 0303 0330 
0033 OOOO 0130 3030 
1300 3333 3003 0101 
0013 OOOO 3003 0301 
1300 3333 0330 3030 
0033 3333 0330 0303 
3300 OOOO 3003 3030 
0033 1333 3003 3030 
3300 OOOO 0330 0303 
0303 0110 OOOO 1100 
1010 1001 1113 0033 
0301 0310 3313 0033 
3030 3001 OOOO 3300 
0301 3001 OOOO 0033 
1030 0330 3111 1100 

oioi iooi nn noo 

3030 0330 OOOO 0033 
0330 0303 0033 3333 
3003 3030 3300 OOOO 
0310 0301 1100 OOOO 
1003 1030 0033 3313 
0330 3030 0013 OOOO 
1001 0101 1300 3333 
0330 3030 3300 3133 
3003 0103 0013 OOOO 



All of these code words provide a combination of prop- 
erties corresponding to that described earlier, i.e.: 55 

the code's minimum Hamming distance is 4 (=2 m " a with 
m=4); an error-correcting decoder can detect up to three 
(2 m ~ 2 -l) errors affecting the bits constituting a code 
word, and also correct one (2 m_3 -l=l with m=4) such 6Q 
error; 

a multi-carrier signal modulated in accordance with the 
16-bit code words in Table 3 has a PMEPR of no more 
than 3 dB; 

a decoder for use with this code can operate using 65 
analytical techniques and can be implemented, for 
example, using combinatorial logic at least in part; 



there are sufficient code words available to enable an 
acceptable data transfer rate to be attained. 

For the preferred case (using a number of code words 
which is an integer power of two, and selecting code words 
in groups of thirty-two having the same coset representative) 
the highest power of two no greater than twelve (the number 
of coset representatives) is 2 3 , that is eight. Thus in the 
preferred case eight of the cosets of 32 code words each 
would be chosen, providing 256 code words, corresponding 
to a binary data word length of 8 bits. 

FIG. 9 shows an 8B/16B encoder 100 for encoding 8-bit 
data words as 16-bit code words in accordance with this 
invention, using code words from Table 3. This encoder 
includes five binary multipliers 102 to 110 for multiplying 
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the five least significant bits of the input data word by 
respective ones of the 16-bit generator matrix rows X 4 , X 3 , 
X2, Xj and X'. The three most significant bits of the input 
data word are supplied to a selector 112 which selects one of 
eight possible coset representative values in accordance with 5 
the combined value of those three bits. In the example 
shown in FIG. 9 the first eight of the twelve possible coset 
representatives listed in Table 3 have been chosen for use in 
the selector 112. The outputs of the five binary multipliers 
102 to 110 and of the selector 112 are combined by a bitwise 10 
exclusive- OR circuit 114 to produce the 16-bit code word. 

FIG. 10 shows a first example corresponding decoder 120 
for converting the 16-bit code words back to the correspond- 
ing data words. The code words are supplied to a 16-bit 
2nd-order RM decoder 122, corresponding to the 8-bit RM 15 
decoder 82 in FIG. 8 but for a larger value of m. Five 
first-order signals 124 produced by the RM decoder 122 
directly provide the decoded values for the five least sig- 
nificant bits of the data word output by the overall decoder 
120; these signals 124 correspond to the coefficients a 4) 20 
a 2 , a a and a' for the rows X 4 , X 3 ,X 2 ,X 1 and X' of the 16-bit 
generator matrix. Six second-order signals 126 from the RM 
decoder 122, corresponding to the six additional rows 
(Xj'Xa), (X/Xa), (X/X4), (X^), (X^X 4 ) and (X 3 *X 4 ) 
of the generator array, are connected to a selector 128 which 25 
generates the values of the three most significant bits of the 
overall output data word. To this end the selector 128 
compares the values of the six second-order signals 126, 
taken together, with the eight 6 -bit values shown in FIG. 10; 
these 6 -bit values could be stored, or derived by an appro- 30 
priate procedure implemented in software or firmware as 
required. Each 6-bit value is associated with a specific 
combination of values for the three most significant data 
bits, and indicates which triplet of the six additional rows is 
present in the coset representative assigned (in the encoder 35 
100) to that combination of three data bit values. If the 
second-order signals 126 have one of the eight stored values, 
the selector 128 provides the corresponding three data bit 
values for the most significant output bits. If the signals 126 
have any other value, then the received 16-bit code word has 40 
been treated by the RM decoder 122 as being other than one 
in the first eight cosets in Table 3, and this is indicated by a 
signal at a data error output 130 of the selector 128. 

Any decoder for an error-correcting code of which the 
minimum Hamming distance is d, such as the decoder 120 45 
(for which d»2 m " 2 ), correctly decodes a received code word 
whose distance from the transmitted code word is less than 
d/2. However, the decoder 120 is not a * minimum distance* 
decoder; a minimum distance decoder has the additional 
property that it decodes to whichever of the possible trans- 50 
mitted code words has the least distance from the received 
code word, even if this distance is d/2 or more. Furthermore, 
occurrence of a signal at the data error output 130 indicates 
that the second-order Reed-Muller code word output by the 
RM decoder 122 is not one that can actually be output by the 55 
encoder 100, which is an indeterminate result 

FIG. 11 shows an alternative 16-bit decoder 140 which 
only outputs valid code words (i.e. code words which can be 
output by the encoder 100) and furthermore is a minimum 
distance decoder with respect to this set of code words. 60 

The method used in the decoder 140 is an example of the 
"Supercode" decoding method described in 'Soft decoding 
techniques for codes and lattices, including the Golay code 
and the Leech lattice' by J. H. Conway and N. J. A. Sloane, 
IEEE Trans. Inform. Theory, 1986, vol.IT-32, no. 1, pp. 65 
41-50. For each coset representative in turn the received 
code word is exclusive-ORed with that coset representative, 
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and the result is decoded as an element of the first-order 
Reed-Muller code of length sixteen, RM(1,4), at the same 
time assigning a "score" to that decoding. The decision as to 
which code word to decode to is made by selecting a coset 
whose representative attains the highest score (arbitrarily if 
more than one attains that score), and within that coset 
selecting the element of RM(1,4) previously identified for 
that coset. 

The decoding of an element of RM(1,4) uses the fast 
Hadamard transform, as described in 'The Theory of Error- 
Correcting Codes' by F. J. MacWilliams and N. J. A. Sloane, 
North-Holland, Amsterdam, 1986. The fast Hadamard trans- 
form on a binary sequence of length 2 m produces 2 m 
transform components, each associated with a pair of ele- 
ments of the first-order Reed-Muller code of length 2 m , 
RM(l,m). The largest magnitude of all the components is the 
"score" for the decoding, and the decoding algorithm 
chooses the element of RM(1 ,m) corresponding to a com- 
ponent of largest magnitude (with the sign of that component 
determining which of the associated pair of code words is 
used). 

Referring to FIG. 11 the decoder 140 receives each 16-bit 
code word in an exclusive-OR gate 142 which performs 
bitwise subtraction modulo two of a value in a coset store 
144 from the code word. (In the case of modulo-two 
arithmetic, bitwise subtraction and addition are equivalent, 
and a bitwise exclusive-OR gate may be used for either 
operation; the function of the gate 142 is denned as sub- 
traction for consistency with subsequent description.) The 
coset store 144 contains the values of the same eight coset 
representatives from which selection is made by the selector 
112 in the encoder 100. The store 144 also has a three -bit 
counter 146 which cycles through all eight possible three-bit 
values successively in response to an increment signal at an 
input 148. The current three -bit value in this counter controls 
which of the eight coset representative values is supplied to 
the exclusive-OR gate 142, the correspondence between 
three-bit values and coset representatives being the same as 
in the encoder 100. Thus as the counter 146 cycles through 
its eight possible values, each of the eight coset represen- 
tatives in turn is provided to the exclusive-OR gate 142. 

The 16-digit binary word output by the gate 142 for the 
current coset representative is provided to a substitution 
circuit 150 which substitutes the values +1 and -1 for the 
binary values 0 and 1 respectively in that word, to match the 
element values in terms of which the fast Hadamard trans- 
form is customarily applied. 

The 16-symbol word from the circuit 150 is supplied to a 
so-called 'Green machine' 152 for performing the fast 
Hadamard transform (see the book by MacWilliams and 
Sloane referenced above, pp. 423-5), or a modification 
thereof (see, for example, 'Fast decoding algorithms for first 
order Reed-Muller and related codes'by A. E. Ashikhmin 
and S. N. Litsyn, Designs, Codes and Cryptography, vol.7, 
1996, pp. 187-214). The output of the Green machine 152 
comprises sixteen transform components labelled Eq to E, s 
(in the general case of a length 2 m code word the subscript 
15 is replaced by 2 m -l); the component E 0 corresponds to 
the left-most symbol of a row of the code's generator matrix, 
and the component E as corresponds to the right-most sym- 
bol. These components are supplied to a comparator 154 
which compares the absolute values of the transform 
components, and identifies the largest absolute value |Ej and 
the corresponding value of z (from 0 to 15); if more than one 
component has the same largest absolute value, one of them 
is chosen arbitrarily. The comparator then outputs the values 
E z and z to a latch 156, which also receives the correspond- 
ing three-bit index value supplied by the counter 146 in the 
store 144. 
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For the first coset representative the latch 156 simply 

stores the values for z and the coset index which it (0000 0000 0000 0000 nil nil till nil) (Xj) 

receives. As the counter 146 is incremented (for the same (oooo 0000 nil nil 0000 0000 nil nil) {X 2 ) 

16-bit input code word), each successive coset representa- (oooo mi oooo nn 0000 nil 0000 nil) (* 3 ) 

tive is supplied to the exclusive-OR gate 142, and respective s (0 on ooil 0011 0011 0011 0011 0011 ooil) [X 4 ) 

values of E z and z are provided to the latch 156. If such a (010J 0]01 0101 0101 0101 0101 0l01 0101) {X$) 

transform component E z has a larger absolute value than the (1UJ 11U Jin im im uu an 
component currently latched, the latch selects the newly 
received component, together with its associated values of z 

and the coset index. 10 The 32-bit base code comprises all sixty-four linear combi- 

When the counter 146 has cycled through all eight pos- nations of rows of this generator matrix, that is combinations 

sible values, the contents of the latch 156 are used to of me form 
assemble the decoded output data at 158. The value of z 

directly determines the values of four decoded row coeffi- 15 

cients a 4 , a 3 , a^ and a 3 when z is expressed as a four-bit where a a , a^, s^ t a 4 , a 5 

binary number a 1 a 2 a 3 a 4 ; a fifth decoded row coefficient a' is The corresponding generator array is created by combin- 

zero if the final latched value E z is positive and is one if this ing the 32-bit generator matrix with additional rows (in this 

value is negative; and the coset index comprises the three case ten) derived by bitwise multiplication of all possible 
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Although the implementation of the fast Hadamard trans- 
























form has been described above as taking place for each coset 
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case where decoding is carried out in hardware, the required 
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can be calculated in parallel, by providing a separate Green 
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The decoder 120 described previously has an advantage in 
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general, relative to the decoder 140, of faster operation. 


30 


(0000 


0011 


oooo 


0011 


oooo 


0011 


oooo 


0011) 


(*3- 




Accordingly it may be desired to use the decoder 120 as a 
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primary decoder, to benefit from its speed; if the signals 126 
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from the RM decoder 122 are found to correspond to a code 
word which cannot be output by the encoder 100, as 

indicated at the data error output 130, the received code and these additional rows are combined by bitwise 

word can then be supplied to the decoder 140 whose output exclusive-OR operations to produce (in the case of code 

is used instead. Such choice between the decoders 120 and words of lcn S th thirty-two) sixty second-order coset repre- 

140 can be adaptive, based on current noise levels as ™ ives - ™™ a fta c^repiese ntative is Produced by 

indicated by the frequency of occurrence of signals at the ^ombming rows (X/XJ, (X, X 3 ), (X 3 *X 4 ) and (X« *X£ 

, A * * I™ rr L • 1 i.-i c the remaining coset representatives are produced by ldenti- 

data error output 130. If such signals occur relatively often, 40 r_ . t1 4 4 . c.t. i_ • . c 

. mi . . . . A . . , rying all possible permutations 01 the subscripts or the row- 

indicating a currently noisy communications path, the identifiers x ^ ^ x ^ x in the express io n (X/X,) 

decoder 140 can be selected for immediate use for each code ©(X 2 *X 3 )©(X 3 *X 4 )S(X 4 *X 5 ) which yield distinct results, 

word, whde continuing to use the decoder 120 m parallel to ^ representatives are eac h combined by bitwise 

monitor the communications quality. When the signals at the exclusive-OR operations with all code words in the 32-bit 

data error output 130 decline, the decoder 120 can be 45 base code to produce a tota ] 0 f 3340 possible code words, 

selected again as the primary decoder. constituting a set of sixty cosets of the base code. These 

It has not been previously recognized that the particular particular code words provide the following combination of 

384 16-bit code words listed in Table 3, out of the 65536 properties: 

possible 16-bit code words, have simultaneously all of the tne code's minimum Hamming distance is 8 (=2 m " 2 with 

different properties set out above (i.e. large minimum Ham- 50 m =5) ; ao error-correcting decoder can detect up to 

ming distance, low value of PMEPR, convenience of prac- seven (2 m-2 -l) errors affecting the bits constituting a 

tical implementation, and useful number of possible code code word, and also correct up to three (2 m-3 -l=3 with 

words) and thus provide significant and highly advantageous m»5) such errors; 

utility. Indeed, it has been considered that finding a set of ss a multi-carrier signal modulated in accordance with these 

code words, especially for longer word lengths, which particular 32-bit code words has a PMEPR of no more 

combined no more than desired error control properties with than 3 dB; 

low PMEPR, constituted a problem which was at best a decoder for use with this code can operate using 

intractable and for which there was not even any certainty analytical techniques and can be implemented, for 

that a solution existed. The present invention provides not 60 example, using combinatorial logic at least in part; 

only this combination, but also convenience of implemen- there are sufficient code words available to enable an 

tation and a code word population of adequate size to enable acceptable data transfer rate to be attained, 

a satisfactory rate of data transfer to be achieved. It has not previously been recognized that these particular 

Code sets with the corresponding properties for longer 3840 32-bit code words (out of the more than 4,000,000,000 

code word lengths can be readily identified by extension of 65 possible code words of this length), have this wide range of 

the procedure described above to produce Table 3. Thus, for differing properties simultaneously, with their associated 

32-bit code words the generator matrix is: highly advantageous practical benefits. 
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For the preferred case (using a number of code words 
which is an integer power of two, and selecting groups of 
sixty-four code words each having the same cosete 
representative) the highest power of two no greater than 
sixty (the number of coset representatives) is 2 5 , that is 5 
thirty-two. Thus in the preferred case thirty-two of the cosets 
of sixty-four code words each would be chosen, providing 
2048 code words, corresponding to a binary data word 
length of 11 bits. In this case, six input data bits a 5 , a 4 , ^ 
a 2 , a a and a 1 are multiplied in an encoder by respective ones 10 
of the 32-bit generator matrix rows X 5 , X4, X 3 , X^, X x and 
X'; and the remaining five bits of the input data word are 
supplied to a selector for selection of one of the thirty-two 
chosen coset representative values, in accordance with the 
combined value of those five bits. For decoding, a 32-bit 15 
2nd-order RM decoder produces six first-order signals 
(comprising six bits of the decoded data word), and ten 
second -order signals coupled to a selector which stores 
thirty-two possible 10-bit values for those ten signals taken 
together and which generates the remaining five bits of the 20 
decoded 11 -bit data word. Alternatively, a decoder analo- 
gous to the decoder 140 of FIG. 11, using the fast Hadamard 
transform, may be used. 

In general terms the generator matrix for a particular code 
word length 2 m is: is 

((2.0 + l)x2 ff, )/2" ((2.1 + l)x2 m )/2 OT 
((2.0-»-l)x2 m - 1 )/2 m " 1 ((Zl-t-l)x2 m_1 )/2 ffl " 1 ••• 

30 

((2.0 + l)x2J')/2' ((2.1 + l)x2')/2J 
((2.0 + l)x2 l )/2 I ((ll + l)x2 1 )/2 l 

(C2.0 + l)x2°)/2° ((ll + l)x2°)/2° - ^ 

({2.z-H)x2 m )/2 OT — ((2.(2 m -l)+l)x2 m )/2 m 
((2z-H)x2"** l )/2 m ~ 1 •■• ((2.(2*- 1) -i-l) xir- 1 )/?- 1 

40 

((%z+\)xl y )!V •■■ ((2.(2 m -l)+Ux2 J ')/2* 

((2.z + l)x2 1 )/2 1 .-. ((2.(2** - 1) + L)x2 l J/2 1 

((2z + l)x2°)/2° ... ((2.(2 m -l)+l)x2°)/2° 45 

where z has integer values from 0 to 2™-l, y has integer 
values from m to 0, and x indicates bitwise multiplication 
(bitwise logical AND) of two factors (such as (2.Z+1) and 2 y ) 5Q 
both expressed as (m+l)-bit binary numbers, and the result 
of division is expressed as a single digit. The first line of this 
matrix (y=m) gives the row X l9 the second line (y=m-l) 
gives the row Xj, the penultimate line (y=l) gives the row 
X^ and the final line (y=0) gives the row X'. 55 

Thus, for example, for a code word length of sixty-four 
bits (m=6), the first symbol (z=0) on the second row X% 
(y-m-1-5) is 

((2.0+l)x2 M - 1 )/2'"- 1 

60 

that is 

(lx32)/32 

Expressing the factors involved in the bitwise multiplication 
as (m+l)-bit binary numbers results in 65 

(0000001x0100000)/32 
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which is (0000000)/32 or (expressed as a single digit) 0. The 
second symbol (z-1) on the second row is ((2.1+l)x2 m_1 )/ 
2 m " 1 which is (3x32)/32, which yields (0000011x0100000)/ 
32, or 0 again. The sixteenth symbol (z-1 5) on this row is 
((2.15+l)x2 m - 1 )/2 m * 1 which is (31x32)/32, which also 
yields (0011111 x0100000)/32, or 0. 

The seventeenth symbol (z=16) is ((2,16+l)x2 m " 1 )/2 m " 1 
which is (33x32)/32. This yields (0100001x01 00000)/32, 
which is (0100000)/32 or 1. The same result is obtained for 
the eighteenth through to the thirty-second symbols. 

For the thirty-third symbol (z-32) the matrix gives 
((2.32+l)x2 w - 1 )/2 m - 1 which is (65x32)/32 and this yields 
(1000001x0100000)/32, that is 0. A result of 0 is likewise 
obtained for the following fifteen symbols. 

The forty-ninth symbol (z-4£) is given by the expression 
((2.48+l)x2"- 1 )/2 ,M - 1 which is (97x32)/32 giving 
(1100001x0100000)/32, or 1, which is the result for the final 
fifteen symbols as well. Thus the second row in its entirety 
is 

(GOOD 0000 0000 0000 1111 1111 111 1111 0000 0000 0000 0000 

1111 1111 11111111) 

For the third row X 3 (y=m-2=4) the first symbol is 
((2.0+l)x2 m ~ 2 )72 m - 2 that is (lxl6)/16 which is (OOOOOOlx 
0010000)/16, or 0; the same result is produced for the next 
seven symbols. For the ninth symbol the expression is 
((2.8+l)x2 m - 2 )/2 m - 2 that is (17xl6)/16 and this yields 
(0010001x0010000)/16, or 1, and similarly for the tenth to 
sixteenth symbols. Thereafter the symbols continue in alter- 
nating blocks of eight zeroes and eight ones for the rest of 
the row: 

(0000 0000 1111 1111 0000 0000 1111 1111 0000 0000 1111 1111 
0000 0000 1111 1111) 

The other rows of the generator matrix for 64-bit code words 
would be derived in similar manner. 

The 2 m -bit base code comprises all linear combinations of 
rows of the generator matrix, that is combinations of the 
form 

where a 1P ag, . . . a m and a' each take values of 0 or 1. 

The corresponding generator array is created by combin- 
ing the 2 m -bit generator matrix with additional rows derived 
by bitwise multiplication of all possible pairs of the rows X 1( 
Xs, X 3 , . . . X m . These additional rows are combined by 
bitwise exclusive -OR operations to produce m!/2 second- 
order coset representatives. Thus a first coset representative 
is produced by combining rows (X^X^), (X^*X 3 ) y . . and 
(X m _ 1 *X m ); the remaining coset representatives are pro- 
duced by identifying all possible permutations of the sub- 
scripts of the row identifiers X 2 , X^, X3, . . . X^ in the 
expression (X 1 *X^®Qt 2 *X 3 )® . . . ©(X^/XJ which 
yield distinct results. 

These coset representatives are each combined by bitwise 
exclusive-OR operations with all code words in the 2 m -bit 
base code to produce a total of m!(2 m * 1 )/2 possible code 
words, constituting a set of m!/2 cosets of the base code. 
These code words have the following combination of prop- 
erties: 

the code's minimum Hamming distance is 2 m ~ 2 ; an error- 
correcting decoder can detect up to 2 m_2 -l errors 
affecting the bits constituting a code word, and also 
correct up to 2 m ~ 3 -l such errors; 

a multi-carrier signal modulated in accordance with these 
code words has a PMEPR of no more than 3 dB; 
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a decoder for use with this code can operate using 
analytical techniques and can be implemented, for 
example, using combinatorial logic at least in part; 

there are sufficient code words available to enable an 
acceptable data transfer rate to be attained. 

The total number of binary code words available for each 
of a range of code word lengths is shown below, together 
with the following other parameters: the number of coset 
representatives; the maximum number of data words which 
is an integer power of two and which can be encoded using 
the available number of code words; the data word length in 
that case; the number of lst-order signals and 2nd-order 
signals provided by an RM decoder for the relevant code 
word length; and the number of selector values stored in an 
encoder and decoder such as those shown in FIGS. 8, 9 and 
10: 
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(0000 0022) 2{X l * X 7 ) 

(0000 0202) 2{Xi * X y ) 

(0002 0002) 2{X 2 *Xy) 

The additional rows are combined by symbol -wise addi- 
tion modulo 2 to produce m!/2 second-order coset repre- 
sentatives. In the case of quaternary code words of length 
eight, a first coset representative is produced by combining 
1 rows ^X/Xa) and 2(Xz*X 3 ): 

(0002 0020) 2{Xi * X 2 ) © 2(X 2 * X 3 ) 

1 The remaining coset representatives are produced by iden- 
tifying all possible permutations of the subscripts of the row 
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15 
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For the sake of simplicity and clarity the above descrip- 30 
tion has been directed to coding using just two symbol 
values, i.e. binary encoding. However, the present invention 
is also applicable to coding with any higher number of 
symbol values which is an integer power j of two, such as 
quaternary (j-2; 2 / -4), octary (j-3; 2^-8) and higher orders. 
In such cases the same generator matrix is used as described 
above, e.g. for a code word containing 8 symbols (2 m =8 with 
m=3): 

40 

(0000 1111) (X,) 
(0011 0011) (X 2 ) 
(0101 0101) (Xi) 
(1111 1111) (X') 

45 

An 8-symbol quaternary base code (four possible symbol 
values) comprises all 256 linear combinations of the rows of 
this generator matrix, that is combinations of the form 

50 

where ® now indicates symbol-wise addition modulo 2? (i.e. 
modulo 4 for quaternary) and a 1( aj, a^ and a' are coefficients 
which each take all integer values from 0 to 2/-1 (i.e. 0 to 55 
3 for quaternary). This base code is linear, i.e. the result of 
symbol- wise addition modulo 4 on any two of its code words 
is another code word in the code. 

The corresponding generator array is created by combin- 6Q 
ing the generator matrix with additional rows. These addi- 
tional rows are derived by symbol-wise multiplication of all 
possible pairs of the rows Xj, X^ and X3, that is X 1 *X 2 , 
Xj*X 3 and X^X-j, and then further multiplication by 2'" 1 , 
e.g. by 2 for quaternary coding or by 4 for octary coding. $5 
Thus for 8-symbol quaternary coding the additional rows 
are: 



identifiers X 1? X2 and X3 in the expression 2(X 1 *X 2 )®2 
(X^Xg) which yield distinct results: 

(0002 0200) 2(Xi * Xy ) © 2(X 2 * X 3 ) 

(0000 0220) 2{X i * X 2 ) © %X V * X 3 ) 

Finally these coset representatives are each combined by 
symbol-wise addition modulo 2* with all code words in the 
base code to produce m!(2 Cm+! >)/2 code words, constituting 
a set of m!/2 cosets of the base code. Examples of the 768 
code values thus produced in the 8-symbol quaternary case, 
and of the corresponding values of the coefficients a a , a^ ^ 
and a', are given below, where CR 1 -2(X 1 *X 2 )02(X 2 *X3), 
CR 2 =2(X 1 *X 3 )®2(X 2 *X 3 ) and CR 3 =2(X 1 *X 2 )e2(X 1 *X 3 ): 



TABLE 4 



0002 0020 


cR 1 e(o.x 1 )©(o.x 2 )©(o.x 3 )e(o.x*) 


1113 1131 


cR 1 e(o.x 1 )©(o.x 2 )©(o.x 3 )e(i.x') 


2220 2202 


ra 1 e(o.x 1 )©(o.X2)©(o.x 3 )©(2,xo 


3331 3313 


CR 1 ©(O.X 1 )e(p.x3©(OJC^©(3.X0 


0103 0121 


CR 1 ffi(0.X 1 )©(0.X 2 )ffi(l.X 3 )©(0.X') 


1210 1232 


CR 1 ffi(QJC 1 )©(0.X^©(lJC^©(lJC^ 


2321 2303 


CR t ©(0 JC t )©(0 .x3©(l .X 3 )©(2,X^ 


2303 2321 


CR 1 ©(0.X 1 )©(2,X 2 )©(1.X 3 )©(2,X') 


1203 2332 


CR 1 ©(1.X 1 )©(3.X 2 )©(1.X 3 )©(1.X') 


1232 3032 


CR 1 ©(2,X 1 )©(2TQ©(1.X 3 )©(1.X') 


2112 1023 


CR 1 ©(3.X 1 )©(3.X 2 )©(3.X 3 )©(2,X') 


3223 2130 


CR 1 ©(3.X 1 )©(3.X 2 )©(3.X3)©(3.X0 


0002 0200 


CR 2 ©(O.X 1 )©(p.X 2 )©(O.X 3 )©(O.X0 


1113 1311 


CR 2 ©(0.X 1 )©(0.X 2 )©(0.X 3 )©(1.X') 


2112 1203 


CR 2 ©(3.X 1 )ffi(3.X 2 )©(3.X 3 )ffi(2.X') 


3223 2310 


CR 2 ©(3.X 1 )©(3.X 2 )©(3.X3)©(3.X0 


0000 0220 


CR 3 ©(0.X 1 )©(0.X 2 )ffi(0.X 3 )©(0.X') 
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TABLE 4-continued 



30 



1111 1331 



3221 2330 



CR 3 ©(0.X 1 )©(0.X 2 )©(0.X 3 )©(1 .X*) 
CR 3 ©(3.X 1 )©(3.X 2 )©(3.X 3 )©(3.X') 



The 2 m -symbol quaternary base code comprises all linear 
combinations of rows of the generator matrix, that is com- 
binations of the form 



10 



20 



25 



30 



35 



where a 5> a 2 , a^, . . . a^, and a' each take values of 0, 1, 2 or 
3. The corresponding generator array is created by combin- 
ing the 2 m -symbol generator matrix with additional rows 
derived by symbol-wise multiplication of aU possible pairs 
of the rows X l7 X3, . . . X m . These additional rows are 
combined by symbol-wise addition modulo 4 to produce 
m!/2 second-order coset representatives. Thus a first coset 
representative is produced by combining rows 2(X 1 *X 2 ), 
2(X 2 *X 3 ), . . . and 2(X m _!*X m ); the remaining coset repre- 
sentatives are produced by identifying all possible permu- 
tations of the subscripts of the row identifiers X v X^, X 3 , . . . 
X m in the expression 2(X a *X 2 )®2(X 2 *X 3 )® . . . 
©2(X M _ 1 *X JB ) which yield distinct results. 

These coset representatives are each combined by 
symbol-wise addition modulo 4 with all code words in the 
2 m -symbol base code to produce a total of m!(4 m+1 )/2 
possible code words, constituting a set of m!/2 cosets of the 
base code. 

Similar codes for other code word lengths and/or for 
coding with different numbers of symbols (e.g. octary) can 
be defined in like manner. Thus, where the number of 
symbol values is 2', the base code comprises all linear 
combinations of rows of the generator matrix, that is com- 
binations of the form 

where a 1; a^, ag, . . . a m and a' each take values of 0, 1, ... , 
2'-l and © indicates symbol-wise addition modulo 2 y . This 40 
base code is linear, i.e. the result of a symbol-wise addition 
modulo 2 ; on any two of its code words is another code word 
in the code. The corresponding generator array is created by 
combining the 2 m -symbol generator matrix with additional 
rows, which are derived by symbol- wise multiplication of 45 
all possible pairs of the rows X ly X 2 , X 3 , . . . X m and then 
further multiplication by 2 y ~ 1 . These additional rows are 
combined by symbol-wise addition modulo 2 y to produce 
m!/2 second-order coset representatives. Thus a first coset 
representative is produced by combining rows 2 /-1 (X 1 *X 2 ), 
2 M (X 2 *X 3 ), ... and 2 M (X m _ 1 *X m ); the remaining coset 
representatives are produced by identifying all possible 
permutations of the subscripts of the row identifiers X lt X^ t 
X 3 , . . . X m in the expression 2 / " 1 (X 1 *X 2 )©2 / - 1 
(X 2 *X 3 )S . . . ©2 /_1 (X m _ 1 *X m ) which yield distinct results. 

These coset representatives are each combined by 
symbol-wise addition modulo l! with all code words in the 
2 m -symbol base code to produce a total of m!(2 jCm+1) )/2 
possible code words, constituting a set of m!/2 cosets of the 
base code. 

In discussing the error detecting and error correcting 
properties of quaternary, octary and higher-order codes, the 
minimum Hamming distance alone does not fully describe 
the extent of possible errors affecting a code word. This is 
because an error in a single symbol can involve a transition 
through up to 2 /_1 adjacent symbol values (where the values 
2/-1 and 0 are considered to be adjacent). The practical 



50 



55 



effect may be as serious as, or more serious than, two or 
more errors in different symbols each involving only a single 
transition to an adjacent symbol value. For example, a 
change of a single octary symbol from value 6 to value 1 (a 
transition through 3 adjacent values) may be more serious 
than a change of two octary symbols, one from 6 to 7 and 
another from 7 to 0 (both being transitions to an adjacent 
value). Therefore an additional useful measure is the mini- 
mum Lee distance, defined as the minimum number of 
symbol value transitions required to transform a valid code 
word into some other valid code word. (The minimum Lee 
distance is not discussed above for binary codes, because in 
the binary case the Hamming and Lee distances are 
identical.) 

The code words defined above for the quaternary, octary 
and higher-order cases have the following advantageous 
combination of properties: 

the code's minimum Hamming distance is 2 m ~ 2 ; an error- 
correcting decoder based on Hamming distance can 
detect up to 2 m_2 -l errors arbitrarily affecting the 
symbols constituting a code word, and also correct up 
to 2 m-3 -l such errors; 
the code's minimum Lee distance is 2 7 "" 1 ; an error- 
correcting decoder based on Lee distance can detect 
errors involving up to 2 m_1 -l symbol- value transitions 
(affecting up to the same number of code word 
symbols), and correct errors involving up to 2 m_2 -l 
such transitions; in fact a minimum distance decoder 
based on Lee distance also corrects certain additional 
patterns of errors arbitrarily affecting the symbols con- 
stituting a code word (i.e. it possesses some of the 
properties of a decoder based on Hamming distance); 
a multi-carrier signal modulated in accordance with these 

code words has a PMEPR of no more than 3 dB; 
a decoder based on either Hamming distance or Lee 
distance for use with such a code can operate using 
analytical techniques and can be implemented, for 
example, using combinatorial logic at least in part; 
there are sufficient code words available to enable an 

acceptable data transfer rate to be attained. 
In practice these codes are used to transfer information 
represented by data bits over a path including the encoder 
and decoder. Although for binary codes the assignment of 
data bits to row coefficients of the generator matrix can be, 
as noted above, arbitrary, in all other cases (quaternary, 
octary, etc.) this transfer can be assisted by judicious assign- 
ment of data bits to row coefficients. 

In the case of a quaternary code, pairs of data bits 
(selected as described below) are each mapped to row 
coefficients a^ a^, a 3 , . . . a„ and a 1 as follows: 

TABLE 5 



Data bits 


Coefficient 


00 


0 


01 


1 


11 


2 


10 


3 



60 



65 



It can be seen that this mapping is an example of a Gray 
code, which has the advantageous property that an error 
involving a change between two adjacent coefficient values 
(3 and 0 being considered as adjacent) causes only a single 
data bit to be in error. Thus a relatively minor error in a code 
word is less likely to produce a large error in the data word 
at the end of the transfer. 



04/17/2004, EAST Version: 1.4.1 



US 6,373 : 

31 

As in the binary case, it is preferable to use a number of 
code words which is an integer power of two, thereby 
simplifying division of the incoming data stream into blocks 
for encoding, and furthermore to select code words in groups 
having the same coset representative, to facilitate implemen- 5 
tation. To this end the number of coset representatives 
chosen is the highest power of two (2 1 ) that is no greater than 
the number of coset representatives available. Then r of the 
data bits are used to select a coset representative, and 
successive pairs of the remaining data bits are converted to 10 
quaternary symbols a x , a 2 , . . . a m and a' according to the 
Gray code in Table 5 above. 

FIG. 12 shows a 13B/16Q encoder 160 analogous to the 
encoder 100 of FIG. 9, but for encoding 13-bit data words 
as 16-quaternary symbol code words in accordance with this 15 
invention, using code words derived as described above 
from the generator array 



(0000 


0000 


1111 


mi) 


<*j) 


(0000 


1111 


0000 


1111) 


<*2) 


(00 u 


001J 


0011 


0011) 


(X 3 ) 


(0101 


0101 


0101 
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(A4) 


(1111 


1111 


1111 


1111) 


(*') 


(0000 


0000 


0000 


2222) 


2(X,*X 2 ) 


(0000 


0000 


0022 


0022) 


2(X l *X 3 ) 


(0000 


0000 


0202 


0202) 


2{X X *X A ) 


(0000 


0022 


0000 


0022) 


2[X 2 *X 3 ) 


(0000 


0202 


0000 


0202) 


2(X 2 *X 4 ) 


(0002 


0002 


0002 


0002) 


2(Xi*X 4 ) 



25 



In particular the encoder uses the eight cosets having the 
coset representatives 



(0002 


0020 


0002 


2202) 


2(X l *X 2 )Q2{X 2 * 


x 3 )e2(* 3 


35 

*X 4 ) 


(0002 


0200 


0002 


2022) 


2(X l *X 2 )e2(X 2 * 


x 4 )e2iXi 


*X 4 ) 


(0000 


0220 


0022 


0202) 


2(X l *X 3 )®2(X 2 * 


Xi)e2tx 2 


*X 4 ) 


(0002 


0200 


0020 


0222) 


2(X,*X 3 )e2(X 3 * 


x 4 )<B2(x 2 


*X 4 ) 


(0000 


0220 


0202 


0022) 


2(X l *X 4 )®2(X 2 * 


* 4 )©2(X 2 


*X 3 ) 40 


(0002 


0020 


0200 


0222) 


2{X l *X A )B2{X 3 * 


X 4 )B2{X 2 


*X 3 ) 


(0002 


0002 


0020 


2202) 


2(X i *X 7 )®2(Xi* 


* 3 )©2(X 3 


*X 4 ) 


(0002 


0002 


0200 


2022) 


2(X 1 *X 2 )e2(X 1 * 


*4)©2(X 3 


*X 4 ) 



Thus 2 r has the value eight, so r is equal to 3. 

Referring to FIG. 12, the 13B/16Q encoder 160 includes 
five Gray coders 162 to 170 for encoding successive pairs of 
the ten least significant bits of the input data word into 
corresponding quaternary symbols, in accordance with 50 
Table 5 above. Each coder 162 to 170 supplies its output 
quaternary symbol to a respective one of five multipliers 172 
to 180, for multiplying that symbol by respective ones of the 
16-symbol generator matrix rows X 4 , X 3 , Xj and X'. The 
three most significant bits of the input data word are supplied 55 
to a selector 182 which selects one of the eight possible coset 
representative values in accordance with the combined value 
of those three bits. The outputs of the five multipliers 172 to 
180 and of the selector 182 are combined by a symbol-wise 
modulo-4 adder 184 to produce the 16-symbol quaternary 60 
code word output by the encoder 160. 

FIG. 13 shows a 16Q/13B minimum, distance decoder 200 
based on Lee distance for decoding 16-symbol quaternary 
code words from the encoder 160 to 13-bit data words. The 
method used in the decoder 200 is based on the "Supercode" 65 
decoding method referred to above in connection with the 
decoder 140 of FIG. 11, but with a significant modification 
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to enable its use with quaternary (and higher-order) codes. 
For each coset representative in turn that coset representa- 
tive is subtracted symbol-wise (modulo 4, for quaternary) 
from the received code word, and the result is decoded as an 
element of the quaternary base code of length sixteen, at the 
same time assigning a "score" to that coset representative as 
described below. The decision as to which code word to 
decode to is made by selecting a coset whose representative 
attains the highest score (arbitrarily if more than one attains 
that score), and within that coset selecting the element of the 
base code previously identified for that coset. 

The decoding of an element of the base code uses the fast 
Hadamard transform. In contrast to the decoder 140, 
however, the decoder 200 applies the fast Hadamard trans- 
form to multiple input vectors whose symbols are not 
constrained to have the value +1 or -1 (as was the case in 
the decoder 140), in order collectively to decode the result 
of the above subtraction as if it were a received code word 
of the base code. For each input vector a transform coeffi- 
cient of maximum magnitude is selected. According to the 
sign of that transform coefficient, one or other of an asso- 
ciated pair of elements of the first-order Reed-Muller code 
RM(l,m) could then be output. Although this is analogous to 
the procedure for decoding a received code word of RM(1, 
m), the output for each input vector would not represent a 
decoding of that vector as an element of RM(l,m). 

Referring to FIG. 13 the decoder 200 receives the 
16-symbol quaternary code words at a modulo-4 subtractor 
202, which subtracts from the code word, symbol-wise 
modulo four, a value in a coset store 204. This coset store 
contains the values of the same eight coset representatives 
from which selection is made by the selector 182 in the 
encoder 160. The store 204 also has a three-bit counter 206 
which cycles through all eight possible three-bit values 
successively in response to an increment signal at an input 
208, so that each of the eight coset representatives in turn is 
provided to the subtractor 202. 

The 16-symbol quaternary word S output by the subtrac- 
tor 202 for the current coset representative is provided to a 
vector generator 210, which for that word S supplies a 
sequence of thirty-two 16-symbol input vectors D to a Green 
machine 212 for performing the fast Hadamard transform on 
each of those vectors. To this end, the generator 210 assigns 
values of 0 or 1 to each of five coefficients b 1; b 2 , b 3 , b 4 and 
b 1 , so that the coefficients taken together range through all 
thirty-two possible combinations of values. For each such 
combination, a value B is derived according to the relation 

and then a vector D is derived from 

D=\Tr-((S-B)mod 4)(-r 

In this expression the subtractions involving values such as 
X', S and B are performed symbol-wise, as is the modulus 
operation indicated by | . . . |. 

The output of the Green machine 212 comprises, for each 
vector D, sixteen transform components labelled E 0 to E 15 . 
These components are supplied to a comparator 214 which 
compares the absolute values of the transform components, 
and identifies the largest absolute value |Ej and the corre- 
sponding value of z (from 0 to 15); if more than one 
component has the same largest absolute value, one of them 
is chosen arbitrarily. The comparator 214 then outputs the 
values z and E^ to a latch 216, which also receives the values 
of the coefficients b lt b 2 , b 3 , b 4 and b' from the generator 210 
and the corresponding three -bit index value supplied by the 
counter 206 in the store 204. 
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For the first coset representative and the first combination 
of values of the coefficients b 19 b 2 , b 3 , b 4 and b', the latch 216 
simply stores the values for z, E r , the five coefficients and the 
coset index which it receives. As thirty-one additional vec- 
tors D are derived by the generator 210 for the first coset 5 
representative, and as the counter 206 is incremented and 
thirty-two more vectors D are derived for each successive 
coset representative, respective values of z and E^ are 
provided to the latch 216, together with the respective values 
of the coefficients and the coset index. If such a transform 10 
component E z has a larger absolute value than the compo- 
nent currently latched, the latch 216 selects the newly 
received component, together with its associated values of z, 
the five coefficients b 3 , b 2 , b 3 , b 4 and b\ and the coset index. 

When the counter 206 has cycled through all eight pos- 15 
sible values, the contents of the latch 216 are used to 
assemble the decoded output data at 218. The value of z, the 
sign of E z and the coefficients b 1( b 2 , b 3 , b 4 and b' together 
determine the values of the ten least significant bits of the 
decoded data, as described below, while the coset index 20 
comprises the three most significant output bits. 

To determine the values of the ten least significant output 
bits, a coefficient c' is set to 1 if E z is negative and to 0 if it 
is positive, and the value of z is expressed as a four-bit 
binary number CjC 2 c 3 c 4 . Five decoded row coefficients a 4 , 25 
a 3 , a 2 , aj and a' are determined from the expressions 



a 2 -b 2 +2c 2 
o 3 =b 3 +2c 3 
a 4 =6 4 +2c 4 
a'-b'+2c f 

The decoded data bits can then be derived, in pairs, from 
these row coefficients by reference to the Gray code in Table 
5 above. 

For an octary code, triplets of data bits (selected as 
described below) are each mapped to row coefficients a A , a^ 
a 3 , . . . a m and a' as follows: 

TABLE 6 



Data bits 


Coefficient 


000 


0 


001 


1 


011 


2 


010 


3 


110 


4 


111 


5 


101 


6 


100 


7 



40 



45 



As before the number of coset representatives chosen is 
the highest power of two (2^ that is no greater than the 
number of coset representatives available. Then r of the data 
bits are used to select a coset representative, and successive 
triplets of the remaining data bits are converted to octary 
symbols a 1; a^, a 3 , . . . a m and a' according to the octary Gray 
code in Table 6 above. 

FIG. 14 shows an 18B/160 encoder 230 analogous to the 
encoders 100 of FIG. 9 and 160 of FIG. 12, but for encoding 
18-bit data words as 16-octary symbol code words in 
accordance with this invention, using code words derived as 
described above from the generator array 
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In particular the 
coset representati 

(0040 0040 

(0004 0400 

(0000 0440 

(0004 0400 

(0000 0440 

(0004 0040 

(0004 0004 

(0004 0004 



encoder uses the eight cosets having the 
ves 
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0040 0444) 
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0400 4044) 
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4(Xj* 



X 2 )<&4(X l 
Xi)®4{X 2 
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* 3 )©4(X3 
X 4 )<&MXt 
* 4 )©4(X3 



*X 3 )©4(X 3 *X 4 ) 
*X 4 )©4(X 3 *JV 4 ) 
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*X 4 )©4(X 2 *X 3 ) 
*X 4 )©4(X 2 *X 3 ) 

*X 3 )©4(X 3 *X 4 ) 
*X4)©4(X3*,Y 4 ) 



30 Thus 2 r has the value eight, so r is equal to 3. 

Referring to FIG. 14, the 18B/160 encoder 230 includes 
five Gray coders 232 to 240 for encoding successive triplets 
of the fifteen least significant bits of the input data word into 
corresponding octary symbols, in accordance with Table 6 
35 above. Each coder 232 to 240 supplies its output octary 
symbol to a respective one of five multipliers 242 to 250, for 
multiplying that symbol by respective ones of the 16-symbol 
generator matrix rows X 4 , X3, X 1 and X'. The three most 
significant bits of the input data word are supplied to a 
selector 252 which selects one of the eight possible coset 
representative values in accordance with the combined value 
of those three bits. The outputs of the five multipliers 242 to 
250 and of the selector 252 are combined by a symbol-wise 
modulo-8 adder 254 to produce the 16-symbol octary code 
word output by the encoder 230. 

FIG. 15 shows a 160/18B minimum Lee distance decoder 
260 for decoding 16-symbol octary code words from the 
encoder 230 to 18-bit data words, corresponding to that 
50 shown in FIG. 13 for decoding quaternary code words. 
Referring to FIG. 15 the decoder 260 receives the 
1 6-symbol octary code words at a modulo-8 subtractor 262, 
which subtracts from the code word, symbol-wise modulo 
55 eight, a value in a coset store 264. This coset store contains 
the values of the same eight coset representatives from 
which selection is made by the selector 252 in the encoder 
230, and outputs each coset representative in turn under the 
control of a three-bit counter 266 responding to an increment 
60 signal at an input 268. 

The 16-symbol octary word S output by the subtractor 
262 for the current coset representative is provided to a 
vector generator 270, which for that word S supplies a 
sequence of 1024 16-symbol input vectors D to a Green 
65 machine 272. To this end, the generator 270 assigns values 
of 0, 1, 2 or 3 to each of five coefficients b l9 b 2 , b 3 , b 4 and 
b*, so that the coefficients taken together range through all 
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1024 possible combinations of values. For each such 
combination, a value B is derived according to the relation 

B«b 1 X 1 eb 2 X 2 &>^Bb^C 4 @b'X f 

and then a vector D is derived from 



D«\W-((S-B)mo& 8)\-2X' 

As before the subtractions involving values such as X', S and 
B are performed symbol-wise, as is the modulus operation 10 
indicated by | . . . |. 

The output of the Green machine 272 comprises, for each 
vector D, sixteen transform components labelled E 0 to E 15 . 
These components are supplied to a comparator 274 which 
identifies a largest absolute value |EJ and the corresponding 15 
value of z (from 0 to 15) and outputs these values to a latch 
276. This latch also receives the values of the coefficients b l9 
b 2 , b 3 , b 4 and b' from the generator 270 and the correspond- 
ing three-bit index value supplied by the counter 266 in the 
store 264. 20 

For the first coset representative and the first combination 
of values of the coefficients b lf b 2 , b 3 , b 4 and b', the latch 276 
simply stores the values for z, E^, the five coefficients and the 
coset index. As 1023 additional vectors D are derived by the 
generator 270 for the first coset representative, and as the 25 
counter 266 is incremented and 1024 more vectors D are 
derived for each successive coset representative, respective 
values of z and E x are provided to the latch 276, together 
with the respective values of the coefficients and the coset 
index. If a transform component E^ has a larger absolute 30 
value than the component currently latched, the latch 276 
selects the newly received component, together with its 
associated values of z, the five coefficients b lt b 2 , b 3 , b 4 and 
b', and the coset index. 

When the counter 266 has cycled through all eight pos- 35 
sible values, the contents of the latch 276 are used to 
assemble the decoded output data at 278. The value of z, the 
sign of E x and the coefficients b 1( b 2 , b 3 , b 4 and b* together 
determine the values of the fifteen least significant bits of the 
decoded data, as described below, while the coset index 40 
comprises the three most significant output bits. 

To determine the values of the fifteen least significant 
output bits, a coefficient c' is set to 1 if E^ is negative and to 
0 if it is positive, and the value of z is expressed as a four-bit 
binary number CiC2C 3 c 4 . Five decoded row coefficients a 4 , 45 
a 3> a 2> a i an d a ' a re determined from the expressions 



a 2 ~b 2 +4c 2 
<2 3 «*& 3 +4c 3 
a 4 =6 4 +4c 4 
a'»i>V4c' 



results are added symbol-wise modulo 2 ; to the selected 
coset representative. 

Decoders corresponding to those shown in FIGS. 13 and 
15 can be constructed by analogy for other lengths of code 
word and/or higher-order codes. In the general case of code 
words of length 2 m symbols, each symbol having 2* possible 
values, the subtractor subtracts each coset representative 
from the received code word modulo 2. Coefficients b lf 
b 2 , . . . , b m and b' can each have integer values between 0 
and 2 /-1 -l inclusive, and taken together they range through 
2<j-i)(m+i) possible combinations of values. The value B and 
the vector D are derived from 



50 



55 



The decoded data bits can then be derived, in triplets, from 
these row coefficients by reference to the Gray code in Table 
6 above. 

Encoders corresponding to those shown in FIGS. 12 and 
14 can be constructed by analogy for other lengths of code 60 
word (mM) and/or higher-order codes (j>3). As noted 
above, r data bits are used to select a coset representative, 
and successive groups of size j of the remaining data bits are 
converted to respective symbols a^ a^,, a 3 , . . . a m and a' 
according to a Gray code table. The symbols a 1? a^, a 3 , . . . 65 
a m and a' are multiplied by respective rows X if X 2 , X 3 , . . . 
X^ and X ( of the base code's generator matrix, and the 



B-bJ 1 Gb#& . . . ®bJC„®b7T 

where © indicates symbol-wise addition modulo 2f. In 
determining the values of the respective least significant 
output bits, the value of z is expressed as a binary number 
c i c 2 • ■ * c mt aQ d decoded row coefficients a 1( a 2) . . . , a^ and 
a* are determined from the expressions 

a 1 »6 1 +2 >_1 c 1 
a'-*'+2 / - 1 C 

Each such row coefficient determines the values of j bits of 
the decoded output data by reference to the Gray code table 
used for encoding. 

The decoder as described above is a 'hard decision* 
decoder, i.e. it takes integer values between 0 and 2^-1 
inclusive as input. The same decoder can also be used for 
'soft decision 5 decoding, wherein the symbols of the 
received code word are represented as real numbers between 

0 and 2', according to the magnitude of the received ana- 
logue signal. The value S in this case is obtained by 
subtraction of the (integer-valued) coset representative 
modulo 2'. The above procedure for decoding is followed to 
produce real-valued transform components, of which one 
with largest absolute value is selected. The rest of the 
decoding procedure is the same as described above. 

The decoders 200 and 260 shown in FIGS. 13 and 15 
provide the decoded output binary data directly. Although 
they subtract a coset representative from the received code 
word and then decode the result as an element of the relevant 
base code, they do not explicitly identify a code word in the 
base code as an intermediate step. It is possible to obtain the 
value of this base code word, using the value z of a transform 
component E z having the largest absolute value, over all 
values B and all stored coset representatives, to select one of 
a pair of elements of the first-order Reed-Muller code 
RM(l,m), together with the corresponding value B. The 
column z of the Sylvester-Hadamard matrix H of order 2 m 
(see below) is transposed to produce a vector C (the left- 
most column of the matrix being column 0); if E^ is positive, 
+ symbols in C are converted to 0 digits and - symbols to 

1 digits, and if negative they are converted to 1 and 0 digits 
respectively; C now comprises the required element of 
RM(l,m). The base code word selected by the decoder can 
be derived from the expression 

The decoded row coefficients a,, a^, . . . , a m and a' can be 
obtained from this code word and the rows X l7 Xz, . . . , X m 
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and X' of the generator matrix by applying Gaussian elimi- 
nation in known manner. The decoded code word corre- 
sponding to that output by the encoder can be obtained by 
adding this base code word symbol-wise modulo 2! to the 
coset representative identified by the latched coset index. 

The Sylvester-Hadamard matrix of order 2 is 
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Hl6 = 
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i.e. all the elements are identical, except for that at the lower 
right, which is the negation of the other elements. To derive 
the matrix for the next higher power of 2, that is H 4 , the 
existing matrix (HJ is reproduced to its right and beneath 
itself, and its negation (-HJ is placed at the lower right: 
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This process is iterated to derive Sylvester-Hadamard 
matrixes for successively higher powers of two. For 
example, to derive H 8 , the matrix H 4 is reproduced to its 



Examples of these vectors and the corresponding transform 
components E are given below: 
For (b.^WHOAOAO): 

B-(0000 0000 0000 0000), D-(0,0,-l,l,l,-l,0,l,0,0 -1, 
1,1,-1,0,0), 

E=(1,-1,-1,9,-1,-7,1,-1,1,-1,-1 7 1,-1 7 1,1,-1). 
For (b^ba^b'MOAOAl): 

B=(llll 1111 1111 1111), D=(1,-1,0,0,0,0,1,0,1,-1,0,0, 

0,0,-1,-1), 

E^-1,5,1,3,1,3,-1,5,3,1,-3,-1,-3,-1,3, 1). 
For (b^b^jWHOAOXO): 

B=(0101 0101 0101 0101), D*(0,-1,-1,0,1,0,0,0,0,-1,- 
1,0,1,0,0,-1), 
For (^,0^3,0^X0,1,1,0,0): 

B=(0011 1122 0011 1122), D=(0,0,0,0,0,0,0,-1,0,0,0,0,0, 
0,0,0), 

E-(-l,l,l,-l,l,-l -1,1,-1,1,1 -1,1,-1 -1,1). 



right and beneath itself, and its negation (-H 4 ) is placed at 40 For ( D i, D 2, D 3, b 4, b, )=(°>l .1,0,1): 



the lower right. 

Consider, for example, the 13-bit data word 1 1011 0110 
0001 applied to the input of the 13B/16Q encoder 160. The 
corresponding row coefficients a 4 , a^, a^, a a and a' have 
values of 2, 1, 3, 0 and 1 respectively, giving a code word 
1320 0213 1320 0213 in the base code; the coset represen- 
tative selected will be 0002 0002 0020 2202. The output 
quaternary code word is therefore 1322 0211 1300 2011. 
Assume that during communication from the encoder 160 to 
the decoder 200 the code word is corrupted to a value of 
1322 0212 1300 2031 (corrupted symbols are in bold 
italics), so that the received code word has a Lee distance of 
3 from the transmitted code word. 

The subtractor 202 in the decoder subtracts each coset 
representative in turn from this received code word. In the 
case of the coset numbered 110, the result of this subtraction 
is 1320 0210 1320 0233, which is supplied to the vector 
generator 210. This generator produces thirty-two vectors D 
each of which is combined in the Green machine 212 to 
produce sixteen transform components E 0 to E 15 of a row 
vector E. Each such vector E is the result of multiplying the 
corresponding row vector D by the Sylvester-Hadamard 
matrix of order 16: 



45 



B=(1122 2233 1122 2233), D=(1,-1,1,-1,-1,1,-1,0,1,-1, 

1,-1-1,1,1,1), 
E«(13 -1,-3 -1,13,1,3,-3,-1,3,1,3,1,-3,-1). 
For (b^b^b^XO,!,!,!^): 

B=(0112 1223 0112 1223), D=(0,-1,0,-1,0,1,0,0,0,-1,0,- 

1,0,1,0,1), 

E=(-l,l,l,-l,-7,7,-l,l,-l,l,l,-l,l,-l,-l,l). 
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For (b 1 ,b 2 ,b 3 ,b 4 ,b')-(l,l,l,l,0): 

B=(0112 1223 1223 2330), D«(0,-1,0,-1,0,1A0,1,0,1, 

0,-1,0,1,0), 
£=(1^3,-1,-3,-1,5,13,-3,-1,3,1,-5,1,-3,-1). 
For ^,^3,^^X1,1,1,1,1): 

B=(1223 2330 2330 3001), D=(l ,0,1,0,-1,0,-1,1,0,1,0,1, 

0,-1,0-1), 

The transform coefficient of largest absolute value over 
this coset is E s =13, which occurs for B=(1122 2233 1122 
2233), i.e. z=5; this coefficient is also larger in absolute 
value than the transform coefficients produced for any of the 
other seven coset representatives in the decoders store 204. 
Column 5 of H 16 , transposed, is (+-+ — +-++-+ — +-+). As 
E 5 is positive, + symbols are replaced by 0 digits and - 
symbols by 1 digits, so C=(0101 1010 0101 1010). Hence 
the base code word selected by the decoder is 2C®B=(1320 
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0213 1320 0213), which can be seen to be the base code 
word incorporated in the output of the encoder 160, With the 
communication errors corrected. 

Minimum distance decoders based on Lee distance for 
quaternary, octary and higher-order codes, such as the 
decoders 200 and 260 shown in FIGS. 13 and 15, may be 
implemented in special -purpose hardware; in this case, and 
as described previously in relation to the decoder 140 of 
FIG. 11, some or all of the fast Hadamard transforms can be 
calculated in parallel, for example by providing a Green 
machine for each one. Alternatively a general-purpose pro- 
cessor may be used, executing an appropriate software 
procedure incorporating functionality illustrated by the pro- 
gram presented below and expressed in the C programming 
language. This program takes as input: 

a received code word of length determined by the param- 
eter M (»m in the above description), each symbol 



40 



having a number of values determined by the parameter 
J (=j in the above description); 
a list of coset representatives. 
It produces as output: 

the position in that list of a coset representative for which 
a transform coefficient of largest absolute value was 
found; 

the base code word identified from, in part, the position z 

of that transform coefficient; and 
the values of the corresponding row coefficients a a , 

a^ . . . , a m and a'. 

Although primarily intended for use with quaternary or 
higher-order codes, this program is also usable with binary 
codes. 



/* Decoding of an arbitrary set of coscta of W(M) using fast 

Hadamard transforms. W(M) is the 2 J-phasc linear code having the 
same forma] generator matrix as RM(1,M). 

7 

^define M 4 /• Integer M >= 2 7 

#denne J 2 /• Integer J >= 1 

Program assumes (J-l) (M+l) < 31 */ 
#defrae NON_BINARY (3-3) /• Flag for all cases except binary (J=l) •/ 

#define L (1 « M) /• Sequence length L « 2 M 7 

#define P (1 « J) /* Allowed sequence symbols are 0 f l ( 2,...P-l */ 

^define ABS(x) ((x) > 0 ? (x) : (-(x))) 

int genlM+lJlj; /* Generator matrix for W(M) "/ 

#indude <stdio.h> 

mainO 

{ 

int coseL_rep[L], recvdL_seq[L], seq[Lj fhLJnput[L]; 

int coefijM+1 J /* coeff of row i of W(M) V 

int a[M+ll b[M+l]; /" coeff[i] = (P/2)a[i] + b[i] V 

int b_sum[L]; /" codeword of W(M) given by the b[i] */ 

int score, posn, rep_id, best_score, best_posn, best_rep_id; 

int i, j, k, best_Jc; 

/* Initialise gen[J] */ 

initialiseO; 

/* Input received codeword */ 

printf(*lnput received codeword as a length %d sequence L); 
printf("of symbols,\neach of which takes values from 0 to %d 

. 

printf("inchisive\n\n"); 
if (input(recvd_seq, L) — EOF) 
exitO; 

/* Decode each coset rep of W(M) in turn 7 
rep_Jd - -1; 
best_score - 0; 

printf("\nlnput a list of coset reps in the same form, "); 
printf^crminatcd by EOF\n"); 
while(input(coset_rep, L) !» EOF) 

rcp_id ++; 

I* Subtract coset rep from received sequence componentwise 

to give sequence seq[] 7 
for 0" - 0; j <: L; j++) 

seq[j] » (recvd_seq[j] - coset_rep[j] + P) % P; 
/• Set b_sum|> sum_{i=0} "M b[i] gen[i][] for all b[i] 

satisfying 0 <= b[i] < P/2 •/ 
for (k = 0; k < (1 « (J-l)*(M+3)); k++) 
{ 

extiact_symbols(k, J-l, M+l, b); 
sum(b_sum, b); 

/• Calculate input vector fht_input[] to fast 

Hadamard transform. In the binary case multiply 
by 2 to give integer values 7 

for (j - 0; j < U j++) 
fht_input(j] - 

#if NON_BINARY 

ABS(P/2 -((seq[j] - b_sum[j] + P) % P)) - P/4; 

#else 
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ABS(P -2'((scq[j] - b_sum[j) + P) % P)) - P/2; 

#endif 

/* Apply fast Hadamard transform to fht_input(] "/ 
fast__hadamard_transfonnC£ht__input 1 &scorc, &posn); 
if (ABS(scorc) > ABS(best_score)) 
{ 

bcst_scorc = score; 
bcst__posn = posn; 
bcst_Jc - k; 
bcst__rcp_id - rcp_id; 

} 

} 

} 

/* Calculate row coefficients for decoded codeword in W(M) */ 
40] - (besL_score < 0 7 1 : 0); 

/* coeff of all- ones codeword in RM(1,M) 7 
extract__symbols(besL_posn ) 1, M, &a[lD; 

r the other M first-order coeffs in RM(1,M) 7 
extract_symbols(best_Jc, J-l, M+l, b); 

/" recover best coeffis b[] *f 
for (i - 0; i < M+l; i-H-) 

coefi^i] - (P/2) * a[i] + b[i]; 
/* Calculate best codeword in W(M) 7 
sum(seq, coeff); 

/* Output number of best coset rep (the coeet reps are indexed by 

0, 1, . . . rep_number) 7 
printf("\nNumbering coset reps as 0,1,... the decoded coset rep "); 

printf("number is %2d.\n", best iep_id); 

/* For this coset of W(M), output best codeword and its row 

coefficients */ 

printf("Foi this coset of the base code W(M), best codeword is\n"); 
for (] - 0; j < L; j++) 

print£("%2d" secjj]); 
printf(",\nand corresponding row coeff s of generator matrix "); 
printf("for W(M) are\n")i 
for (i - 1; i < M+l; i++) 

printf(" for rowX%d: %2d\n", i, coeff[iD; 
printf(" for row X': %2d\n\n", coefflOD; 

} 

initialiseO 
{ 

int a[Ml i, j; 

/* Set gen(il] to be i-th row of generator matrix for W(M) 

(where row 0 corresponds to X' and row i corresponds to X_i) 7 
for (j = 0; j < U j ++) 
{ 

gen[0ljj=l; 

extract_symbols(j, 1, M, a); 
for (i = 1; i <= M; i++) 
gen[i][j] - 

} 

} 

int inpuu^seq, length) 
int *seq, length; 

{ 

int j; 

if (scanf( t< %d" ( &seq[0]) — EOF) /* Check for end of file 7 

return(EOF); 
for (j - 1; j < length; j-M-) 

scanf(" %d", %seq[jD; 

return(O); 

} 

/* Set seq[] to sum_{i-0} {M} b[i] gen[i][] mod P */ 
sum(seq, b) 

int seojLl b[M+l]; 
{ 

mt i, j; 
/• Initialise 7 

for G • 0; j < U 

seq[j] - 0; 
for (i = 0; i <= M; i++) 
for (j - 0; j < L; j++) 

seq[j] - (seqlj] + b[i] • genlilj]) % P; 

} 

/• Store the s least sig (2 u)-ary symbols of the integeT c, in order 

from most sig to least sig, in seq[0], seq[l], seq[s-l] V 
int extract_symbols(c t u t s, seq) 
int c, u, s, *seq; 
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unsigned int mask; 
int i; 

mask » (1 « u) -1; /* mask for u least sig bit positions 7 

for (i = s-1; i >= 0; i — , c >>= u) 
scq[i] » c & mask; 

} 

fast_hadamard_transform(seq, scorc_ptr, posn_ptx) 

int scqlH *score_ptr, *posn_ptr; 

{ 

int vecfM+lflL]; 

int best_score, bcst_posn; 

int p, i, j, k, t; 

/* Copy scq[] to vecfOl] •/ 

for 0 - 0; j < L; 

vec(0lj] o seq[jj 
for (i » 1; i M; i++) 
{ 

/' Set vecfij] to the result of applying the first i matrices 
of fast Hadatnard transform to seq[] */ 
t = 1 « (i-1); 

for (j - 0; j < (1 « (M-Q); j++) 

for (k - 0; k < t; k++) 

{ 

p-k + j*(l«0; 

vecfijp] - vec[i-llp] + vec[i-l][p+t]; 
vec[i]Ip+t] - vec[i-ljp] - vecfi-ljp+t]; 

} 

} 

/* Find entry of vec[M][] of maximal magnitude */ 

best__soore - 0; 

for (j - 0; j < L; j++) 

{ 

if (ABS(vec(Mlj]) > ABS(best_soorc)) 
{ 

best_posn = j; 
best_score = vec[Mlj]; 

} 

} 

*posn__ptr = best_posQ; 
*score_ptr = best_score; 



The 2od-order binary RM decoders 82 and 122 shown in 
FIGS. 8 and 10 may be implemented in special-purpose 40 
hardware, or by using a general-purpose processor executing 
an appropriate software procedure such as that presented 
below and expressed in the C programming language. 



/* Decoder for 2ad order binary Recd-Mullcr code RM(2,M) of length 2*M, 
using majority logic Reed decoding algorithm. 

V 

#define M 4 

tfdefine L (1 « M) /• Sequence length L - 2*M 7 

int gen[M+llL]; 

#include <stdio.h> 

mainO 

{ 

int i, k; 

int b[M+ll seq[L]; 
int prodfL], prod_comp[L]; 
/* Initialise genffl] 7 
initial iseO; 

/* Input received codeword •/ 

printf("\nlnput received codeword as a length %d sequence*',L); 
printf(" of symbols,\neach of which takes values 0 oi l\n"); 
if (input(seq f L) — EOF) 
exitO; 

printf("\nDecoded codeword is given by the following terms\n"); 
/* Determine coefficients of order 2, 1 and 0 7 
for(k = 2; k >= 1; k— ) 
{ 

/* Initialise set b[] which indexes combinations 7 
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for (i = 1; i < k; i++) 

b[i]-ii 
b[k] - k-1; 

/* Select all combinations of k objects from 

{1,V.,M} 7 
while(aext_comb(b, M, k) -» 1) 
{ 

r Calculate product of rows indexed by 
entries of b[], and those by their 
complement 7 
product(prod, prod„comp, b, k); 
/• Use majority logic decoding to determine if 
coefficient corresponding to b[] is set. If 
so, print coefficient and update sequence */ 
maj ]ogic(prod, prod_comp, b, seq, k); 

} 

} 

!* Treat k « 0 separately */ 
product(prod, prod_comp, b, 0); 
maj_logic(prod, prod__comp, b, seq, 0); 

} 

initialiseO 
{ 

int a[M], i, j; 

/* Set gen[il] to be i-th row of generator matrix for RM(1,M) 

(where row 0 corresponds to X' and row i corresponds to X_i) 7 
for (j - 0; j < L; j ++) 
{ 

gen[0lj]-l; 

extract_symbols(j, 1, M, a); 
for (i - 1; i <- M; i++) 

gen[i]0] - a[i-lj; 

, 1 

/* Store the s least sig (2 u)-ary symbols of the integer c, in order 
from most sig to least sig, in seq[0], seq[l], seq[s-l] 

7 

int extract_symbo]s(c, u, s, seq) 
int c, u, s, ""seq; 

{ 

' unsigned int mask; 
int i; 

mask » (1 « u) -1; /* mask for u least sig bit positions 7 
for (i » s-l; i >= 0; i — , c >>«= u) 
seq[i] - c & mask; 

} 

product(prod, procLcomp, b, k) 
int prod[L], prod_comp[L]; 
int b[M+l], k; 
{ 

int i, u; 
int a[M+l]; 

/* Convert combination represented by b(] into 1-0 form a[] 7 
for (i = 1; i <= M; i++) 

a[i]oO, 
for (i o 1; i <= k; i++) 

a[b[il]-l; 

/* Initialise product of rows indexed by a[], and product of 

rows indexed by complement of a[] */ 
for (u - 0; u < L; u++) 

prod{u] - prod comp[u] - 1; 

for (i - 1; i <- M; i++) 

if(a[i]~ 1) 

{ 

for (u - 0; u < 1^ u-m-) 

prod[u] "= gen[iju]; 

} 

else 
{ 

for (u = 0; u < L; u++) 

prod_comp[u] *= gen[iju]; 

} 

} 

maj_logic(prod, prod_comp, b, seq, k) 
int prod(L], prod_comp[L], b[M4l \ 
int seaJL], k; 
{ 

int zeroes, ones, sum; 
int i, u, v; 
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zeroes » ones ™ 0; 
for (u - 0; u < L; u++) 

if (prod[L-l-u] — 1) 

{ 

sum d 0; 

for (v ■ 0; v < L; v++) 

if (prod_comp[v] = 1) 

sum += seq[v-ut 

if (sum % 2 — 0) 
zeroes ++; 

else 

ones ++; 

} 

if (ones >= zeroes) 
{ 

if (k > 0) 
{ 

printf (" X%d'\ b[l]); 
for (i - 2; i <» k; i++) 

printf(" * X%d", b(i]); 
printf("\n"); 

} 

else 

printf(" XV); 
for (u - 0; u < L; u++) 

eeq[u] - (seq[u] + prod[u]) % 2; 

} 

} 

/* Lexicographically next k-subset from M-element set */ 
nexL_comb(b, n, r) 
int b[M+l], n, r; 
{ 

int i, h; 

for (i = r; i >= 1 ; i — ) 
{ 

if (b[i] < n-r+i) 
{ 

b[i]-«i]+l; 

for (h = t+1; h <= r, h++) 
b[h] = b[h-l] + 1; 

return (1); 

) 

} 

retum(O); 

> 

int input(seq, length) 
int *seq, length; 

{ 

int j; 

if (scanf("%d", %seq[0]) = EOF) /* Check for end of file V 

return(EOF); 
for (j = 1; j < length; j++) 

scanf(" %d", %seq[j]); 

retum(O); 

} 

output(seq, length) 
int *seq, length; 

{ 

int j; 

for (j - 0; j < length; j++) 

printf(" %d", seojj]); 
printf( u \n"); 



In the case of quaternary codes of length 2 m symbols, 
decoding can also be accomplished by using a binary 
decoder for a code of length 2 m+1 bits; this may be desirable 
for ease of implementation (by avoiding the need for cir- 
cuitry to operate on quaternary signals) and/or for speed of 
operation. Such decoding will be illustrated by reference to 
the 16-symbol quaternary code word 1322 0211 1300 2011 
(already used as an example above), derived from a base 
code word 1320 0213 1320 0213 (row coefficients a 4 , a3, 
a 1 and a' with values of 2, 1 , 3, 0 and 1 respectively) and a 
coset representative 0002 0002 0020 2202. This code word, 
representing the binary data word 1 1011 0110 0001, is 



assumed to be received with errors as 1322 0212 1300 2031. 
Each symbol of the received code word is first converted to 
a pair of binary symbols, according to the Gray code in Table 
5, and these binary symbols are assigned to respective halves 

60 of a 32-bit word. Thus, for the first quaternary symbol (1), 
the corresponding pair of binary symbols is 01; the first 
symbol (0) of this pair comprises the first symbol of the 
32-bit word, and the second symbol (1) of the pair comprises 
the seventeenth symbol of that word. The second quaternary 

65 symbol is 3, so the corresponding binary pair is 10; the first 
of these symbols provides the second symbol of the 32-bit 
word, and the second symbol (0) of the pair provides its 
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eighteenth symbol. Continuing this procedure for the rest of ten least significant output bits are derived, in pairs, from the 
the quaternary code word produces the 32-bit word decoded row coefficients a 1? a^, a 4 and a' by reference to 

the Gray code in Table 5 above. 

It is noted that decoding quaternary code words contain- 
This binary word is supplied to a 32-bit 2nd-order binary 5 m S 2 ™ s)™ 0015 ^8 a b ™vy decoder for a code of length 
RM decoder, similar to the decoder 122 in FIG. 10 but 2 m+1 bits in this manner involves reference to the Gray code 
extended to handle additional bits; this produces two sets of table for ^ ^ncl purposes: one purpose is to assist the 
binary output signals which together comprise the output transfer of data, with detection and correction of errors, over 
value of the RM decoder: six first-order signals representing me complete communication path, as already described with 
coefficients u 1? a,, u 3 , u 4 , u 5 and u', and ten second-order 10 reference to the encoders 160 and 230 of FIGS. 12 and 14; 
signals representing coefficients u 12 , u 13 , u 14 , u 15 , u^, u^, second Purpose is to permit the use of binary decoder 

"25, u 34 , u 35 and u 45- In thc casc of 32-bit word ° oted circuitrv . thereby avoiding the expense and effort of imple- 
above, the output signals of the RM decoder would have the menting circuitry for operating directly on quaternary sig- 
values: na ^ s m g ener al, enhancing speed of operation. 

15 To apply a binary decoder in the general case of a 
quaternary code word containing 2 m symbols, the code word 
is converted, symbol by symbol, to a corresponding 2 m+a -bit 
binary code word as described above and using the Gray 
code in Table 5, respective symbols in each pair of Gray 
20 code binary symbols being assigned to respective halves of 
that binary code word. This binary code word is decoded as 
an element of the second-order Reed-Muller code of length 
2 m * 1 > RM(2, m+1), using a binary decoder to give coeffi- 

m , 1 , cients u 1f u,, . . . u m .,, u\ u,,, u 13 , . . . u--.,. These 

These first- and second-order signals are supplied to a a coeffidents can ^ be ttulsfi)rined int0 ^ re ^ uired qua . 

binary to quaternary coefficient converter which provides , e Qut t coefficients according to the relationships 

row coefficients a 1? a 2 , ag, a 4 and a' corresponding to the 

input quaternary data word according to the following q'-Iu'm^ 

relationships: 

fl / =(2« /+1 +w 1 J+1 (l+2u 1 ))mod 4 for 1 %i^m 

30 

ai = (2*2 +«I2(1 +2«i)) mod4 [=01 . N , . c . 

J ff,X«, + i r *M+«vM»i* + i)™d2forl£i<*g/» 

a 2 - +«i 3 (l +2ui)) mod4 [=3] 

a 3 = (2u 4 +«i 4 (l +2«i)) mod4 [= 1] The values a lF a^, , . . a m , and a' give 2(m+l) decoded data 

04 = (2w 5 +«i5(i + 2« 1 )) mod4 [ = 2] bits via the Gray code in Table 5; the remaining bits are 

a'=2U+ui [=1] 35 recovered by reference to a look-up table of the coset 

representatives chosen for use, according to the values of the 
coefficients a 12 , a 13 , . . . a^,.^. 
The values in brackets [ ] indicate the results of these The assignment of respective symbols in each pair of 
relationships for the specific values of the first- and second- Gray code binary symbols to respective halves of the 
order signals given above. It can be seen that they produce 4 q 2 m+1 -bit binary code word, as described above, is not the 
the appropriate quaternary value 2, 1, 3, 0 and 1 for the five only possibility. Alternative assignments of the binary sym- 
decoded row coefficients a 4 , a 3 , a 2 , a j, and a 1 . bols to respective locations in the binary code word may be 

The coefficient converter also produces six signals repre- used, with corresponding alteration of the relationships for 
senting coefficients a 12 , a 13 , a 14 , a^ 3 , a^ and a 34 according transforming the coefficients u^u^ . , . u,^, u', u 12 , u 13 , . . . 
to the following relationships: 45 u mj/ri+ l into quaternary output coefficients. For example, the 

quaternary code word to be decoded could be expanded by 
ai2 = («23 +■ "ii«i3) mod 2 [ = l] concatenating the pairs of Gray code binary symbols in the 

a 13 = (W24 -»-«]2"i4) mod2 [ = l] same order as the corresponding quaternary symbols in the 

a l4 = (u^ + u 12 «. s) mod 2 [ = 0] quaternary code word. Thus the quaternary word 1322 0212 

a* = <«» + «»«i4) mod 2 [ = 0] 50 1300 2031 would be converted to the binary word 

<*2A = («3S +«13«13) mod2 [ = 0] ono lln 001] oin ono 0000 n00 10Q1 



a* = («4S +■ «H«i3) mod 2 [ = 1] 



In this case the relationships for coefficient transformation 
would be 



These additional signals together identify which of the six 55 

additional rows 2(X 1 *X 2 ), 2(X 1 *X 3 ), 2(X 1 *X 4 ), 2(X^X 3 ), «'-2«^.i 
2(X,*X 4 ) and ^X^X^) respectively of the generator array 

were combined to form the coset representative from which M(2«^«^,Xi+2« m+ ,))inod 4 for 1 Zi^m 

the 16-symbol quaternary code word was produced. In the a^t S»,^«^mod 2 for i*i<ksm 
case of the data word 1 10110110 0001 the coset represen- 60 

tative in question is 2(X 1 *X 2 )©2(X 1 *X 3 )©2(X 3 *X 4 ), and it The substitution of binary processing for quaternary pro- 
can be seen that the relevant coefficients a 12 , a 13 and z^ A cessing is also possible with decoders having the form of the 
have values of 1 to indicate this. 16QA3B decoder 200 of FIG. 13. This may be 
The combination of the six coefficient signals a 12 , a 13 , a 14 , accomplished, in the case of the decoder 200, by modifying 
a 23> ^4 and a 34 can be used to derive the three most 65 it so that the received quaternary code word and the coset 
significant bits of the data word by means of a selector representative output by the coset store 204 are both con- 
similar to the selector 128 in the decoder 120 of FIG. 10. The verted to 32-bit code words, symbol by symbol and using the 
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Gray code in Table 5 as already described. The subtracter to a 2 M * J -bLt binary word. For determining the values of the 

202 is correspondingly modified for operation modulo 2 2(m+l) least significant output bits from row coefficients, 

rather than modulo 4. the value of z is expressed as a binary number c a c 2 . . . c,^, 

The 32-bit word S output by the subtracter 202 for the and me deco ded row coefficients a 1( . . . , a^ and a' are 

current coset representative is provided to the vector gen- s determined from ^ expr essions 
erator 210, which is also modified to supply a sequence of 
sixteen 32-bit input vectors D to the Green machine 212 for 

performing the fast Hadamard transform on each of those a^ic^+b^i+ic ,))mod 4 for i<-i<«m 
vectors. To this end, the modified generator 210 assigns 
values of 0 or 1 to each of four coefficients b it b 2 , b 3 and M , 10 
so that the coefficients taken together range through all 
sixteen possible combinations of values. For each such 

combination a quaternary 16-symbol value is derived ^ above of quat ernary codes have been 

accor ing to e re ation ^ expressed in terms of quaternary symbols having possible 

Q> l x 1 +b 2 x 2 +biX :s +bjc 4 )mod 4 values of 0, 1, 2 and 3. For actual transmission, for example 

by quadrature phase-shift keying (QPSK) modulation, these 

and this value is also converted to a 32-bit word B, in the values would correspond to four different values of phase 

same manner as the input quaternary code word and the snift Qne ssibility ^ phase shifts of ^ 90 °, 1SQ « and 

current coset representaUve It should be n^ ^ ^ ( ivalent to lex vectors j u _ x ^ _ { 

the above relation is similar to that for the value B in the case ).^. . ** . 

of the unmodified decoder 200, it is not identical; in par- respectively, where i-V-1). However, other variants are 

ticular the modified relation does not involve a coefficient b\ P ossible m whlch P hasc smfts are ' for sample, 45 , 

Each vector D is derived by bitwise addition modulo 2 of 135°, 225° and 315° (equivalent to complex vectors (l+i)/ 

the 32-bit word S and the respective 32-bit word B, and the V2, (-Ui)/V2, (-l-i)/V2 and (l-i)/V2 respectively). This 

values +1 and -1 are substituted for the binary values 0 and principle applies for codes with more symbol values (j>2), 

1 respectively in that vector D. Then the vector D is supplied and indeed the symbol values may be represented by any 2/ 

to the Green machine 212, modified to implement the fast distinguishable signal modulations; their amplitudes need 

Hadamard transform of length thirty -two. not be equal, and, in the case of phase -shift keying, the phase 

The output of the Green machine 212 comprises, for each 3Q differences between adjacent pairs need not all be equal, 
vector D, thirty- two transform components, which are sup- 
plied to the comparator 214 (modified to receive thirty-two The methods and apparatus described above for encoding 
components) and the latch 216, for selection, over all the and decoding COFDM signals provide a highly advanta- 
values of B for each of the coset representatives, of the gcous combination of error detection and correction 
transform component E x with the largest absolute value, and 35 properties, low PMEPR, convenience of implementation 
the corresponding values of z (from 0 to 31), the four and useful rate of data transfer. In particular where long code 
coefficients b lf b 2 , b 3 and b 4 , and the coset index. The words (; e many carr i ers ^ a multi-carrier system) are 
selected value of z, the sign of and the coefficients b lf b 2 , requ i re d to provide a high data rate these properties can be 
b 3 and b 4 together determine the values of the ten least obtained without undertaking an extensive (and probably 
significant bits of the decoded data as described below, 4Q ^ actical) cxhalistivc sca rch through all possible code 
while the coset index comprises the three most significant words Iq circumstances it may be desired to modify 

ou Jut bits. the selection of code words to alter the trade-offs between 

To determine the values of the ten least significant output 

bits, a coefficient c' is set to 1 if E, is negative and lo 0 if it ^ dlff6rent P«»Perties. 

ispositive^ndthev^eofzisexpressedasafive-bitbinary 4j For e j be desired t£) ^ avaikbu 

number c.c^c c 5 . Five decoded row coefficients a 4 , a„ a,, ^ on of code word in order t0 increase the data 

a, and a' are determined from the expressions f r it *u r *• u*u 

1 r transfer rate, at the expense of accepting a higher maximum 

fl 1 =(2c 2 +& 1 (i+2c 1 ))mod4 va l ue of PMEPR. Considering binary codes of length six- 
teen (m»4), for example, it is possible to construct every 

a?={icj-b 1 $+2ci))mo& 4 so possible coset of the base code whose coset representative is 

a 3 =(2c 4 +b,(i +2c 1 ))mod 4 a linear combination of the additional rows in the generator 

array, and to derive the maximum peak envelope power 

a 4 -(2c 5 +ii 4 (i+2c 1 ))mod4 (p£p) resu iting from use of code words within that coset. 

a'mic'+c Tb^ enables the cosets to be listed in ascending order of 

55 maximum PEP over each coset, as shown in Table 7. In this 

The decoded data bits can then be derived, in pairs, from table the first column contains six coefficients which con- 

these row coefficients by reference to the Gray code in Table stitute multipliers for the respective additional rows of the 

5- code's generator array which are added symbol-wise 

To use a modified decoder 200 in this way for the general mo dulo 2' to create the coset representative; thus, for 

case of a quaternary code word containing 2 m symbols, the examp i ej a i under the heading '12' indicates that the 

code word and the coset representative are converted sym- addilional row (x/XJ is multiplied by 1 for inclusion in 

bol by symbol, to a corresponding binary code word. Hie ^ co$ct cscntativC( whcreas a 0 it fa multi . 

binary coefficients b lf b , b taken together range ed Q effectivel omitted> ^ middle column 

trough 2 possible combinauons of values. The value B is * * ^ ntalive ^ the ^ cohmm 

derived by converting the value of the quaternary expression . , . * n „ n ' c , 

shows the maximum PEP resulting from use of any code 

(^+^2+ . . . +fcjrjmod 4 word in that coset. 
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TABLE 7 



54 



Binary length 16 
'Additional row' coefficients 



PEP over 



12 13 1 4 23 24 34 Coset representative coset (W) 



0001 0010 0001 1101 
0001 0001 0010 1101 

0000 0110 0011 0101 

0001 0100 0001 1011 
0001 0100 0010 0111 
0001 0010 0100 0111 

0000 0101 0110 0011 

0001 0001 0100 1011 

oooo oioi con looi 

0000 0011 0101 1001 
0000 0110 0101 0011 
0000 0011 0110 0101 



0001 0001 0001 1110 
0000 0011 0101 0110 

0000 0101 0011 0110 

0001 0100 0111 0010 
0001 0100 0100 1110 

0000 0110 0011 1010 

0001 0111 0001 1000 
0001 0111 0100 0010 
0000 0110 0101 1100 
0000 0101 0110 1100 

0000 0011 0110 1010 

0001 0010 0111 0100 
0001 0001 0111 1000 
0001 0111 0010 0100 
0001 0010 0010 1110 



0001 0100 0001 0100 
0001 0010 0001 0010 

0000 0110 0000 0110 

0001 0001 0100 0100 

0000 0101 0101 oooo 

0001 0100 0100 OO01 
0001 0001 0010 001 0 

0000 0011 0011 oooo 

0001 0010 0010 0001 
0000 0000 0110 0110 

0000 0110 0110 oooo 

0001 0111 0111 0001 
0000 0101 OOOO 1010 
0000 0011 OOOO 1100 
0000 0110 0000 1001 

0000 0000 0101 1010 

0001 0010 0100 1000 
0003 0111 0100 1101 

0000 0000 0011 1100 

0001 0100 0010 1000 
0001 0111 0010 1011 
OOOO OOOO 0110 1001 
0003 0100 0111 1101 
0003 0010 0111 1011 

0000 0110 0110 1111 

0001 0111 0111 1110 

0000 OOOO 0011 0011 

0001 0001 0111 0111 
OOOO 0011 0011 1111 
OOOO 0011 OOOO 0011 

0000 OOOO OOOO 1111 

0001 0111 0001 0111 
OOOO 0101 OOOO 0101 

0000 0101 0101 1111 

0001 0001 0001 0001 
OOOO OOOO 0101 0101 

oooo oooo oooo oooo 



31.59 
31.94 
31.95 
31.98 
31.98 
31.98 
31.98 
31.98 
31.99 
31.99 
32.00 
32.00 

-- ( 

49.82 
49.87 
49.98 
50.88 
51.10 
51.12 
51.65 
51.76 
51.81 
52.87 
52.90 
53.47 
53.56 
53.82 
53.99 

— ( 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 
64.00 

— ( 
98.95 
99.72 

101.43 
101.55 
105.60 
105.85 
106.22 
106.41 
106.69 
109.48 
109.75 
256.00 



10 



(1) 



(2) 



20 



25 



30 



35 



45 



interval between the time values for successive samples 
affects the accuracy of the result: a shorter time interval 
increases the accuracy, but at the expense of increasing the 
number of samples required for the entire waveform and 
thus increasing the amount of calculation. Typically the 
derivation of the PEP is repeated for successively shorter 
intervals until the values of PEP obtained are stable to some 
predetermined number of significant figures or decimal 
places. In the case of Table 7, the calculation has been 
repeated in this way until the second decimal place, after 
rounding, is stable. 

As explained above, the 384 code words associated with 
the first twelve coset representatives listed in Table 7, above 
dashed line (1), provide a PMEPR of no more than 3 dB, i.e. 
a factor of two (the mean envelope power for code words of 
length 2 m is 2 m ; hence for code words of length sixteen the 
mean envelope power is 16, so a PEP of up to 32 gives a 
PMEPR of up to two). 

If the code words used are those associated with the first 
27 coset representatives in the first half of Table 7, i.e. above 
dashed line (2), plus any five of the next 25 coset 
representatives, i.e. above dashed line (3), it becomes pos- 
sible to encode ten data bits while maintaining the minimum 
Hamming distance unchanged at 2 m ~ 2 (»4). This increases 
the number of data bits which can be represented with 16 -bit 
code words from eight to ten, while the maximum PMEPR 
increases from two to only four (rather than sixteen which 
would be the case if any code word were permitted). These 
extra data bits can be accommodated by the selector 112 in 
FIG. 9, which in this case can be modified to select from any 
of the chosen 32 (=2 S ) coset representatives, in accordance 
with the values of five most significant bits of a data word 
(instead of only three such bits). 

Tables similar to Table 7 can be prepared for other lengths 
of binary code word. For a code word length of 2 m , any code 
comprising a set of cosets from such a table has a minimum 
Hamming distance of 2 m ~ z (as for the binary codes previ- 
ously described with PMEPR of at most 3 dB), and any 
appropriate decoder described above can be used. 

Similar trade-offs between data transfer rate and PMEPR 
can be made with quaternary, octary and higher-order codes. 
Table 8 shows coset representatives for quaternary codes of 
length 16 symbols, ordered by maximum PEP over each 
coset, in the same format as Table 7: 

TABLE 8 



(3) 



Quaternary length 16 
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60 



'Additional row' coefficients 



PEP over 



It should be noted that the derivation of the PEP over a coset 
involves calculating the amplitude of the envelope of the 65 
composite waveform produced by a code word at successive 
sample points along the envelope in the time domain. The 



12 


13 


14 


23 


24 


34 


Coset representative 


coset 0* 


0 


0 


1 


1 


0 


1 


0002 0020 0200 0222 


3X00 


0 


1 


0 


0 


1 


1 


0002 0200 0020 0222 


3X00 


0 


1 


0 


1 


1 


0 


0000 0220 0022 0202 


3X00 


0 


1 


1 


1 


0 


0 


0000 0022 0220 0202 


3Z00 




0 


0 


0 


1 


1 


0002 0200 0002 2022 


3X00 




0 


0 


1 


0 


1 


0002 0020 0002 2202 


3X00 




0 


1 


0 


0 


1 


0002 0002 0200 2022 


3X00 




0 


1 


1 


0 


0 


0000 0022 0202 2002 


3X00 




1 


0 


0 


0 


1 


0002 0002 0020 2202 


3X00 




1 


0 


0 


1 


0 


0000 0202 0022 2002 


3X00 


0 


0 


1 


1 


1 


0 


0000 0220 0202 0022 


3X00 


0 


1 


1 


0 


1 


0 


0000 0202 0220 0022 


3X00 


0 


0 


0 


1 


0 


1 


0002 0020 00O2 0020 


64.00 


0 


0 


0 


1 


1 


0 


0000 0220 0000 0220 


64.00 


0 


0 


1 


0 


1 


1 


0002 0200 0200 0002 


64.00 


0 


0 


1 


1 


0 


0 


0000 0022 0202 0220 


64.00 


0 


0 


1 


1 


1 


1 


0002 0222 0200 0020 


64.00 



(1) 
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TABLE 8 -continued 
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Quaternary length 16 



PEP over 



12 


13 


14 


23 


24 


34 


Coset representative 


coset (W) 


0 


1 


0 


0 


0 


1 


Anno nnno n/v^n fw>n 

0002 0002 0020 0020 


64.00 


0 


1 


0 


0 


1 


0 


0000 0202 0022 0220 


64.00 


0 


1 


0 


1 


0 


1 


nrt/vi nrt^n **wv>a nnn*> 

0002 0020 0020 0002 


64.00 


0 


1 


0 


1 


1 


1 


0002 0222 0020 0200 


64.00 


0 


1 


1 


0 


0 


0 


0000 0000 0220 0220 


64.00 


0 


1 


1 


0 


1 


1 


0002 0200 0222 0020 


64.00 


0 


1 


1 


1 


0 


1 


nAftfl f\(\^t\ rtnoi (VinA 

0002 0020 0222 0200 


64.00 


0 


1 


1 


1 


1 


1 


AAAO rt^^l nflfl^ fW\*> 

0002 0222 0222 0002 


64.00 


1 


0 


0 


0 


0 


1 


0002 0002 0002 2220 


64.00 


1 


0 


0 


0 


1 


0 


nftrtn nino nnnn ono n 

0000 0202 0000 2020 


64.00 


1 


0 


0 


1 


0 


0 


r\r\f\f\ niv>o nnnn oonn 

0000 0022 DU0U 2200 


64.00 


1 


0 


0 


1 


1 


0 


0000 0220 0000 2002 


64.00 




0 


0 


1 


1 


1 


0002 0222 0002 2000 


64.00 


1 


0 


1 


0 


0 


0 


0000 0000 0202 2020 


64.00 


1 


0 


1 


0 


1 


1 


0002 0200 0200 2220 


64.00 




0 


1 


1 


1 


0 


nnnn /1T>/1 Aini ^Vlrtrt 

0000 0220 0202 2200 


64.00 


I 


0 


1 


1 


1 


1 


0002 0222 0200 2202 


64.00 






o 


Q 


Q 


0 


nnnn nnnn nrv>'? omn 
uuuu uuuu uuzz zzuu 


64 00 


1 


1 


0 


1 


0 


1 


nnnn nivin r\fV%t\ *VT1A 

0002 0020 0020 2220 


64.00 




1 


0 


1 


1 


0 


itnnn ntvyn nn^i '*wv»n 
0000 0220 0022 2020 


64.00 


1 


1 


1 


0 


0 


0 


nnnn nnnn (M/in ortrt'* 

0000 0000 0220 2002 


64.00 


1 


1 


1 


0 


0 


1 


nrtrti nr\Ai innrt 

0002 0002 0222 2000 


64.00 




1 


1 


0 


1 


0 


nnnn n^ni /vion wnn 

0000 0202 0220 2200 


64.00 




1 


1 


0 


1 


1 


0002 0200 0222 2202 


64.00 




1 


1 


1 


0 


0 


nnnn nnoi M<^n nnon 

0000 0022 0220 2020 


64.00 


0 


0 


0 


0 


1 


1 


0002 0200 0002 0200 


64.00 


0 


0 


1 


0 


0 


1 


0002 0002 0200 0200 


64.00 


0 


0 


1 


0 


1 


0 


0000 0202 0202 0000 


64.00 


0 


1 


0 


1 


0 


0 


0000 0022 0022 0000 


64.00 


0 


1 


1 


1 


1 


0 


0000 0220 0220 0000 


64.00 




0 


1 


1 


0 


1 


0002 0020 0200 2000 


64.00 




1 


0 


0 


1 


1 


0002 0200 0020 2000 


64.00 




1 


0 


1 


1 


1 


0002 0222 0020 2022 


64.00 




1 


1 


1 


0 


1 


0002 0020 0222 2022 


64.00 




1 


1 


1 


1 


0 


0000 0220 0220 2222 


64.00 
















(2) 


1 


1 


0 


1 


0 


0 


0000 0022 0022 2222 


128.00 


0 


0 


0 


0 


0 


1 


0002 0002 0002 0002 


128.00 


0 


0 


0 


0 


1 


0 


0000 0202 0000 0202 


128.00 


0 


0 


0 


1 


0 


0 


0000 0022 0000 0022 


128.00 


0 


0 


0 


1 


1 


1 


0002 0222 0002 0222 


128.00 


0 


1 


1 


0 


0 


1 


0002 0002 0222 0222 


128.00 


1 


0 


1 


0 


1 


0 


0000 0202 0202 2222 


128.00 


0 


0 


1 


0 


0 


0 


0000 0000 0202 0202 


128.00 


0 


1 


0 


0 


0 


0 


0000 0000 0022 0022 


128.00 


1 


0 


0 


0 


0 


0 


0000 0000 0000 2222 


128.00 


1 


1 


1 


1 


1 


1 


0002 0222 0222 2220 


128.00 


0 


0 


0 


0 


0 


0 


0000 0000 0000 0000 


256.00 



10 



15 



20 



25 



live generator array. Other coset representatives for the base 
codes could be included to generate larger tables in similar 
manner, but it is generally less desirable to select code words 
for use from such tables solely on the basis of tabulated PEP, 
because such selections of code words would be likely to 
have inferior error detection and correction properties. It is 
also noted that the generator arrays used in driving Tables 7 
and 8 are the same generator arrays used in deriving the 
respective minimum PMEPR binary and quaternary codes 
described earlier. 

Table 9 shows coset representatives for octary codes of 
length 16 symbols, ordered by maximum PEP over each 
coset. In this case the first column contains six coefficients 
constituting multipliers for 2 m -symbol words which are of 
the form 2(X ( >*X Jt ) (whereas the additional rows of the octary 
generator array are of the form 4(X,*X^, and which are 
added symbol-wise modulo 2/ to create the coset represen- 
tative; thus, for example, a 3 under the heading '12' indicates 
that the symbols of the word 2(K 1 *X 2 ) are each multiplied 
by 3 for inclusion in the coset representative. Thus in this 
case the coset representatives listed are not limited to linear 
combinations of additional rows in the octary generator 
array. In view of the length of this table, only selected rows 
are given here: 

TABLE 9 



30 



35 



40 



45 



This table can be used in a similar manner to Table 7. Thus, 
for example, the code words associated with the first twelve 
cosets in Table 8, above dashed line (1), together with those 
associated with any twenty of the next forty cosets, between 
dashed lines (1) and (2), can be used to encode up to fifteen 
data bits as quaternary code words sixteen symbols long; the 
maximum PMEPR of a multi-carrier signal modulated in 
accordance with those code words is four (6 dB). 

Tables similar to Table 8 can be prepared for other lengths 
of quaternary code word. For a code word length of 2 m , any 
code comprising a set of cosets from such a table has a 
minimum Hamming distance of 2 m ~ 2 and a minimum Lee 
distance of 2 m ~ x (as for the quaternary codes previously 
described producing a PMEPR of at most 3 dB); any 
appropriate decoder described above can be used with such 
a code, and retains the associated properties already dis- 
cussed. 

Each of the coset representatives listed in Tables 7 and 8 
is a linear combination of the additional rows in the respec- 



55 



60 



65 



Octarv length 16 



"Additional row' coefficients 



PEP over 



32 


13 


14 


23 


24 


34 


Coset representative 


coset (W) 


0 


0 


2 


2 


0 


2 


0004 0040 0400 0444 


32.00 


2 


0 


0 


0 


2 


2 


0004 0400 0004 4044 


32.00 


2 


0 


0 


2 


0 


2 


0004 0040 0004 4404 


32.00 


2 


0 


2 


0 


0 


2 


0004 0004 0400 4044 


3100 


2 


0 


2 


2 


0 


0 


0000 0044 0404 4004 


32.00 


2 


2 


0 


0 


0 


2 


0004 0004 0040 4404 


32.00 


2 


2 


0 


0 


2 


0 


0000 0404 0044 4004 


32.00 


0 


0 


2 


2 


2 


0 


0000 0440 0404 0044 


32.00 


0 


2 


0 


2 


2 


0 


0000 0440 0044 0404 


3Z00 


0 


2 


0 


0 


2 


2 


0004 0400 0040 0444 


3200 


0 


2 


2 


0 


2 


0 


0000 0404 0440 0044 


3X00 


0 


2 


2 


2 


0 


0 


0000 0044 0440 0404 


3X00 


1 


0 


2 


2 


1 


0 


0000 0246 0404 2064 


48.00 


1 


2 


1 


0 


2 


0 


0000 0404 0246 2064 


48.00 


0 


0 


2 


2 


1 


1 


0002 0240 0406 0644 


48.00 


0 


1 


2 


2 


0 


1 


0002 0046 0420 0464 


48.00 


1 


2 


0 


1 


2 


0 


0000 0426 0044 2604 


48.00 


0 


2 


0 


1 


2 


1 


0002 0420 0046 0464 


48.00 


0 


2 


1 


0 


2 


1 


0002 0406 0240 0644 


48.00 


2 


0 


0 


1 


1 


2 


0004 0220 0004 4664 


48.00 


2 


0 


1 


0 


1 


2 


0004 0206 0206 4044 


48.00 


2 


1 


0 


1 


0 


2 


0004 0026 0026 4404 


48.00 


1 


1 


2 


2 


0 


0 


0000 0044 0426 2604 


48.00 


2 


1 


1 


0 


0 


2 


0004 0004 0220 4664 


48.00 


0 


2 


3 


0 


2 


3 


0006 0402 0640 0244 


48.00 


2 


0 


0 


3 


3 


2 


0004 0660 0004 4224 


48.00 


2 


0 


3 


0 


3 


2 


0004 0602 0602 4044 


48.00 


2 


3 


0 


3 


0 


2 


0004 0062 0062 4404 


48.00 


2 


3 


3 


0 


0 


2 


0004 0004 0660 4224 


48.00 


3 


2 


0 


3 


2 


0 


0000 0462 0044 6204 


48.00 


3 


2 


3 


0 


2 


0 


0000 0404 0642 6024 


48.00 


3 


3 


2 


2 


0 


0 


0000 0044 0462 6204 


48.00 


0 


0 


2 


2 


1 


3 


0006 0244 0402 0640 


48.00 


0 


0 


2 


2 


3 


1 


0002 0644 0406 0240 


48.00 


0 


0 


2 


2 


3 


3 


0006 0640 0402 0244 


48.00 


0 


2 


0 


3 


2 


3 


0006 0460 0042 0424 


48.00 


0 


2 


1 


0 


2 


3 


0006 0402 0244 0640 


48.00 


0 


3 


2 


2 


0 


3 


0006 0042 0460 0424 


48.00 


1 


3 


2 


2 


0 


0 


0000 0044 0462 2640 


48.00 


3 


0 


2 


2 


3 


0 


0000 0642 0404 6024 


48.00 
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TABLE 9-continued 



TABLE 9-contimied 



Oclarv length 16 



PEP over 



12 


13 


14 


23 


24 


34 


Coset representative 


coset (W) 


0 


2 


0 


3 


2 


1 


nnn/i f\A£.A f\f\A £ fV/tOn 

0002 0404 UU46 0420 


48.00 


0 


1 


2 


2 


0 


3 


f\(\f\A nCMI A/fly! f\A Cf\ 

000 0 0042 U424 0460 


48.00 


2 


0 


0 


1 


3 


2 


rifvi/i nfi'JA nr\(\A A"i^r\ 
UUU^ *tZOU 


48.00 


2 


0 


1 


0 


3 


2 


firtrtj* fiAnf mft£ aaac\ 
0UU4 UOUi UZUO *W4U 


48.00 


2 


1 


0 


3 


0 


2 


f\r\t\A nn&i nni*c /i/i/tr» 
0U04 UUOi LHJZO 


48.00 


] 


0 


2 


2 


3 


0 


nnnn ruC^II fY/tfl/t 

0UU0 U04Z 


48.00 


1 


2 


3 


0 


2 


0 


nnon n/tfiA nfiAi i/iytn 
UUUU U Wt UO*ti Z40U 


/is m 


0 


2 


0 


1 


2 


3 


nnrtri (\aoa ruvn rvi/tn 
UUUO U4£4 lAWx WOU 


40, uu 


0 


2 


3 


0 


2 


1 


nnni n/inc n£/i>t m An 
UUU2 U4O0 UD44 UZ4U 


48.00 


0 


3 


2 


2 


0 


1 


0UU2 UU40 U404 U4ZU 


48.00 


2 


0 


0 


3 


1 


2 


nnn/i nn;>i nnnyi /i^in 
UUU*4 U/04 UUU4 


48.00 


2 


0 


3 


0 


1 


2 


nrtrivi rnnit n^no a a An 
UUU4 UZUO UOUZ 444U 


48.00 


2 


1 


3 


0 


0 


2 


UUU4 UUU4 UoZ4 4ZoU 


48.00 


2 


3 


1 


0 


0 


2 


n nrt A n nn a ni a a a c.*) t\ 
UUU4 UUU4 UZt>4 4oZU 


48.00 


1 


2 






2 


0 


0000 0462 0044 2640 


48 00 


3 


0 


2 


2 


1 


0 


0000 0246 0404 6420 


48.00 


3 


1 


2 


2 


0 


0 


0000 0044 0426 6240 


48.00 


2 


3 


0 


1 


0 


2 


0004 0026 0062 4440 


48.00 


3 


2 


0 


1 


2 


0 


0000 0426 0044 6240 


48.00 


3 


" 


1 


0 


2 


0 


0000 0404 0246 6420 


48.00 


0 


0 


2 


2 


0 


1 


0002 0046 0406 0442 


54.63 


0 


2 


0 


0 


2 


1 


0002 0406 0046 0442 


54 63 


0 


0 


2 


2 


2 


3 


0006 0446 0402 0042 


54.63 


0 


2 


0 


2 


2 


3 


0006 0446 0042 0402 




0 


2 


2 


0 


2 


3 


0006 0402 0446 0042 


54.63 


0 


2 


2 


2 


0 


3 


0006 0042 0446 0402 


54.63 


0 


0 


2 


2 


3 


2 


0004 0646 0400 0242 


54.63 


0 


0 


2 


2 


1 


0 


0000 0246 0404 0642 


54.63 


2 


0 


0 


2 


3 


2 


0004 0646 0004 4202 


54.63 


2 


0 


0 


0 


1 


2 


0004 0206 0004 4642 


54.63 


2 


0 


2 


0 


3 


2 


0004 0602 0400 4642 


54.63 


2 


0 


2 


2 


3 


0 


0000 0642 0404 4602 


54.63 










(132 rows omitted) 




3 


2 


2 


2 


0 


3 


0006 0042 0446 6260 


58.90 


1 


2 


2 


2 


0 


1 


0002 0046 0442 2620 


58.90 


2 


2 


3 


3 


0 


2 


0004 0062 0646 4260 


58.91 


2 


2 


1 


1 


0 


2 


0004 0026 0242 4620 


58.91 


3 


2 


2 


0 


2 


3 


0006 0402 0446 6620 


58.97 


1 


2 


2 


0 


2 


1 


0002 0406 0442 2260 


53.97 


2 


3 


2 


0 


3 


2 


0004 0602 0466 4620 


58.97 


2 


1 


2 


0 


1 


2 


0004 0206 0422 4260 


58.97 


2 


1 


2 


2 


3 


0 


0000 0642 0426 4624 


58.98 


2 


3 


2 


2 


1 


0 


0000 0246 0462 4264 


58.98 


2 


2 


3 


1 


2 


0 


0000 0426 0642 4624 


58.98 


2 


2 


1 


3 


2 


0 


0000 0462 0246 4264 


58.98 


2 


1 


2 


3 


0 


0 


0000 0066 0426 4040 


58.99 


2 


3 


2 


1 


0 


0 


0000 0022 0462 4040 


58.99 


2 


3 


0 


3 


2 


0 


0000 0462 0066 4004 


59.00 


2 


1 


0 


1 


2 


0 


0000 0426 0022 4004 


59.00 


2 


2 


1 


0 


3 


0 


0000 0606 0246 4400 


59.00 


2 


2 


3 


0 


1 


0 


0000 0202 0642 4400 


59.00 


2 


2 


1 


1 


2 


0 


0000 0426 0246 4220 


59.03 


2 


2 


3 


3 


2 


0 


0000 0462 0642 4660 


59.03 


2 


1 


2 


2 


1 


0 


0000 024 6 0426 4220 


59.03 


2 


3 


2 


2 


3 


0 


0000 0642 0462 4660 


59.03 


3 


0 


2 


0 


3 


2 


0004 0602 0400 6064 


59.05 


1 


0 


2 


0 


1 


2 


0004 0206 0400 2024 


59.05 


2 


0 


3 


2 


3 


0 


0000 0642 0606 4004 


59.05 


2 


0 


1 


2 


1 


0 


0000 0246 0202 4004 


59.05 


0 


1 


2 


0 


2 


1 


0002 0406 0420 0024 


59.06 


0 


3 


2 


0 


2 


3 


0006 0402 0460 0064 


59.06 


2 


1 


2 


1 


0 


0 


0000 0022 0426 4004 


59.09 


2 


3 


2 


3 


0 


0 


0000 0066 0462 4004 


59.09 


1 


2 


0 


1 


0 


2 


0004 0026 0040 2204 


59.12 


3 


2 


0 


3 


0 


2 


0004 0062 0040 6604 


59.12 










(3,850 


rows omitted) 




0 


0 


0 


0 


0 


1 


0002 0002 0002 0002 


218.51 


0 


0 


0 


1 


0 


0 


0000 0022 0000 0022 


218.51 


0 


0 


3 


0 


0 


0 


0000 0000 0606 0606 


218,51 


0 


3 


0 


0 


0 


0 


0000 0000 0066 0066 


218.51 


3 


0 


0 


0 


0 


0 


0000 0000 0000 6666 


218.51 


0 


0 


0 


0 


1 


0 


0000 0202 0000 0202 


218.51 


0 


0 


0 


0 


3 


0 


0000 0606 0000 0606 


218,51 



20 



25 



35 



45 



55 



60 
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Octary length 16 
'Addjtjopal row' coeflcjepts 



12 13 14 23 24 34 



PEP over 
Coset representative coset (W) 



0 0 0 3 0 0 0000 0066 0000 0066 218.51 
0 0 0 0 0 3 0006 0006 0006 0006 218.51 
10 0 0 0 0 0 0 0000 0000 0000 0000 256.00 



In producing Table 9 the calculation of PEP has been 
repeated for successively shorter intervals between envelope 
sample time values until at least the first decimal place, after 
rounding, is stable. 

It can be seen from Table 9 that, in the case of octary 
codes, there are a great many cosets available for which the 
PMEPR is not even as high as four. 

Tables similar to Table 9 can be prepared for other lengths 
of octary code word. For a code word length of 2 m , any code 
comprising a set of cosets from such a table has a minimum 
Hamming distance of 2 m ~ 2 and a minimum Lee distance of 
2 m ~ 1 (as for the octary codes previously described for which 
the PMEPR is at most 3 dB). This is so even though, in the 
octary case, the coset representatives listed are not limited to 
linear combinations of additional rows in the octary genera- 
tor array. This unexpected feature contributes significantly to 
30 the very large number of octary cosets associated with low 
PMEPR. 

Each of the coset representatives listed in Table 9 is a 
linear combination of 2 m -symbol words of the form 
2(X I *X A ). Even more octary cosets can be made available by 
allowing linear combinations of 2 m -symbol words of the 
form (X^Xjt); in this case the minimum Hamming distance 
of the codes remains unchanged although the minimum Lee 
distance may decrease. In practice this option would be 
40 exploited by adding rows to the selector 252 in the encoder 
230 (FIG. 14), thereby enabling the selector to accept 
additional input data bits for determining selection of a coset 
representative. The octary decoder described above can be 
used with any such code, and retains the associated prop- 
erties discussed above. 

Still more coset representatives for the base code could be 
included to generate larger tables in similar manner, but it is 
generally less desirable to select code words for use from 
50 such tables solely on the basis of tabulated PEP, because 
such selections of code words would be likely to have 
inferior error detection and correction properties. 

Tables similar to Table 9 can be produced for higher-order 
codes, such that each of the coset representatives is a linear 
combination of 2 m -symbol words of the form 2(X ( *X Jt ). For 
a code word length of 2 ra , any code comprising a set of 
cosets from such a table has a minimum Hamming distance 
of 2 m " 2 and a minimum Lee distance of 2 m_1 . Even more 
cosets can be made available by allowing linear combina- 
tions of 2 m -symbol words of the form (X^X*); in this case 
the minimum Hamming distance of the codes remains 
unchanged although the minimum Lee distance may 
decrease. The appropriate decoder described above can be 
used with any such code. 

Ordered lists of cosets, such as in Tables 7, 8 and 9, can 
be produced using the following program in the C program- 
ming language: 
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/• Orders the coscts of W(M) in Y(M) according to the maximum peak power 
over the coset. W(M) is the 2 J-phase linear code having the same 
formal generator matrix as RM(1,M). For J=l (binary), Y(M) is RM(2,M). 
For J>1 (quaternary, octary, etc), Y(M) is the 2" J-phase linear code 
having the same formal generator matrix as RM(2,M), except that the 
entries of the rows corresponding to quadratic terms are multiplied by 
2. 

V 

^define M 4 /* Integer M >= 2 */ 

#define J 1 /* Integer J >- 1. Program assumes MJ < 32. 

Program also assumes M(M-l)/2 < 32 for J - 1, 
and (J-l)M(M-l)/2 < 32 for J > 1 7 
^define OVER 356 /* Ovcrsampling multiple 7 
tfdeftne NON_BINARY (M)/* Flag for all cases except binary (J=l) 7 
#define L (1 <:< M) /* Sequence length L ■ 2 M 7 

^define P (1 « J) /* Allowed sequence symbols are 0,1,2^-1 7 

#if NON_BINARY 

^define REPS ((unsigned long)l « ((J-l)*M*(M-l)/2))/* # of coset reps 7 
^define QUAD^MULT 2 /* multiple for quadratic terms 7 

#else 

#define REPS ((unsigned long)l « M*(M-l)/2) 

^define QUAD.34ULT 1 

#endif 

#define PI 3.14159 

^include <studio.h> 
#include <math.h> 
double power( ); 

double max„power( ); 

double max_power_coset( ); 

int gen[M+llL]; /* Generator matrix for W(M) 7 

typedef struct 

{ 

int a[M+l!M+l]; /* Coefficients of quadratic terms7 

int seojLj /* Sequence comprising quadratic terms 7 

double max_power, /* Maximum power over coset */ 
} coset_jep; 

coset__rep list(REPS]; /* List of coset reps, ordered by max power */ 

main( ) 

< 

coset_rep s; 

unsigned long mask; 
unsigned long c; 
int i, j, k; 

/* Initialise gen[ J ] 7 
initialise( ); 

for (c = 0; c < REPS; C++) 
{ 

/* Set up coefficients of quadratic terms 7 
mask - 1; 

for (i = M; i >= 1; i — ) 
for (] = M; j >= i+1; j — ) 
{ 

s.a[i][j] - (c & mask 7 1:0); 

mask «» 1; 

for (k - 1; k < J-l; k++) 

{ 

s.afilj] - ^ * SA[i][j] + (c & mask ? 1 : 0); 
mask «= 1; 

} 

/* Each s.a[i][j] takes values 0 or 1 for binary, 
and in range 0,l,2,...,P/2 -1 otherwise */ 

} 

/* Set up sequence as 

6um_{j>i} QUAD_MULTs.a[iIj] gen(il ] gen[fl ] 7 
for (k - 0; k < L; k++) 

s.seq[k] = 0; 
for (i - M; i >- 1; i — ) 
for (j - M; j s— i+1; j — ) 
{ 

for (k = 0; k < L; k++) 

s.seqlk] - (s.seq[k] + QUAD_MULT * 6.a[i][j] * 
gen[ilk] - genljflkD % P; 

} 
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/* Calculate maximum power over the coset of W(M) whose 
representative is s.scq^ ]. Coset members have form 
s.sesq[] + sum_{i=0} {M} a[i] gen(il ] for 0 <- a[i < P. 

7 

s.max power = max_power_coset(s.seq); 

/* Update records */ 
update(&fi, c); 

} 

/• Print ordered list of coset reps 7 
print_Jist( ); 

} 

initialise( ) 
{ 

int af^l i> j; 

/" Set gen[ij[ ] to be i-th row of generator matrix for W(M) 

(where row 0 corresponds to X' and row i corresponds to X__i) */ 
for (j = 0; j < L; j ++) 
{ 

gen[01jj-l; 

extract_jsymbols(), 1, M, a); 
for (i - 1; i <- M; 
gen[ilj'] = a[i-l]; 

} 

} 

/" Store the s least sig (2 u)-ary symbols of the integer c, in order 
from most sig to least sig, in seq[0], seq[l], seq[s-l] 

7 

int extract_symbols(c, u, s, scq) 

int c, u, s, *seq; 

{ 

unsigned int mask; 
int i; 

mask - ( « u) -1; /* mask for u least sig bit positions 7 
for (i - s-1; i >- 0; i — , c »- u) 
seoji] =» c & mask; 

} 

double max_power_coset(seq) 

int seoJL]; 

{ 

unsigned long c, mask; 
double y, ooset__max; 

int a[M+ll d(Ll i, k; 
coset_max = 0.0; 

for (c m 0; c < ((unsigned long)l « M*J); C++) 
{ 

/* Set up coefficients e[ ] of linear terms */ 
mask = P-3; /* mask for J least sig bits 7 
for (i = M; i >~ 1; i— ) 
{ 

a[i] - (c & mask) » J*(M-i); 

/* Each a[i] takes values in 0,1,2,.. ,,P-1 */ 
mask <<- J; 

} ... 

/* Set up sequence as d[ ] = seq[ ] + sum__{i=l} {M} a[i] genfil ] 
No need to include row gen[0l ] in power calculations 7 

for (k - 0; k < L; k++) 

d[k] - seqtkj 
for (i - M; i >- 1; i — ) 
{ 

for (k = 0; k < L; k++) 

d[k] = (d[k] + a[i] * gen[ilk]) % P; 

} 

f* Do max power calculations 7 
y " max_poweT(d); 
if (y > coset_max) 
coset_max - y\ 

} 

retum(coset_max); 

} 

double max_power(d) 

int d[Lj 

{ 

double x, y, seq_max; 
seq_max = 0.0; 

for (x«0.0; x<1.0; x+-1.0/(OVER*L)) /* Sample (OVER)L points */ 
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{ 

y » power(d, x); 
if (y > seq_max) 
seq_jnax - y; 

} 

retuin(seq_max); 

} 

double power(d, x) 
int dlLj 
double x; 
{ 

double cO, cl, arg; 
int n; 

cO - cl = 0.0; 

for (n=l; n<=L; n++) 

{ 

arg = 2*PI* (n*x + (double)d[n-l]/P); 
cO cos(arg); 
cl +-= sLn(arg); 

} 

retuxn(cO*cO +cl*cl); 

} 

update(s_ptr, c) 
coset_rep *s_ptr; 
unsigned long c; 

int u, m; 

for (u - 0; u < c; u++) 

{ 

if (s ptr -> max_power < list[u].max power) 

{ 

/* Shove up all subsequent entries of list[ ] */ 
for (m » c; m >- u+1; m — ) 

listjm] = Ust[m-li 
/* Replace listfuj 7 
list[u] o * s _ptr; 
return; 

} 

} 

/* Place *s_pu at end of list( ] V 
Iist(u] «• *s_ptr, 

} 

print listf ) 

{ 

int i, j, k, u; 
for (i - 1; i <= M; 
for (j - i+1; i <° M; 

printf( M %2d%d" i, J); 
printff" %d times oversampling\n\", OVER); 
for (u = 0; u < REPS; u++) 
{ 

for (i = 1; i <= M; i++) 
for (j o i+1; j <b M; 

{ 

printf("%3d M , list[u].a(ilj]); 

} 

printf( w "); 

for (k - 0; k < L; k++) 

printf("%2d", list[u]j>eq[k]); 
printf("%10.2f\n", list [u>max_power); 

} 

} 



As well as choosing trade-offs between data transfer rate 
and PMEPR, it is possible to trade off, for example, mini- 
mum Hamming or Lee distance against data transfer rate. 
For example, by choosing a subset of four of the cosets 
defined by the twelve top rows in Table 7, the minimum 
distance can be increased from 4 to 6 in return for a 
reduction in number of data bits from 8 to 7, while preserv- 
ing a maximum PMEPR of 3 dB. An example of such a 
subset comprises the cosets defined by the second, third, 
fourth and sixth rows of Table 7. Similar trade-offs are 
possible for other lengths of code word and higher-order 
codes. In the extreme case a single coset of the base code 



may be used, providing a minimum Hamming and Lee 
distance of 2 m ~ l for any order code, but encoding only 
j(m+l) data bits. 

Identifying possible trade-offs between data transfer rate 

60 and PMEPR by constructing a list of coset representatives 
ordered by maximum PEP over each coset, as described 
above, has the disadvantage, especially for longer code 
words and/or higher-order codes, of requiring extensive 
calculation. The following technique identifies a set of coset 

65 representatives which can be combined with the base code 
to produce code words associated with a PMEPR no greater 
than a known value, without calculating the maximum PEP 
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over each coset (although this technique does not provide 
any ordering by PMEPR for the cosets in the set). The 
technique will be illustrated first in the case of a binary code 
of length 16 (m-4), that is with a generator array 
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10 



As a first step towards producing the set of coset 20 
representatives, an initial set of eight (2 n, ~ 1 ) expressions is 
constructed, each having the form 



where v l9 v 2 and V 3 are coefficients having values of 0 or 1 
and ® indicates bitwise addition modulo 2. Thus the initial 
set is 
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TABLE 10-continued 
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A code whose code words are obtained by combining any 
25 of these coset representatives with the base code has the 
following properties: 

the code's minimum Hamming distance is 4 (2 m ~ 2 ); 
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a multi-carrier signal modulated in accordance with these 
code words has a PMEPR of no more than 6 dB; 



This initial set is then extended by applying all possible 
permutations of the subscripts of the row identifiers X J( 
X 3 and X 4 in each expression. A final set is produced by 

identifying the distinct expressions within the extended set; 4S 

for this purpose, terms with the same pair of subscripts but a decoder based on Hamming distance for use with such 

in reverse order (such as Xj*X2 and X 2 *Xj) are regarded as a code can operate using analytical techniques and can 

being the same, and the ordering of such terms within an be implemented, for example, using combinatorial 

expression is disregarded. logic at least in part; 

TTiis final set of expressions, forty-nine in all, denning the so ^ m suffident code available fe enab , e m 

required set of coset representatives, is given in Table 10, in accepta ble data transfer rate to be attained. 

terms of six coefficients which constitute multipliers for the 

respective additional rows of the code's generator array In the S eneral case of a bma y of len S m 2 " blts 

which are added modulo 2 to create the corresponding coset whcre m = 3 > thc sct of 2 expressions is constructed 

representative; thus, for example, a 1 under the heading '12* 55 with me f° rm 
indicates that the additional row (X^XJ is multiplied by 1 

for inclusion in the coset representative, whereas a 0 inch- QC l m X^^C 2 *^® • * • ©(^0*^0^1(^1 '^M^aC^^m)® 

cates it is multiplied by 0 (i.e. effectively omitted). * ' ' ^-iP^-i*^) 



TABLE 10 



60 



12 13 14 23 24 34 12 1 3 14 23 24 34 



where v u v^ . . . v m-1 are coefficients having values of 0 or 
1, and the permutation of subscripts and identification of 
1 00100 0 10101 distinct expressions performed as described above. The 

q 10100 1 11000 65 number of cosets and the total number of code words 
0 01010 1 01110 identified by this procedure for various code word lengths 

(in terms of the parameter m) are given below: 
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m 


Cosets 


No. of code words 


3 


7 


112 


4 


49 


1568 


5 


552 


35328 


6 


7800 


998400 


7 


126360 


32348160 



A modification of this technique, applicable for word 
lengths of 64 or more bits (ml 6), enables a large subset of 
these cosets to be identified, without the need for examining 
the extended set of expressions to find and discard expres- 
sions which are not distinct. In other words, this modified 
technique can be guaranteed to produce an extended set of 
expressions which are inherently distinct. To this end, the 
combinations of values of the coefficients v lt v 2 , . . . v m _j are 
limited to combinations such that the binary words of the 
form v a v 2 v 3 . , . v m _ 2 v m-1 have the following two properties: 

each word v a v 2 v 3 , . . v m _ 2 v m _ 1 has Hamming weight of 
at least 4; and 

LCn-3V2J l(m-3V2J 

Constructing an initial set of expressions with the form 
given above but with this limited set of combinations of 
coefficients, and applying all possible permutations of the 
subscripts of the row identifiers in each expression, directly 
yields a final set of distinct expressions defining the required 
subset of coset representatives. The code words produced by 
combining the base code with any coset representative in 
this subset result in a PMEPR of at most 6 dB. At least m!.b 
2 m ~ 4 coset representatives are identified by this procedure 
when ml 7. 

Similar techniques can be used for identifying cosets of 
quaternary and higher-order code words which result in a 
known maximum PMEPR. For a quaternary code of length 
2 m where ml 3, an initial set of 2 m_1 expressions is con- 
structed of the form 

2(x l *x 2 )®2<x 2 *x 2 )® . . . ®iQt^ 2 x m _ 
i)©2v 1 pr 1 -Aj©2v 2 (y 2 *j£'j© . . . wv^&^rxj 

where v l7 v 2 , . . . v m-1 are coefficients having values of 0 or 
1 and © indicates symbol -wise addition modulo 4. This 
initial set is then extended by applying alt possible permu- 
tations of the subscripts of the row identifiers X lt X 3 , . . . 
X m in each expression. A final set of expressions, defining 
the required set of coset representatives, is produced by 
identifying the distinct expressions within the extended set, 
as previously described for the binary case. 

A quaternary code whose code words are obtained by 
combining any of these coset representatives with the qua- 
ternary base code has the following properties: 

the code's minimum Hamming distance is 2 m ~ 2 and its 

minimum Lee distance is 2 m_1 ; 
a multi-carrier signal modulated in accordance with these 

code words has a PMEPR of no more than 6 dB; 
a decoder for use with such a code can operate using 
analytical techniques and can be implemented, for 
example, using combinatorial logic at least in part; 
there are sufficient code words available to enable an 
acceptable data transfer rate to be attained; in general, 
this procedure identifies at least m!.^"* distinct coset 
representatives when ml 7. 
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The number of cosets and the total number of quaternary 
code words identified by this procedure for various code 
word lengths (in terms of the parameter m) are given below: 



5 



10 



m 


Cosets 


No. of code words 


3 


7 


1792 


4 


49 


50176 


5 


552 


2260992 


6 


7800 


127795200 


7 


126360 


8281128960 



As in the binary case, for word lengths of 64 or more 
15 symbols (ml 6) a large subset of these cosets can be 
identified, without the need for examining the extended set 
of expressions to find and discard expressions which are not 
distinct. This is accomplished by limiting the combinations 
of values of the coefficients v lt v 2 , . . . v m _ a to combinations 
20 such that the binary words of the form v^Vj . . . v TO , 2 v m _ 1 
have the two properties specified in the discussion above for 
the binary case. For ml 7 at least m!.2 m " 4 distinct coset 
representatives are identified by this procedure, and the code 
words produced by combining the quaternary base code with 
25 any coset representative in this subset yield a PMEPR of at 
most 6 dB. 

For quaternary code words of length 2 m where m!3, a 
modification is possible, by allowing the coefficients v a , 
v 2> • * * v m-i to nave anv °f me f° ur values 0, 1, 2 or 3 and 
30 modifying the form of the expressions in the initial set to be 

i(x x *x£®i{xsxi)® . . . ®nx„_ 2 *x m _ 

(i.e. without a multiplier of 2 before each coefficient). This 
35 initial set is then extended by applying all possible permu- 
tations of the subscripts of the row identifiers X lt X^ Xg, . 
. . X m in each expression. Afinal set of expressions, defining 
the required set of coset representatives, is produced by 
identifying the distinct expressions within the extended set, 
40 as previously described for the binary case. The number of 
cosets and the total number of quaternary code words 
obtained with this modification for various code word 
lengths (in terms of the parameter m) are given below: 

45 



m 


Cosets 


No. of code words 


3 


37 


9472 


4 


661 


676864 


5 


14472 


59277312 


6 


360600 


5908070400 


7 


£1290240 


^84557168640 



It can be seen that many more cosets are available in this 
55 case, so there is a potential for a higher data transfer rate. 
The code words produced by combining the quaternary base 
code with any coset representative identified by this modi- 
fied procedure yield a PMEPR of at most 6 dB; however, the 
minimum Lee distance with this modification is reduced to 
60 2 m ~ 2 . 

An additional modification produces an extended set of 
expressions which are guaranteed to be distinct after per- 
muting the subscripts of the row identifiers, by fixing the 
values of the coefficients v 1 and v m _ 3 to be 1 and 3 
65 respectively while allowing the other coefficients to have 
any value from 0 to 3. For ml 4, at least m!.4 m ~ 3 distinct 
coset representatives are obtained with this approach. 
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For an octary code of length 2 m where m^3, an initial set 
of 4 m ~ 5 expressions is constructed of the form 



4(xsx^®4(x 2 *x 3 )® . . . eyx^ 2 *x m _ 

1 )©2v 1 (X l -Xj©2v 2 (X 2 *X m )© . . . ©2v„_ 1 (x, nH 



where v l7 v 2 , . . . v ro-1 are coefficients having values of 0, 1, 
2 or 3 and © indicates symbol-wise addition modulo 8. This 
initial set is then extended by applying all possible permu- 
tations of the subscripts of the row identifiers X 1( X^ t X 3 , . . . 

in each expression. A final set of expressions, defining 10 
the required set of coset representatives, is produced by 
identifying the distinct expressions within the extended set, 
as previously described for the binary case. 

Octary code words obtained by combining any of these 
coset representatives with the base code have the following 15 
properties: 

the code's minimum Hamming distance is 2 m ~ 2 and its 

minimum Lee distance is 2 m ~ 1 ; 
a multi-carrier signal modulated in accordance with these 

code words has a PMEPR of no more than 6 dB; 
a decoder for use with such a code can operate using 

analytical techniques and can be implemented, for 

example, using combinatorial logic at least in part; 
there are sufficient code words available to enable an 

acceptable data transfer rate to be attained; in general, 

this procedure identifies at least m!.4"" 3 distinct coset 

representatives when m>4. 
The number of cosets and the total number of octary code 
words for various code word lengths (in terms of the 30 
parameter m) are given below: 



m 


Co sets 


No. of code words 


3 


37 


151552 


4 


661 


21659648 


S 


14472 


3793747968 


6 


360600 


756233011200 


7 


* 1290240 


>2.16xl0 13 



20 



25 



70 



m 


Cosets 


No. of code words 


3 


169 


692224 


4 


5917 


193888256 


5 


243912 


63940067328 


6 


>5.52xl0* 


>1.15xl0" 


7 


>3.09xl0 8 


>5.19xl0 15 



It can be seen that many more cosets are available in this 
case, so there is a potential for a higher data transfer rate. 
The code words produced by combining the octary base 
code with any coset representative identified by this modi- 
fied procedure yield a PMEPR of at most 6 dB; however, the 
minimum Lee distance with this modification is reduced to 

2 m-2 

An additional modification produces an extended set of 
expressions which are guaranteed to be distinct after per- 
muting the subscripts of the row identifiers, by fixing the 
values of the coefficients v x and v m-1 to be any one of the 
following ordered pairs of values, where the value of v a is 
given first: 1,2; 1,3; 1,5; 1,6; 1,7; 2,3; 2,5; 2,6; 2,7; 3,5; 3,6; 
3,7; 5,6; 5,7; 6,7, The other coefficients may have any value 
from 0 to 7. For m^4, at least 15.m!.8 m_3 distinct coset 
representatives are obtained with this approach. 

These techniques for identifying cosets defining code 
words producing a PMEPR of no more than 6 dB can be 
extended to higher-order codes. More generally, large sets of 
coset representatives can be identified for codes containing 
2 m symbols (where m^3) each with 2f possible values, for 
which the PMEPR is guaranteed to be at most 2 P+1 , for any 
fixed integer p where 1 =p=m~2. To this end, p sets of 
coefficients are defined as follows 



35 



40 



As in the quaternary case, for word lengths of 16 or more 
symbols (m^4) a large subset of these cosets can be 
identified, without the need for examining the extended set 
of expressions to find and discard expressions which are not 45 
distinct. This is accomplished by fixing the values of the 
coefficients v x and v m _ a to be 1 and 3 respectively while 
allowing the other coefficients to have any value from 0 to 
3. Exactly m!.4 m " 3 distinct coset representatives are 
obtained with this approach. 50 

For octary code words of length 2 m where m^3, another 
modification is possible, by allowing the coefficients v 3 , 
v 2> • • • v m-i *° Dave anv °f me eight values 0 to 7 and 
modifying the form of the expressions in the initial set to be 

4(A r 1 *^ 2 )©4(Y,'^r 3 )© . . . ®Mx^_ 2 *x m _ 

(Le. without a multiplier of 2 before each coefficient). This 
initial set is then extended by applying all possible permu- 
tations of the subscripts of the row identifiers X lf X 2 , X 3 , . . . 60 
X m in each expression. A final set of expressions, defining 
the required set of coset representatives, is produced by 
identifying the distinct expressions within the extended set, 
as previously described for the binary case. The number of 
cosets and the total number of octary code words obtained 65 
with this modification for various code word lengths (in 
terms of the parameter m) are given below: 



J 



55 



where each coefficient can have an integer value from 0 to 
2/-1. A further p(p-l)/2 coefficients w jy , where l^i<j§p; 
are defined where each vt-j can have an integer value from 
0 to 2 y -l. For an initial choice of the coefficients v, y and w t> - 
an expression is constructed of the form 

2'- 1 cr/x 2 )©2>- 1 (* 2 -x ) )© . . . mt-'ix^sx^ 

evx^x^ev^ix^x^® . . . ®v hm . p (x m ./x m _^ l ) 

Qv^txsx^^ev^txsx^je . . . <&> Xrn - p (x m - P *x„_^ 

®»,*Px m Xj&>&Qr 2 m XM . - . ®v p ^<X m -p*X m ) 

® W l,j(X m -p+i *X m _ p+3 )® ^2,30^ - p< -2 *X m _p+^ 

® W \.*SXm-p+\ 'X m - p ^)(3w 2 ^ 4 fi n _ pf2 m X m _ ( ^^)(Bv>'3,4(X m -p- t .3*X m _ p+ 
*) 

where ® indicates symbol-wise addition modulo 2 ; ; further 
expressions of this form are constructed for the remaining 
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possible combinations of choices of the coefficients and 
w i/f to obtain an initial set of 2 /><m **' )+;>< ^ 1 > /2 expressions. 
This initial set is then extended by applying all possible 
permutations of the subscripts of the row identifiers X a , X^, 
X 3 , . . . X m in each expression. A final set of expressions, 5 
defining the required set of coset representatives, is pro- 
duced by identifying the distinct expressions within the 
extended set, as previously described for the binary case. 
The set of code words produced by combining the base code 
of order 2 J with any coset representative identified by this 10 
procedure have a minimum Lee distance and a minimum 
Hamming distance both of 2 m ~ 2 . In the case where j^2, a 
further modification is possible in which the coefficients v, y 
and w,y are restricted to have even values in the ranges 
specified above. In this case fewer cosets will be identified, 15 
but the resulting code will have a minimum Lee distance of 

2 m-l 

Additional trade -o fife between data transfer rate, PMEPR 
and error detection and correction capabilities can be 
obtained in other ways than those described above. For 20 
example, instead of using a 32-bit code word to encode 11 
data bits directly as described above, successive blocks of 8 
data bits could be encoded as 16-bit code words using a 
8B/16B code derived from Table 3, and then pairs of these 
code words concatenated to produce a composite 32-bit 25 
word for transmission. In comparison with direct use of the 
11B/32B code, such concatenation would have the following 
properties: 

a better data transfer rate (16 data bits represented by 32 
transmitted bits); and 30 

a PMEPR of, at worst, 4(=6 dB), i.e. double the maximum 
PMEPR for the codes used for the component code 
words. 

Thus concatenation in this manner could be useful in cir- 
cumstances where a good data transfer rate is desired, 35 
together with controlled, though not necessarily optimum, 
PMEPR. 

In decoding a concatenated code, the code word would be 
split into two successive component 16-bit code words, 
which would then be individually processed by a 16B/8B 40 
decoder 120 or 140 to retrieve respective 8-bit data words. 

Shorter-length code words can be combined to form 
longer composite code words for transmission in other ways 
than by concatenation of pairs of code words. Thus, for 
example, four component code words can be concatenated, 45 
and the received composite code word split into four por- 
tions for decoding. Alternatively, it is envisaged that the 
component bits or symbols of individual shorter code words 
can be interleaved, one by one or in pairs, quadruples, etc. 
The component code words of a composite code word need 50 
not be members of the same coset. Furthermore the choice 
of coset representative for a particular component code word 
can range over multiple values. In this case the identity of 
the coset to which each selected component code word 
belongs can itself be used to transfer information, e.g. by 55 
defining one or more additional bits of the data word, in a 
manner analogous to the function of, for example, the 
selector 70 in the encoder 60 of FIG. 8. 

The composite transmitted code word created by concat- 
enation or interleaving in this manner is not directly derived 60 
as described above from a generator matrix having rows of 
the same length, but from a generator matrix having shorter 
(e.g. half-length) rows. Nonetheless, it can be shown that the 
code constituted by such composite code words in fact 
comprises, and therefore can be selected from, a set of 65 
complete cosets of a code having a generator matrix with 
rows of the same length as the composite code words. 
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In another modification a 16-bit code word and a 32-bit 
code word, each selected from respective codes derived as 
described above, could be concatenated to produce a com- 
posite 48-bit word for transmission. This 48-bit word could 
be combined with 16 additional zero bits to produce a 64-bit 
word for processing with the IFFT for transmission; the 
additional 16 zero bits enable oversampling to be used upon 
reception, while providing the optimum choice of a power of 
two for the number of samples. Another possibility is to 
concatenate three 16-bit code words, each selected from a 
16-bit code, to produce the composite 48-bit code word. 

It is noted that the binary, quaternary and octary codes of 
a given length presented above are all derived from a single, 
common generator matrix, which can indeed also be used to 
derive similar higher-order codes. It is envisaged that this 
property would be valuable in the implementation of adap- 
tive systems for changing the order of code used according 
to prevailing characteristics of a communications path. 
Thus, when a high-quality path is in use, a higher-order 
code, such as octary, could be used to obtain a high rate of 
data transfer. If the quality of the path deteriorates, leading 
to increased error rates, it could be advantageous to switch 
to use of a lower-order code, such as quaternary or even 
binary. With the present invention this can be readily accom- 
plished without the need for complete replication of circuitry 
for each different-order code available. The same basic 
circuitry could be used, incorporating or based on the 
common generator matrix, together with switchable param- 
eters corresponding for example to the values in the addi- 
tional rows of each code's generator array. Selection of the 
appropriate values of these parameters would cause the same 
basic circuitry to generate the particular order of code 
required. 

Coset representatives derived as described above provide 
particularly advantageous properties. In the most general 
case a coset representative can be any value combined with 
code words in a base code by symbol- wise addition modulo 
2f t to shift their values in a systematic manner. 

The code words in the codes defined above also have 
useful properties for establishing synchronization of timing 
of circuit operation with a received signal. Accordingly, the 
circuitry for decoding these codes as described above can 
also be used to identify synchronization sequences in the 
received data, thereby providing economies in circuit imple- 
mentation. It is also envisaged that these codes have advan- 
tageous properties for communication by techniques other 
than COFDM; thus, for example, they may be used as codes 
for code division multiple access (CDMA) techniques or for 
modem transmission of digital data over analogue commu- 
nication paths. 

What is claimed is: 

1. A method of encoding data words for transmission 
using multi-carrier frequency division multiplexing com- 
prising: 

selecting a set of one or more cosets of a linear sub -code 
of a code having a generator matrix as follows: 
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((20 + 1) x T- y )}2 m ~ i ((Z 1 + 1) x 2 m " 1 )/ 2 mH 



74 



uzo + Dxaw 

((ZO + Dx^)/^ 
((2.0 + l)x2°)/2° 



((2.H-l)x2 J ')/2 jr 

((lH-Dx^)/^ 
((2.1 + l)x2 0 )/2 0 



((2.*4-l)x2 ff, )/2 iB 
((2.z^l)x2 m - l )/2 m - 1 

((2.z+-l)x2>)/2> 

((2.z+l)x2 1 )/2 1 
((2.z+l)x2°)/2° 



((l(2 B, -l) + i)x2 m )/2 m 
((2.(2 m -l) + I)x2 B, - l )/2 m - 1 

((2.(2" - I) + I) x 20/2' 

((2.(2" -1)+ i)x2 v )l2 l 
((2.(2* , -i)+l)x2 0 )/2° 



where z has integer values from 0 to 2 m -l, y has integer 
values from m to 0, x indicates bitwise multiplication 
of factors expressed as (m+l)-bit binary numbers, and 25 
the result of division is expressed as a single digit; and 

encoding the data words as 2™ -symbol code words in 
accordance with their value and code words selected in 
complete cosets from the set of one or more cosets to 
achieve a desired peak-to -mean envelope power ratio 30 
characteristic of the transmission. 

2. The method of claim 1, wherein m is at least 3. 

3. The method of claim 2, wherein m is at least 4. 

4. The method of claim 1, wherein m is four and said 
generator matrix is: 35 

(0000 0000 1111 1111) 
(0000 1111 0000 1111) 

(0011 0011 0011 0011) 40 
(0101 0101 0101 0101) 

(1111 1111 1111 1111). 

5. The method of claim 1, wherein m is five and said 45 
generator matrix is: 

(0000 0000 0000 0000 1111 1111 1111 1111) 
(0000 0000 1111 1111 0000 0000 1111 1111) 

(oooo 1111 0000 1111 oooo 1111 0000 1111) 

(0011 0011 0011 0011 0011 0011 0011 0011) 
(0101 0101 0101 0101 0101 0101 0101 0101) 

(1111 1111 mi 1111 1111 1111 1111 1111). 



are derived by symbol-wise multiplication of all possible 
pairs of the first m rows of the generator matrix and are 
further multiplied by a positive integer s, where s is a power 
of two less than 
5 12. The method of claim 11, wherein s has the value 1. 

13. The method of claim 11, wherein j is at least 2 and s 
has the value 2. 

14. The method of claim 11, wherein j is at least 3 and s 
has the value Tf" 1 . 

10 15. The method of claim 11, wherein m is four and said 
generator array is: 



(0000 0000 1111 1111) 
(0000 1111 0000 1111) 
(0011 0011 0011 0011) 
(0101 0101 0101 0101) 

(1111 1111 1111 1111) 

(0000 0000 0000 ssss) 
(0000 0000 OOss OOss) 
(0000 0000 OsOs OsOs) 
(0000 OOss 0000 OOss) 
(0000 OsOs 0000 OsOs) 
(000s OOOs OOOs 000s). 



20 



50 



55 



16. The method of claim U, wherein m is five and said 
generator array is: 

(oooo oooo oooo oooo 1111 1111 1111 1111) 

(oooo oooo 1111 1111 oooo oooo 1111 1111) 

(0000 1111 OOOO 1111 OOOO 1111 0000 1111) 
(0011 0011 0011 0011 0011 0011 0011 0011) 
(0101 0101 0101 0101 0101 0101 0101 0101) 

(1111 1111 1111 1111 1111 1111 1111 1111) 

(0000 0000 0000 0000 0000 0000 ssss ssss) 
(0000 0000 0000 0000 0000 ssss 0000 ssss) 
(0000 0000 0000 0000 OOss OOss OOss OOss) 
(0000 0000 0000 0000 OsOs OsOs OsOs OsOs) 
(0000 0000 0000 ssss 0000 0000 0000 ssss) 
(0000 OOOO OOss OOss 0000 OOOO OOss OOss) 
(0000 0000 OsOs OsOs 0000 0000 OsOs OsOs) 
(0000 OOss 0000 OOss 0000 OOss OOOO OOss) 
(0000 OsOs OOOO OsOs OOOO OsOs 0000 OsOs) 
(OOOs OOOs OOOs OOOs OOOs OOOs OOOs OOOs). 



6. The method of claim 1, wherein the number of code 
words is an integer power of two. 

7. The method of claim 1, wherein each symbol has 2 J 60 
different possible values, where j is a positive integer. 17. The method of claim 11, wherein all coset represen- 

8. The method of claim 7, wherein j has the value 1. tatives of cosets in said set of cosets have the form 2 7 " 1 

9. The method of claim 7, wherein j has the value 2. (X 1 *X 2 )©2 / * 1 (X 2 *X 3 )© . . . ©2 >_1 (X m _ 1 *X m ), or a permu- 

10. The method of claim 7, wherein j has the value 3. tation of those subscripts, where X 2 to X m are row labels of 

11 . The method of claim 7, wherein all said code words 65 the first m rows in the generator matrix, * indicates symbol - 
also belong to a second code having a generator array wise multiplication and © indicates symbol -wise addition 
comprising the generator matrix plus additional rows which modulo 7f. 
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((2.0 + l)x2" , )/2" ((Xl + l)x2" , )/2 ,n 
((2.0-t-l)x2 m_1 )/2 m ~ l ((2.1 •Kl)x2 n, " I )/2 m " 1 



10 



15 



18. The method of claim 1, wherein said code words are 
selected from a list of cosets of a linear sub -code of said 
code, said list being ordered to indicate maximum peak 
envelope power over any code word in a coset relative to 
that for other cosets. 

19. The method of claim 18, wherein at least a predeter- 
mined number of cosets is selected such that the maximum 
peak-to-mean envelope power ratio over all the code words 
in the selected cosets is minimized. 

20. The method of claim 18, wherein the code words are 
selected from cosets for which the peak-to-mean envelope 
power ratio is no greater than a predetermined value. 

21. The method of claim 20, wherein said code words are 
binary or quaternary code words and said predetermined 
value is 4, or said code words are octary code words and said 
predetermined value is 3. 

22. The method of claim 18, wherein the code words are 
selected from cosets whose representatives are defined by 
linear combinations of rows derived by symbol-wise mul- 
tiplication of all possible pairs of the first m rows of the 
generator matrix and are further multiplied by a positive 20 
integer s, where s is a power of two less than the number of 
different possible values for each symbol. 

23. The method of claim 1, including the steps of modu- 
lating multiple carrier signals in accordance with the code 
words encoding said data words, and combining the modu- 25 
lated carrier signals, to produce a composite signal having a 
peak-to-mean envelope power ratio reduced to that for a 
composite signal modulated directly in accordance with said 
data words. 

24. A method for transmitting data words comprising: 
selecting a set of one or more cosets of a linear sub-code 

of a code having a generator matrix as follows: 
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selecting a set of one or more complete cosets of a linear 
sub -code of a code having a generator matrix as fol- 
lows: 



((2.0+ l)x2 fl, )/2 m ((11 + l)x2* , )/2" 
((2.0+ l)x2 m " 1 )/2 m - 1 ((2.1 + [)x2 n - l )fZ m - 1 



((2.0+l)x2')/2* 

((2.0+ I)x2 , )/2 I 
((Z0+l)x2°)/2° 



((Zl + l)x2')/2J f 

((2.1 + l)x2 1 )/2 l 
(ai + l)x2°)/2° 



((2^ + l)x2 ffl )/2 M 
((U + l)x2 m - 1 )f2 m 

«2* + l)x2')/2* 

«2^+l)x2 l )/2 l 
((2.z + l)x2°)/2 0 



((2.(2 m -l)+l)x2 m )/2 ffl 
((2.(2" -1) + l)x2 M - l )/2'"- 

((2.(2 m -l) + l)x2J')/2-' 

U2.(2"-l) + l)x2 l )/2 1 
((2.(2"-l)+-l)x2°)/2° 



((2.0 + l)x 2W 



((Zl + l)x2*)/2' 



where z has integer values from 0 to 2 m -l, y has integer 
values from m to 0, x indicates bitwise multiplication 
of factors expressed as (m+l)-bit binary numbers, and 
the result of division is expressed as a single digit; 

encoding the data words as 2 m -symbol code words, 
wherein m is at least 4 and the data words are encoded 
in accordance with their value and code words selected 
from the set of one or more complete cosets; and 

transmitting the code words using multi-carrier frequency 
division multiplexing. 

26. A method of encoding data words comprising: 

selecting a set of one or more complete cosets of a code 
having a generator matrix as follows: 



((Z0 + l)x2 1 )/2 1 
((Z0 + l)x2°)/2° 



((2.1 + l)x2 1 )/2 1 
((2.1 + l)x2°)/2° 



((2.? + l)x2 m )/2 J " 
((2,^ + l)x2 w - t )/2 ff, - 1 



((2.(2'" -1)+ l)x2 m )/2 m 
((2.(2f* , -l) + l)x2r- | )/2™- 1 

((2.(2"-l)+l)x2 y )/2 y 
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((2.0+ l)x2 M )/2 w ' ((11 + l)x2T)/2 m 
((2.0+ l)x2 m - 1 )/2 m_1 ((2.1 + i)x2T' l )/2 m - L 



<(2.0+l)x2>)/2> 

((2.0+ l)x2')/2 1 
((2.0+l)x2°)/2° 



((ll + l)x2')/2' 

((Zl + Dx^)/^ 
((2.1 + l)x2°)/2° 



«2.z + l)x2 l )/2 l 
((2.* + l)x2°)/2° 



((2.(2 i "-l)+l)x2 , )/2 1 50 
((2.(2 w -l)+l)x2°)/2° 



where z has integer values from 0 to 2 m -l, y has integer 
values from m to 0, x indicates bitwise multiplication 
of factors expressed as (m+l)-bit binary numbers, and 
the result of division is expressed as a single digit; 

encoding each data word in accordance with its value and 
code words selected in complete cosets from the set of 
one or more cosets for desired peak-to-mean envelope 
power ratio characteristics of transmission; 

transmitting the encoded data words as 2™ -symbol code 
words using multi-carrier frequency division multi- 
plexing; and 

decoding each code word using an analytical decoder. 
25. A method of encoding and transmitting data words 
comprising: 
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((l2+l)x2 m - 1 )/2 m 

«2* + l)x2')/2" 

((2.2+ I)x2 1 )/2 I 
((2.z+l)x2°)/2° 



((2.(2"- 0+1) x2 m )(2 m 
((2.(2 w, -l)+l)x2 m - 1 )/2'"- 1 

((2.(2 m -l) + l)x2 J ')/2' 

((2.(r-l) + l)x2 1 )/2 1 
((2.(2 m -l) + l)x2°)/2° 



where z has integer values from 0 to 2 m -l, y has integer 
values from m to 0, x indicates bitwise multiplication 
of factors expressed as (m+l)-bit binary numbers, and 
the result of division is expressed as a single digit; and 

encoding the data words as 2 m -symbol code words for 
transfer over a communication path, each code word 
symbol having 2 / different possible values, wherein m 
is at least 4 and j is at least 3, and the data words are 



04/17/2004, EAST Version: 1.4.1 



US 6,373,859 Bl 



77 



78 



encoded in accordance with their value and code words 

selected from the set of one or more complete cosets. 
27. A method of encoding data comprising: 
selecting code words from a set of one or more cosets of 

a linear sub-code of a code having a generator matrix 

as follows: 



((2.0 + l)x2 m )/2 ,B 
((2.0 + l)x2 m ~ 1 )/2 m_l 



((2.1 + \)xT)/2 m 
((2.1 + l)x2'"- 1 )/2 m ~ 1 



((2.0 + DX2W ((2.1 + l)x 2^/2' 



((2.0 + l)x2 1 )/2 1 
((2.0 + l)x2°)/2° 



((2.1 + I)x2 l )f2 l 
((2.1 + l)x2°)/2° 



((2. Z + l)x2 ffl )/2' n 
((Zz + i)x2 w-1 )/2 w 

((2.1 + 0x2*)/ 2* 

((2. z + l)x2 1 )/2 1 
((Zz + l)x2°)/2° 



((2.(2 m -1)+ I)x2 m )/2 B1 
((2,(2" - l)-Hl)x2 w " 1 )/2 w_t 

((2.(2 M -l)+l)x2 y )/2 J ' 



10 



15 



20 



of factors expressed as (m+l)-bit binary numbers, and 
the result of division is expressed as a single digit, and 
each code word symbol has 2 1 different possible values 
where j>l, the method comprising: 

receiving a set of code words; 

deriving coefficients for rows of the matrix by applying a 
plurality of fast Hadamard transforms or equivalents 
thereof to a 2 m -symbol input word; and 

combining results of the transforms or their equivalents to 
identify the 2 m -symbol code word. 

31. The method of claim 30, wherein a total of 2 (/ ~ lXm+1) 
fast Hadamard transforms or equivalents thereof are applied 
in respect of said input word. 

32. The method of claim 31, wherein coefficients for rows 
of said matrix are derived by: 

generating 2 tf_1 ^ m+1 ^ values B according to the relation 



((2.(2" 
((2.(2" 



-1)+L)x2 l )/2 1 
-l) + l)x2°)/2° 



where z has integer values from 0 to 2 m -l, y has integer 
values from m to 0, x indicates bitwise multiplication 
of factors expressed as (m+l)-bit binary numbers, and 
the result of division is expressed as a single digit; and 

encoding the data words as 2 m -symbol code words for 
transfer over a communication path, each code word 
symbol having 2! different possible values, wherein j is 
a variable positive integer, and the data words are 
encoded in accordance with their value and a code 
having, for any value of j, the selected code words. 

28. The method of claim 1, wherein said set of cosets 
contains a plurality of cosets and said linear sub-code is the 
same for each coset. 

29. The method of claim 28, wherein said linear sub-code 
comprises all of said code. 

30. A method of identifying a 2 m -symbol code word 
within a set of code words having a generator matrix 



for all possible combinations of values of coefficients b 1( 
b 2) . . . , b m and b', each of which can have integer values 
25 between 0 and 2 / ~ 1 -l inclusive, where X, to X m are row 
labels of the first m rows in the matrix, X' is a row label of 
the final row of said matrix and 0 indicates symbol-wise 
addition modulo 2 J ; 

30 generating for each value B a vector D according to the 
relation 
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((2.0 + l)x2 m )/2 M 



((2.l4.1)x2" , )/2" 



((2.0 ■t-l)x2 m " 1 )/2 m - 1 ((2.1 + l)x2 w_1 )/2 m 



(aO + l)x2J')/2J r 

((2.0 + l)x2 1 )/2 1 
((2.0 + l)x2°)/2° 



((2.H-l)x2 J ')/2J r 



((2.1 
((2.1 



-ljx^)/^ 
-l)x2°)/2° 
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((2.z + l)x2 m )/2" 
((2.z + l)x2 w_l J/2' w 

((2. Z + l)x2^)/2 J ' 

((2.z + l)x2 1 )/2 1 
((2.z + l)x2°)/2° 



((2.(2' 
02.(2" - 



-1)+ l)x2 m )/2 m 
■i)^i)x2T l )/2 m - i 



D^lT^X-^S-B) mod l^-l^X' 

where S is said input word and subtractions involving values 
X', S and B, and the modulus operation | ... | are performed 
symbol-wise; 

applying a fast Hadamard transform or an equivalent 
thereof to each vector D to derive 2 m transform com- 
ponents; 

identifying a transform component of largest absolute 
value; and 

deriving row coefficients corresponding to the value of B 
and the result of the fast Hadamard transform or 
equivalent thereof associated with said identified trans- 
form component. 

33. The method of claim 30, wherein said symbols of said 
input word have integer values. 

34. A coder for encoding data words comprising: 

first circuitry for selecting code words in complete cosets 
from a set of one or more cosets of a linear sub-code of 
a code having a generator matrix as follows: 



((2.(2"- 1) + Ijx2 y )/2 J ' 



((2.0+ l)x2 m )/2 m 
((2,0+ l)x2 m - 1 )/2 m - 



((2.1 + I)x2")/2 m 
((2,l + l)x2* , - 1 )/2 n 
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((2.(2" 
((2.(2" 



-l)x2°}/2° 



((2.0+ l)x 2^/2* ((2.1 + l)x2 J ')/2^ 
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((2.0 ■ 
((2.0 ■ 



-l)x2 1 )/2 1 
- l)x2°)/2° 



{{2S + i)x2 l )/2 l 
((2.1 + l)x2°)/2 0 



where z has integer values from 0 to 2 m -l, y has integer 
values from m to 0, x indicates bitwise multiplication 
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-continued 



((2.n-l)x2 m )/2 ,n 
«2.z + l)x2°)/2° 



((2.(2" , -l) + l)x2 m )/2 m 
((2.(2 V7, -l) + l)x2 OT - l )/2 m - 1 

((2.(2"-l)+l)x2>')/2> 

((2.(2*'-l)+l)x2 l )/2 l 
((2.(2"-l)+l)x2°)/2° 
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where z has integer values from 0 to 2 m -l, y has integer 
values from m to 0, x indicates bitwise multiplication 
of factors expressed as (m+l)-bit binary numbers, and 
the result of division is expressed as a single digit; and 

second circuitry for receiving data words and for encod- 
ing them as 2 m -symbol code words in accordance with 
their value and the code words selected in complete 
cosets, to obtain desired peak- to -mean envelope power 
ratio characteristics of transmission using multi-carrier 
frequency division multiplexing. 
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